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Hamilton Harbour Study 197^ > J^ne 1975 

INTRODUCTION 

In 1972-73, the Ontario Ministry of the Environment under- 
took a water quality study of Hamilton Harbour. The study 

obJeGtlves were: to determine existing water quality and 
Its trends over recent years; to Identify water quality 
problems and recommend corrective measures;- and to determine 
the effects of landfllllng on water quality. One of .the 
recommendations of that study was that artificial destratl- 
fication be considered as a means of Improving harbour water 
quality. Artificial" mixing will commence In June 1975 • 

The purpose of this report is to provide further information 
on the water quality, including 197 4 data, to provide a 
literature survey and comment on other artificial mixing 
experiments and to indicate the expected results of the 
artificial destratlfication experiment in Hamilton Harbour. 
The report format is similar to that of the 197 4 report with 
the subject matter broken down into components; each is 
presented in technical detail in appendices and integrated 
in the following text. The sections are as follows: 

A. Water Quality Update - 197^. 

B. Sediment Statistical Analysis. 
C'. Dispersion. 

D. Periodicities and Mass Exchange through Burlington Ship 
Canal . 

E. Sediment Dissolved Oxygen Uptake. 
P. Phytoplankton. 

G. Artificial Mixing and Water Quality in Hamilton Harbour. 

As this is essentially a follow-up to the report released in 

Mayt 197^9 certain relevant sections of that report have 
been reprinted in Section H for ready reference. In Section 
H, one finds the description of the harbour,, purpose of- the 
previous study ^ study outline ^ summary and findings ^ recomm.--ic 
atlons and further studies. Section I is a summary of the 
proposed water quality surveys for 1975-7. 

DISCUSSIO'N 

A. Water Quality Update - 197^' 

The water quality survey program, in 1974 consisted of collect- 
ing surface and bottom, water samples four times at six 

different locations as well as operating six recording 
chemistry and current meters for a period of several weeks. 

While these surveys were limited ^ some obeervatlons can be 
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made. Armrionia (nitrogen) levels continued to decrease from 
1972 to 197^ and are now generally within the Ministry of 
the Environment public surface water supply criteria of 0.5 
mg/l as N. Total phosphorus did not change significantly 
from 1972 to 1974, Although soluble phosphorus and chlorophyll 
a concentrations in July and October appear lower in 1974 
than 1972 5 the available data were not adequate to quantify 
this observation. Analysis of the data from the recording 
chemistry meters confirmed the extent of the diseolved 
oxygen depletion as well as its variation. The calibration 
of the recording chemistry meters was checked with in situ 
volume samples. Severe dissolved oxygen depletion was 
evident in the hypolimnlon of an offshore location (20im 
depth) in the eastern end of the harbour. During a strongly 
stratified period, July 11 to July 21, 197^, daily mean 
values of dissolved oxygen were less than 2.5 rng/1 with a 
minimum value of 0.6 mg/l measured on July 13. The 1972 
survey showed the dissolved oxygen depletion to vary widely 
in the Windermere basin between successive weeks. This 
variation is confirmed by the recording chemistry meter 
results which showed the standard deviation of dissolved 
oxygen to be 20 percent of the mean value in the Windermere 
basin. This variation appears to be a function of the 
Hamilton Pollution Control Plant operation combined with 
variable natural assimilative characteristics of the harbour. 
Counts for total conform and enterococci are generally 
lower in 197^ than in 1972 but the sample size is too small 
to assess the change statistically. Fecal coliform counts 
were about the same In 197^ as in 1972. 

B. Sediment Statistical Analysis 

In the 1974 report 5 It was stated that the heavy metal 
discharges should be identified and action taken to reduce 
these discharges. A statistical correlation matrix method 
oT relating bottom sediments by the ratios' of trace metal 
and PCB. content has been developed and applied to the data 
from the 1972 surveys. Details of the method , including the 
computer program which has. some general applications (e.g. 
analyses which involve diversities) are contained in Section 
B. Results of the analyses indicate that generic relation- 
ships exist for the sediment In the eastern portion of the 
harbour (adjacent to the B^urllngton Skyway) and the northeast 
shore, ThiS' suggests that sediment chemistry is related to 
the same sO'Urcej probably the Hamilton PollutlO'n Control 
Plant outfall. 

C.. Dispersion 

The dispersion coefficients used in the numerical model were 
transposed from other lake experiments, Co^nsequently , a 
series of limited multiple drogue tracking experiments were 
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carried out In 197^ to validate the range of dispersion 
coefficients used. As a function of time, the measured 
dispersion coef fie lent behaved in a manner similar to that 
In other Great Lakes studies; the value of 1 . 2 x 10** cm^/ce; 
used In the numerical model was shown to be satisfactory. 

1* Hamilton Harbour Periodicities and Mass Exchange Through 
the Burlington Ship Canal „^_^ _^_^_ 

In the previous report the harbour periodicities In both 
movement and water chemistry were found to be of prime 
importance In the Interpretation of trends in water quality j, 
process analysis and the application of a numerical model. 
In 1974 J six recording current and chemistry meters were 
operated In the harbour to confirm this observation. The 
periodicities discussed In, the previous report were observed 
In 1974 in both the eplllmnlon and hypollmnlon at the six 
locations surveyed in the eastern portion of the harbour ana 
canal. Periodicities related to lake and harbour oscilla- 
tions 3, semi-diurnal processes and the harbour oscillations 
related to the small canal connection to the lake (so-callea 
Helm_holtz mDdes) were Identified in the water chemistry and 
mDvement records.- These periodicities ranged from, 12 hours 
to 20 minutes and Influenced the Interpretation of the water 
quality data obtained by conventional boat survey techniques. 

The 1974 report indicated the importance of the exchange 
process between the harbour and the lake in maintaining 

water quality In the harbour. Analysis of the data collected 
in 1974 In the Burlington Ship Canal indicated that mass 
exchanges of the order of one percent of the harbour volume 
can occur In one day and persist for p^'rlods in excess of 20 
days. This l;s a very large exchange rate compared to the 
throughput computation based on water balance which indi- 
cated an exchange rate of only 80' percent of the harbour 
volume in one year (mean value of 0.2 percent of the harbour 
volume per day). The Importance of these exchange me'Lhanlsm,L 
In maintaining water quality in the harbour Is obvious. It 
is Important that any modifications to the existing Barlingtc 
Ship Canal or developments on the lake side of the canal be 
carried out in such a way that water quality be maintained 
In the harbour. It is desirable that any works In the 
vicinity of the Burlington Ship canal be preceded by an 
environmental Impact statement which demonstrates that such 
work will not reduce the exchange mechanisms with the lake. 

There is an apparent conflict concerning the Importance of 
the canal exchange mechanism in maintaining water quality in 
the harbour. The results of the nume'rlcal mO'del run in the 

1974 report showed that Increasing these exchange rates did 
not significantly improve the harbour water quality. The 
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Increased exchange rate between the harbour and the lake 
used in the numerical model was approximately 0.03 percent 
of the harbour volume per day. This value was selected as 
being equivalent to the total industrial discharges to the 
harbour. The Influence of this email exchange rate increase 
on the water quality could not be detected with the sensi- 
tivity of the numerical model. However^ if a value of one 
percent of the harbour volume per day or 8lO mVsec was used 
in the numerical model, pronounced changes would occur in 
the harbour water quality prediction. Artificially pumping 
such a large quantity of water is nevertheless not a practi- 
cal consideration Tor enhancing harbour water quality. 

£. Sediment Dissolved Oxygen Uptake 

Low dissolved oxygen levels in the hypolimnion were identi- 
fied as a major water quality concern and artificial mixing 
was proposed as a means of increasing dissolved oxygen 
stocks hy natural aeration in the 197^ report. in this 
reporr, it was estimated that bottom sediments and discharges 
accounted for approximately two-thirds of the dissolved 
oxygen deficit during the critical period in August- While 
dissolved oxygen demand of the principal discharges to the 
iiarbour are known 5 sediment demands v;ere unknown. No 
measurements of sediment dissolved oxygen demand in Hamilton 
Harbour were undertaken since the 1972 surveys. An extensive 
literature survey of the sediment dissolved oxygen uptake 
was undertaken to determine the best manner of measuring 
Chese demands and to permit a better estimate of these 
demands in Hamilton Harbour. This section discusses different 
methods for assessing sediment dissolved oxygen demand and 
recomputes the dissolved oxygen demand of sediments using 
transposed data and the dissolved oxygen budget. Generally 
it is observed that in situ measurements of sediment dissolved 
oxygen uptake must be made before reliable estimates of 
dissolved oxygen demand can be made. Such measurements are 
expensive as" the variability of both the sediments throughout 
the harbour and the supply mechanisms must be represented. 

An estimate of sediment oxygen uptake in Hamilton Harbour for 
a period (April 28 - August 19) in 1972 was obtained by using 
measured results obtained from other studies performed under 
similar conditions. It was estimated that 3 x 10^ kg of 
dissolved oxygen was taken up by the sediments in this period, 
as compared with an estimated deficit of 1.5x10^ kg (computed 
from observed and saturation stocks in August 1972), Of this 
deficit, 0.5 X 10^ kg of dissolved oxygen was assumed due to 
seston decay. Estimated values of dissolved oxygen input to 
the harbour" from the mass exchange throughout the canal for 
the same period ranged from 0.4 to 4.3 x 10^ kg, depending 
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upon conditions of flow and stratification in the canal. 
The importance of the harbour-lake exchange "as an oxygen 
source is apparent when considering stO'Cks. This suggests 
that the estimate in the 197^ report that effluent and 
sediment demand represent apprO'Xlmately two-thirds of the 
deficit is not unreasonable.. In May 1972s dally industi-ial 
and municipal biological and chemical o^xygen demand loadings 
to the harbour were three times the estimated sediment 
oxygen demand value. However , all the discharges are to the 
surface waters with much of the oxygen demand satisfied by 
aeration 5 photosynthesis and lake water. Only a portion of 
the oxygen demanding discharges are particulate matter 
which will be dep-osited in the sediments. 

F. Phytoplankton 

In 1972, measurements were made of chlorophyll a_ as an 
Indication of phytoplankton standing crops. A better 
definition of the phytoplankton community was required to 
predict the effects of artificial (^estratlflcation. .in 
1974 J enumerations of phytoplankton were made at six loca- 
tions on six different occasions. The samples were collected 
by passing a 32 oz. bottle through the euphotle zone. Each 
sample was identified to genera and standing crop expressed 
in areal standard units. Unfortunately, the survey was too 
limited in extent to determine spatial distribution of 
phytoplankton. There was a conspicuous lack of blue-green 
algae (Cyanophyceae ) although the pH (7.0-8.0) and moderate 
alkalinity (80-100 mg/1 as CaCOa) would provide an abundance 
of carbon dioxide to support both green and blue-green 
algae. Genera associated with organically polluted waters 
such as Euglena , Phacus , Cryptomonas and Rhodomonas pre- 
dominated. Production is relatively low compared to the 
levels in Cootes Paradise (connected to the harbour at the 
western end). The available nutrients In the harbour are 
sufficient to support large algal populations. The low 
production In 1974 is likely a result of the lack of staMl- 
ity in the vjater column due to natural vertical mlxiU g 
and/or some unidentified toxicity. 

Limited zooplankton data show moderate levels which are 
unlikely to affect phytoplankton levels. One interesting 
observation was the absence of Daphnla which is co^mmon to 
mO'St fresh water co^mmunitles , 

G. Artificial Mixing and Water Quality in Hamilton Harbour ' 

Artificial mixing has been used many time^s to Improve water 
quality and manage biO'loglcal productivity in aquatic systems, 
This section presents a literature survey of artificial 
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mixing experiments, discusses the results and indicates the 
expected results In Hamilton HarbO'Ur. 

The general water quality ifflprO'Vements associated with 

artificial mixing are: 

1. Elimination of taste and odour 

2. Improved transparency 

3. Larger habitat for aquatic life ^ 

4. Increased dissolved oxygen ** 

5. Reduced algae blooms 

6. More neutral pH as carbon dioxide Is released from 

hypolimnlon '"^! 

7. Decreased dissolved metal content particularly iro^n 
and manganese 

The success of artificial destratif Icatlon is dependent on **' 

the establishment of an oxidized mlcrozone at the sediment- 
water interface. 

The main goal of the destratif icatlon in Hamilton Harbour Is 

to relieve anoxic conditions in the hypolimnlon. If the 

anoxic condition is the factor limiting the establishment of ^ 

aquatic life a more extensive and diverse habitat should 

develop capable of supporting an improved fish community. 

It is expected that primary production will not increase 

significantly as it will be limited by turbidity. It is -^• 

also not expected that algal composition will change to 

blue-greens as control of pH, alkalinity and free carbon 

dioxide will not favour their development. General water « 

quality should also improve although it is not possible to 

predict the extent as the chemical and biological responses 

to artificial destratif icatlon are not known. ^ 
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SUMMARY 

Water Quality 

Wltli the exception of the ammonia levels, there was no 
general Improvement In water quality in the harbour between 
1972 and 197^. Dissolved oxygen levels in the hypolimnion 
appear to have declined since 1972 indicating the need for 
corrective action. The importance of dissolved oxygen 
supplied by the exchange between the harbour and lake was 
identified by estimation of sediment oxygen uptake 3 as well 
as water chemistry and current measurements. Phytoplankton 
densities are generally low for the level of nutrients 
available. Limited production could be caused by physical 
processes preventing the establishment of a stabilized water 
column for extended periods of time. The predominant 
species are members of the Chlorophyceae . Peculiarities in 
zooplankton composition might be related to low dissolved 
oxygen levels. 

Sediment AnalysliS 

A statistical technique for relating the heavy metal contents 
of sediment was applied to the survey data of 1972. It was 
found that heavy me'tal ratios of the sediments in the 
eastern and northeastern portion of the harbour were related 
to a similar source, probably the Hamilton sewage treatment 
plant effluent- The sediment distribution was related to 
the predominant circulation patterns in the harbour. 

Artificial Mixing 

A compilation of data obtained in other artificial mixing 
experiments and comparison with conditions in HamlltO'n 
Harbour indicate that such a procedure should eliminate the 
dissolved oxygen deficits in the hypolimnion. None of the 
other experiments we're conducted on such a large body of 
■water receiving large Industrial discharges. Similarly , it 
is not known what effect artificial mixing will have on the 
exchange mechanisms between the harbour and lake which in 
the past have been Important in providing dissolved oxygen 
to the harbour stocks. 

Plankton surveys Indicate that, the communities are diverse 
enough to support fish provided the vertical habitat can be 
extended. Artificial mixing should extend this habitat if 
the limitation to present production Is not related to some 
factor other than low dissolved oxygen levels. It is not 
expected that the present phytoplankton species (greens) 
will increase greatly in, density. Blue-green species dens- 
ities will be limited by the disruption of stratification. 
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RECOMMENDATIONS 

Water quality surveys should be conducted in support of the 
artificial mixing to determine the effect on plankton 

communities and productivity ^ fish abundance and distri- 
bution as well as harbour water quality. 

Further sediment analysis should be undertaken to identify 
sO'Urces O'f contamination. 

Proposed works In the region of the ship eanal must be 
assessed in terms of their Influence on the exchange mechanisms 

between the harbour and the lake. 
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HAMILTON HARBOUR STUDY 

SECTION A 
WATER QUALITY UPDATE 

SUMMARY 

RESULTS OF WATER QUALITY SURVEYS CONDUCTED DURING 1974 IN 
HAMILTON HARBOUR ARE PRESENTED. FOR MOST PARAMETERS/ INSUF- 
FICIENT DATA WERE COLLECTED TO DRAW FIRM CONCLUSIONS. 

NUTRIENT DATA INDICATED DECREASES FROM 1972 TO 1974 IN 
AMMONIA AND TOTAL KJELDAHL NITROGEN/ AND AN INCREASE IN 
NITRATE DURING THE SAME PERIOD. 

DISSOLVED OXYGEN DATA COLLECTED BY AUTOMATIC DIGITALLY 
RECORDING WATER QUALITY METERS INDICATED THAT A SEVERE 
DISSOLVED OXYGEN DEPLETION PROBLEM CONTINUES TO EXIST 
WITHIN THE HYPOLIMNION. LARGE VARIATIONS IN DISSOLVED 
OXYGEN CONTENT/ WHICH OCCUR WITHIN HOURS NEAR THE HAMILTON 
SEWAGE TREATMENT PLANT/ INDICATE THE VARIABILITY OF LOAD- 
INGS AND DYNAMICS IN THE HARBOUR. 

BACTERIAL DATA IN 1974 ARE IN GENERAL NOT SIGNIFICANTLY 
DIFFERENT FROM RESULTS OBTAINED IN 1970-72. 



HAMILTON HARBOUR STUDY 

WATER QUALITY UPDATE 

To assess the continuing changes in water quality in Hamilton Harbour 
and the effects of reduced nutrient loadings, a water quality survey was 
conducted in 1974 by the Ontario Ministry of the Environment. Surface 
and bottom samples were collected at the six locations indicated in Figure 
1 . In addition, automatic recording water quality meter records of temp- 
erature, dissolved oxygen, pH and conductivity were obtained in June 
and July near the Windermere Basin (location M8) , at the mouth of the 
Ottawa Street slip (location M5) , and in the epilimnion and hypollmnion 
of an offshore location at the eastern end of the harbour about 1 km south 
of the Burlington ship canal (locations M6 and M7) . The results of the 
latter are discussed more fully in Section D, Hamilton Harbour Periodic- 
ities and Mass Exchange. 

Unfortunately, the six station set was sampled on only four occasions 
during 1974, and only two depths were obtained at each station. As a 
result detailed statistical analysis of the data is impossible. Enough 
data were collected to provide a rough estimate of continuing trends 
in nutrients, dissolved oxygen and microbiology. These will be dis- 
cussed briefly in turn. 

Nutrients 

Harbour-wide means for nutrients and chlorophyll a are given in Table 1. 
It is apparent that a further large decrease in ammonia nitrogen has occur- 
red from 1972 to 1974, with 1974 values being largely within the M.O.E. 
public surface water criterion of 0. 5 mg/l . Accompanying this is a de- 
creased total Kjeldahl nitrogen value. This provides continuing evidence 
of the decreased industrial loadings previously stated (M.O.E., 1974), and 
particularly the diversion by Steico of ammonia-bearing wastes to the Hamilton 
W.P.C.P. early in 1974. Increased bacterial nitrification may also be important, 
as the results indicate an increase in nitrate from 1972 to 1974. The 1974 
figures are consistent with the previously demonstrated occurrence of nitrifi- 
cation during early summer followed by nitrate reduction (M.O.E., 1974), 
but no definite conclusion can be drawn from the small data set. 
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Total phosphorus did not change significantly from 1972 to 1974. An apparent 
decrease in soluble phosphorus is evident, except possibly in Ocotober. 
However, the samples for soluble phosphorus were not frozen in the field, 
as is desirable for proper preservation. Uptake of dissolved phosphorus 
by microplankton, or adsorption by particulate matter, could have occurred 
between sampling and analysis. IJC studies indicated that confidence limits 
were about 0.003 mg/l, comparable to the July-August 1974 means. This Is 
comparable to deviations recorded here. 

The July and October chlorophyll a results are considerably lower In 1974 
than in 1972. It is difficult to draw a conclusion regarding this parameter, 
as its seasonal variation Is known to be large (individual weekly means in 
July-October 1972 varied from 11 to 28 ug/l), and a period of high phyto- 
plankton growth may easily have occurred at some other date in 1974. 

Presently secondary treatment is reducing effluent total phosphorus levels 
from the Hamilton water pollution control plant to less than 1 .5 mg/l total 
phosphorus. By the end of 1975 the plant will be meeting provincial objectives 
of 1 mg/l . 

During 1975, weekly sampling at 4 stations and monthly sampling at 8 stations 
is to be done to further confirm the existence of seasonal and year-to-year 
trends. 

Dissolved Oxygen 

Dissolved oxygen was measured at the six station grid on only two dates in 
1974. Values obtained are compared with 1972 data in Table 2. Although 
no firm conclusion can be drawn from such meager data, they suggest that 
the DO depletion in 1974 was as great as, if not greater than, that observed 
in 1972. 

Definite evidence of the magnitude of the oxygen depletion problem in 1974 
was provided by the results of the automatic digitally-recording water chem- 
istry meters deployed in June to July at locations M6, M7 and M8. Daily means 
and standard deviations of temperature, dissolved oxygen, pH and conductivity 
are indicated in Tables 3 to 5, and a portion of the record obtained in the off- 
shore location hypolimnlon (M7) is shown In Figure 2. That stratification 
existed during the study is evidenced by the temperature means in Tables 3 
and 4, as well as the temperature profile recorded on July 16, 1974 at 1215 hrs, 
(Figure 3) . In the hypolimnlon, all daily mean dissolved oxygen values were 
below 2.5 mg/l, with the minimum daily value recorded on July 13, when the 
mean was 0.63 mg/l and the maximum reached 0.96 mg/l. This indicates that 
the hypolimnlon of Hamilton Harbour was seriously depleted of oxygen during 
the entire recorded period. It is hoped that the proposed mixing program in 
1975 will alleviate most of this depletion. 
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The extent of DO depletion at station 20 as recorded in 1972 varied tremend- 
ously between successive weekly surveys. Surface values in the July to 
October 1972 period varied from a minimum of 0.6 mg/l on October 24 to 
a maximum of 9.6 mg/l on August 6. That large between-day variations of 
dissolved oxygen do occur is confirmed by the data in Table 5, which were 
checked and calibrated with in situ field calibration checks. This indicates 
the variable effect of the Hamilton W.P.C.P. discharge, combined with bio- 
logical mechanisms of photosynthesis, respiration and sediment oxygen 
demand, and atmospheric reaeration, on the dissolved oxygen on this portion 
of the harbour. Because of the remote location relative to the proposed location 
of the aerator, and shallow depth of this area, combined with the effect of the 
W.P.C.P. outfall, it is not known if the proposed mixing program will have 
much effect on dissolved oxygen at this location. 

Microbiology 

The sampling locations are shown in Figure 1, and the bacteriological water 
quality data for the 1974 surveys of the Ontario Ministry of the Environment 
appear in Table 6 (Total Coliform), Table 7 (Fecal Coliform) and Table 8 
(Enterococci) . 

Table 9 indicates the geometric means of fecal conforms per 100 ml. for the 
years 1970, 1971 , 1972 and 1974. Because of the low number of samples taken 
in 1974, means are being compared with single samples. The year 1971 still 
appears as the year with the lowest counts. The significance of the differences 
between the years was tested statistically with a t-test for some sampling 
locations with an adequate number of determinations and the results appear 
in Table 10. The differences between locations are not significantly differ- 
ent at the 95% significance level; consequently, there are not significant 
statistical differences between 1970 and 1971, 1971 and 1972, 1972 and 1974, 
and 1970 and 1974. 

By comparison of the related tables of 1972 and 1974, there is an improvement 
in total coliform and enterococci counts. Fecal coliform counts remained about 
the same. However, many of the 1974 values are single samples providing at 
best crude qualitative indicators. 
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NH3 (N) 

Total Kjeldahl 
Nitrogen (N) 

Nitrate (N) 

Total 
Phosphorus (P) 

Soluble 
Phosphorus (P) 

Chlorophyll a 



TABLE 1 
NUTRIENT SURFACE MEANS, 1972-4 







JULY 








AUGUST 








OCTOBER 




1972 




1974 




1972 




1974 




1972 




1974 










Number 






Number 






Number 


Mean 


Mean 


S. D. c 


>f Samples 


Mean 


Mean 


S. D. c 


>f Samples 


Mean 


Mean 


S. D. c 


>f Samples 


2.4 


0.016 


0.026 


6 


1.9 


0.30 


0.14 


5 


1.5 


0.42 


0.35 


6 


3.4 


2.13 


0.78 


6 


3.0 


0.92 


0.30 


11 


2.2 


1.30 


0.71 


6 


1.89 


2.6 


1.8 


6 


1.77 


3.5 


0.4 


6 


1.58 


2.5* 


0.7 


6 



0.087 0.067 0.017 

0.016 0.002 0.002 
23 7.8 2.3 



0.082 0.102 0.048 11 

0.024 0.004 0.006 6 
20 19.3 2.7 6 



0.093 0.057 0.027 

0.021 0.017 0.014 
18 4.6 0.9 



NOTE : All concentrations are in mg/1, except chlorophyll a in Mg/1. 

* This result is (NO3 + NO2) - N, NO2 - N results were unreliably 
high at several points; NO3 - N is affected as this is obtained 
by difference. 
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TABLE 2 



Date (1974) 



August 21 (surface) 
(bottom) 

October 2 (surface) 
(bottom) 



DISSOLVED OXYGEN DATA, 1972-74 
Station Mean D.O 



20 


253 


254 


256 


262 


265 


7.5 


_ 


8.7 


7.8 


7.2 


9.8 


3.3 


- 


0.0 


3.0 


5.2 


4.5 


6.2 


6.4 


6.9 


5.7 


7.5 


- 


5.6 


4.4 


4.5 


7.2 


7.2 


- 



mg/l 

8.2 
3.2 

6.5 
5.8 



% Sat, 

96 

34 

62 

55 



1972 Mean D.O 



mg/l 

7.1 
3.7 

7.3 

5.4 



% Sat. 

84 
41 

77 
55 






Note : Individual figures are in mg/l. 1972 data are 

8 - station means for August 29 and October 4-5. 



TABLE 3 
WATER QUALTIY METER OPERATED AT OFFSHORE LOCATION EPILIMNION (M6) 

D/>ILY MEANS STARTING AT 16.25 HR ON GIVEN DATE 
(HALF - HOURLY READINGS) 



DATE 


V 1 1 / »•_■ 


TEMP(C) 


D.O.(MG/Li 


0.0. i%) 
SATURATION) 


PH 


C0N0(25 


C) 

1 


74192. 


MEAN 
ST DEV 


21.93 
0.61 


6.97 

1.04 


79.14 
12.61 


7.84 

0.12 


527.48 
7.91 




74193. 


MEAN 
ST DEV 


23.18 
0.57 


9.66 
1.15 


112.32 
14.35 


8.27 
0.20 


517. <i7 
10.87 




74194. 


MEAN 
ST OEV 


23.65 
0.22 


9.58 
0.80 


112.29 
9.35 


8.31 
0.16 


515.40 
5.66 


{ 


74195. 


MEAN 
ST DEV 


24.23 

0.86 


7.84 
0.69 


92.9^ 

9.09 


7.75 

0.40 


521.33 
4.42 


1 


74196. 


MEAN 
ST OEV 


22.18 
0.46 


6.96 
0.58 


79.35 
6.81 


7.88 
0.19 


518.42 

10.71 


I 


74197. 


MEAN 
ST DEV 


21.45 
0.84 


6.37 
2.14 


71.93 
24.97 


7.R1 
0.29 


520.52 
5.58 


■■^ 


74198. 


MEAN 
ST OEV 


22.47 
0.13 


8.09 
I . 1 1 


92.68 
12.63 


7.95 
0.26 


517.65 
5.00 


•"^ 


74199. 


MEAN 
ST DEV 


22.55 
0.22 


7.57 
0.33 


86.92 
3.87 


7.86 
0.15 


517.35 
4.38 


-^ 


74200. 


MEAN 
ST DEV 


21.98 
0.55 


7.11 
0.65 


80.69 
7.32 


7.67 
0.28 


518.40 
9.39 


.-« 


74201. 


MEAN 
ST DEV 


20.51 

0.54 


6.16 

0.88 


68.03 
10.36 


7.69 

0.21 


511.18 
9.74 


-« 




GRAND MEAN 
ST DEV 


22.47 

1.14 


7.68 
1.54 


88.26 

18.87 


7.91 
0.32 


518.76 

e.36 


^ 



Note: Above dates are stated in terms of 365-day year. Day 74192 is 
July 11, 1974, etc. 
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TABLE 4 
WATER QUALTIY METER OPERATED AT OFFSHORE LOCATION EPILIMNION (M7) 

DAILY MEANS STAWTING AT 16,25 HR ON GIVEN DATE 
(HALF - HOURLY READINGS) 



DATE 




TEMP(C) 


D.0.(MG/L1 


D.O. {% 
SATURATION) 


PH 


CONO(25 C 


74192. 


MEAN 


12.87 


0.86 


8.08 


7.43 


471.23 




ST OEV 


0.54 


0.49 


4.54 


0.10 


9.33 


74193. 


MEAN 


12.82 


0.67 


6.35 


7.41 


471.31 




ST OEV 


0.30 


0.35 


3.24 


0.09 


7. 18 


74194. 


MEAN 


12.99 


0.63 


5.91 


7.40 


475. S4 




ST OEV 


0.23 


0.42 


3.95 


0.10 


5.38 


74195. 


MEAN 


12.55 


2.00 


18.63 


7.45 


461. 19 




ST OEV 


0.64 


1.09 


10.01 


0.11 


11.88 


74196. 


MEAN 


12.35 


1.17 


10.87 


7.45 


460.13 




ST OEV 


0.34 


0.79 


7.28 


0.11 


7.18 


74197. 


MEAN 


12.50 


1.27 


11.81 


7.45 


457.79 




ST DEV 


0.48 


0.72 


6.65 


0.11 


6. 18 


74198. 


MEAN 


12.17 


1.21 


11.25 


7.41 


454.50 




ST OEV 


0.36 


0.63 


5.80 


0.11 


6.70 


74199. 


MEAN 


12.84 


1.02 


9.56 


7.43 


463.06 




ST DEV 


1.11 


0.55 


5.07 


0.11 


14.19 


74200. 


MEAN 


12.52 


1.56 


14.61 


7.44 


450.94 




ST DEV 


0.51 


0.71 


6.57 


0.10 


8.40 


74201. 


MEAN 


12.84 


3.01 


28.43 


7.51 


447.19 




ST OEV 


0.47 


1.09 


10.41 


0.11 


5.72 




GRAND MEAN 


12.64 


1.28 


12.02 


7.44 


461.81 




ST DEV 


0.61 


0.93 


8.72 


0.11 


12.09 



Note: Above dates are stated in terms of 365-day year. Day 74192 is 
July 11, 1974, etc. 
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TABLE 5 
WATER QUALITY METER OPERATED NEAR WINDERMERE BASIN (M8) 

DAILY MEANS STARTING AT 18.11 HR QN GIVFN DATE 
(10 - MINUTE READINGS) 





TEMP(C) 


D.n.(MG/L) 


D.O. (% 
SATURATION) 


PH 


C0ND(25 C) 


MEAN 


18.58 


2.49 


26.37 


0.0 


575.65 


ST DEV 


0.55 


0.83 


8.87 


0.0 


19.69 


MEAN 


18.88 


2.79 


29.67 


0.0 


581.33 


ST DEV 


0.46 


1.07 


11.44 


0.0 


13.45 


MEAN 


19.17 


2.88 


31.01 


0.0 


591.62 


ST DEV 


0.41 


1.58 


17.24 


0.0 


rw 3 I 


MEAN 


19.50 


2.51 


27.14 


0.0 


601 .77 


ST DEV 


0.42 


1.17 


12.79 


0.0 


16.53 


MEAN 


21.10 


7.23 


80.23 


0.0 


608. Ofl 


ST DEV 


0.25 


0.84 


9.43 


0.0 


34.4 3 


MEAN 


21.32 


6.41 


71.58 


0.0 


654.80 


ST DEV 


0.42 


0.94 


10.64 


0.0 


38.99 


MEAN 


22.02 


6.09 


69.17 


0.0 


601.58 


ST DEV 


0.25 


l.Ol 


11.51 


0.0 


10.70 


MEAN 


21.83 


6.10 


69.30 


0.0 


604.34 


ST DEV 


0.19 


1.01 


11.58 


0.0 


17.87 


MEAN 


21.09 


4.05 


44.90 


0.0 


586.44 


ST DEV 


0.37 


0.69 


7.71 


0.0 


25.68 


MEAN 


22.02 


4.61 


52.29 


0.0 


58J. 10 


ST DEV 


0.59 


0.64 


7.61 


0.0 


26.40 


MEAN 


22.02 


3.16 


36.09 


0.0 


584.44 


ST DEV 


1.04 


1.64 


19.04 


0.0 


20.08 


MEAN 


20.85 


2.23 


25.10 


0.0 


543.81 


ST DEV 


1.26 


1.38 


15.94 


0.0 


22.02 


MEAN 


23.16 


4. 88 


56.97 


0.0 


663.89 


ST DEV 


0.41 


0.73 


9.28 


0.0 


56.27 


GRAND MEAN 


20.83 


4.2 5 


47.52 


0.0 


597.71 


ST DEV 


1.49 


2.01 


23.02 


0.0 


30.06 



>te: Above dates are stated in terms of 365-day year. Day 74176 is 
June 25, 1974, etc. 

Meter was not equipped with pH sensor. 
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TABLES: Hamilton Harbour Survey 1974 Total Coliform/100 ml 



Depth in Meters From Surface 

Station 

Number ]_J TO IJ 16 19 

20 240 775(2) 

253 190(2) 446(2) 

254 192(2) 



256 

262 1470 1800 

265 1689(2) 



580 



When more than one sample was taken at any location, the geometric mean 
is reported and the number of samples is recorded in parenthesis. 
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TABLE 7: Hamilton Harbour Survey 1974 Fecal Coliform /100 ml 

Depth in Meters From Surface 

t3 19 



Station 






Number 


1.5 


to 


20 


260 


243(2) 


253 


63(2) 


48 


254 


114(2) 


452 


256 


500 




262 


364 




265 


18(2) 





70(2) 



121(2) 



When more than one sample was taken at any location, the geometric mean 
is reported and the number of samples is recorded in paranthesis. 
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TABLE 8: Hamilton Harbour Survey 1974 Enterococci/100 ml 

Depth In Meters From Surface 

13 16 19 



Station 






Number 


1.5 


10 


20 


29(2) 


9(2) 


253 


1(2) 


4(2) 


254 


1{3) 




256 


16(2) 




262 


20(3) 


76(3) 


265 


1(2) 





2(3) 



12 27(2) 



61(2) 



When more than one sample was taken at any location, the geometric mean 
is reported and the number of samples is recorded in parenthesis. 
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TABLE 9 

Hamilton Harbour Survey Geometrical Means of Fecal Conforms/ 100 ml 



Station Year Year Year Year 
Number 1970 1971 1972 1974 

20 831(8) 246(5) 319(14) 260(1) 

256 1026(7) 359(5) 488(14) 500(1) 

262 355(8) 286(5) 494(15) 364(1) 

265 51(8) 12(5) 91(14) 18(2) 



In parenthesis are the number of samples used for calculation of the mean. 
1970 and 1971 data are from water pollution surveys of Hamilton Municipal 
Laboratories. 1972 and 1974 data are from Ministry of Environment surveys 
at a depth of 1.5 m. 
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TABLE 10 

Hamilton Harbour Survey t-Test of Fecal Coliforms/100 ml. for 1970-74 



Comparison Between Means for Years 

Station 1970-71 1971-72 1972-74 1970-74 

Number D.F. t D.F. t D.F. t D.F. t 

20 11 0.919 17 1.461 

256 10 0.738 17 1.408 

262 11 1.118 17 1.036 

265 11 0.948 13 1.200 14 0.651 13 0.6911 

(256+20) 38 0.656 24 0.5308 

(262+265) 30 1.007 27 0.5776 
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O 6RAQ SAMPLING LOCATIONS 

® LOCATION OF RECORDING CHEMISTRY MONITORS 



FIG. I 1974 SAMPLING LOCATIONS 

1 \ \ } ) \ } 



^ i 



FIGURE 2 



DISSOLVED OXYGEN RESULTS 
OFFSHORE LOCRTION HYPOLIMNION 
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FIG. 3 - TEMPERATURE PROFILE AT OFFSHORE 
MONITORING SITE, (M6 AND M7) 
JULY 16, 1974 AT 12 .16 HR. 



WATER QUALITY 
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HAMILTON 

HARBOUR 
STUDY 74 

section B 

Sediment statistical analysis 



HArilUnri HARBODP. STUDY 
SFCTinri R 
SEDF-iENT STATI5;TICAL ANALYSIS 



SUMMARY 



THE RELATIONSHIP BETWEEN SEDIMENT SAMPLES HAS BEEN 
STATISTICALLY ANALYZED BY COMPARING THE RATIOS OF 
CONCENTRATIONS OF TRACE POLLUTANTS. THE SENSITIVITY 
OF THE TEST WAS ESTABLISHED BY TESTING RANDOM NUMBERS 
FROM LOG-NORMAL AND UNIFORM DISTRIBUTIONS/ AND THE 
METHOD WAS APPLIED TO J/^7'^ HAMILTON HARBOUR SEDIMENT 
DATA. ALTHOUGH THE P7"> DATA WERE INADEQUATE TO 
PRESENT CLEARLY DEFINED RELATIONSHIPS, THEY SUGGEST 
SIMILAR ORIGINS FOR THE SEDIMENTS ALONG THE EASTERN 
AND NORTHEASTERN EDGES OF THE HARBOUR/ IN AGREEMENT 
WITH OSCILLATORY CURRENTS IN THIS AREA PREDICTED BY 
THE NUMERICAL MODEL. 



Correlation of Hamilton Harbour Sediment Analyses 

Tn July 1972, as part of the Hamilton Harbour Study, a series 
of bottom sediment samples was collected and analyzed for PCD, 
lead, chromium, cadmium, mercury and iron. In the initial 
discussion of these samples (Mor, 1974), it was stated that 
high concentrations of heavy metals occurred, with the worst 
locations near station 262 and the Hamilton l^TPCP outfall (Sta- 
tion 26 8) . It was further recommended that the sources of 
heavy metal discharges be identified and some action be taken 
to reduce these loadings. 

To this end, a statistical technique called "ratio matching" 
(Anders, unpublished) has been modified and adapted for estab- 
lishing potential generic relationships. This technique is 
based on tlie fact that sediment samples of common origin will 
tend to have similar ratios of concentrations of trace pollu- 
tants (heavy metals, PCB, etc.) whereas the absolute concentra- 
tions of these individual pollutants may vary considerably due 
to dilution with inert materials such as Si02 and CaC03. This 
technique could also be used for developing generic relationships 
between individuals such as biological communities based on 
species populations. 

For each sample, Anders divides the concentration of each chemi- 
cal parameter by that of each other parameter to obtain the 
"concentration ratio matrix", (X^ j } , which is a triangular mat- 
rix of size m X m, where m is the number of chemical parameters 
involved. Two samples A and B are compared by means of the ratio 
matrix {Y-- (An) } , which is generated by dividing each element of 
the {Xj^ j }a matrix by the corresponding element of the (Xj^ j }g 
matrix. The ratio matrix of A and D is thus: 



(Y^j (AB) } = 



Xii (A) 
Xii (B) 



X21 (A) 
X21 (B) 



X31 (A) 

X31 (B) 

I 

I 
Xn^i (A) 

^mi (B) 



X22 (A) 
X22 (B) 



3 2 



(A) 



X32 (B) 



Xm2(S) 



Xmm(^) 
^mm^^) 



n-1 



ir 



I 



Here, each X^. . represents an element of an (x^ j } matrix, and 
there are (9) Y matrices, where n is the number of samples i 



nvolved. 



If two sediment samples have similar origin, then the ratios of 
the pollutant concentrations of interest will have similar values, 
and thus the Y. . values will be close to 1. Consequently, Anders 
proposed that the similarity of each pair of samples be estimated 
by calculating the proportion of the matrix elements Yj^j (kl) 
which satisfy the conditions 



1 
M 



Y^^ (IcA) 



M 



where M is a variable parameter which he called the "matching 
criterion." Using this condition, each Y matrix produced a quantify 
called the similarity coefficient. 

These similarity coefficients Zj^j^ form a "correlation matrix" 
of size n x n, as follows: 



^^kO 



1 

Zab 1 



'An 



BC 



'Bn 



Anders suggested that the matching criterion M be chosen such 
that about 10% of the coefficients from random samples would 
exceed 0.33. Certainly the value of the method is limited by 
the need to choose a realistic value for M. If the difference 
of M from 1 was too large or too small, the statistical resolution 
of the procedure was seriously limited. 

Anders established the statistical significance of a given 
similarity coefficient by analyzing sets of randomly selected 
sediment samples from Michigan and California, using values of 
M = 1.5 and 1.3. He thus obtains a probability of obtaining a 
given similarity coefficient with random data as a function of 
the similarity coefficient and of M. However, the nature of the 
samples used to obtain such a probability is subject to questions- 
he observes higher probabilities of random coefficients when 
analyzing only Michigan samples as compared to that obtained with 
both Michigan and California samples. The difference is attributed, 
reasonably enough, to a geological similarity of the sample sites 
in Michigan. 
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Consequently, Anders' procedure bears two shortcomings: 
(a) the choice of M, and (b) the choice of "random" samples. 
As a result, one is really determining qualitatively which co- 
efficients, if any, are significant. 

In the past the main difficulty in the analysis of sediment data 
was the fact that the concentrations of the chemical parameters 
of interest could not be directly employed. Anders has overcome 
this problem by implementing the Y concentration ratio matrix 
for two samples. Lastly he has offered valuable insight towards 
a means of establishing a similarity coefficient from the Y matrix 

In an effort to overcome the difficulties listed above, an alter- 
nate method for analyzing the Y matrix has been devised, and is 
presented here. 

The similarity coefficient is redefined as an evaluation of the 
relative differences of all (Y^ j } (ki) matrix elements from 1. 
For each of the i^) off -diagonal elements in the (Yj^ j } (kl) matrix, 
the absolute value of the logarithm is used as a measure of the 
difference of the element from 1. As the similarity is inversely 
related to the difference of each element from 1, and a value of 
1.0 is desired for a perfect correlation (all (Y^-;} values = 1.0), 
1.0 is added to the logarithms, and the similarity coefficient 
takes the form 



'kl == 



m 


m 


Z 


T. 


i=i 


j-i+1 



£n {Yij} (kZ) 



+ 1 



(") 



In this equation, k and Z vary from 2 to n , where n is the 



number of samples, %j^}^ = 1 and Z, « = Z 



ilk 



Thus a triangular matrix 
of size n x n', similar in form to^Anders' Z matrix, is obtained. 

As previously discussed, the other problem in Anders' method is 
the selection of random samples in order to define a statistical 
significance for given similarity coefficients. We have used a 
random number generator to produce random data from which similarity 
coefficients have been calculated and statistically analyzed. The 
nature of the distribution from which random data was generated, 
the number of "stations" and the number of "chemical parameters" 
were all varied. Following this, sediment data provided by Anders 
was analyzed to provide a comparison with his method, and the 1972 
Hamilton Harbour results are also presented. 

RESULTS 

The results of several sets of random data analyses are presented 
in Table 1. In the first attempts to generate random similarity 
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coefficients, normal distributions with relatively small standard 
deviations (10-2 0% of mean) were used. However, the samples were 
so similar to each other that all similarity coefficients were 
significant; indeed, the means of the coefficients were greater 
than those observed with real data. A better distribution of 
similarity coefficients was obtained with two other distributions: 
a log-normal distribution with the natural log of the mean and 
standard deviations each equal to 1, and a uniform distribution 
with values between and 100 ppm for each parameter. with either 
of these distributions, the mean similarity coefficient varied 
between 0.4 5 and 0.54, independent of the number of stations or 
the number of parameters. However, the standard deviations in- 
creased from about 0.07 to about 0.12 as the number of chemical 
parameters was decreased from 8 to 4 , while remaining independent 
of the number of stations involved. ^s a measure of statistical 
significance, it was decided to use tlio narametric 7 test (Menden- 
hall and Scheaffcr, 1<^71) . For tootiniy the !iyTx>tliosi r. llmt r ^ f n 
the test takes the form 



7. = 0.5 



1-r ) I 1-po 



In this case we are testing whether any similarity coefficient 
is different from zero (po = 0) , v/here n is the number of chemicnl 
parameters involved. Alt hou gh tiie original formula of Mendenhall 
and Scheaffer contained /n*3 as a multiplier, v''n-2 was selected 
here as it produces a better fit v;ith random data ov^r a range of 
values for n. If Z>:^a/2, where a is the significance level selected, 
the similarity coefficient is indicative of a significant relation- 
ship between the two sampling stations involved. 

From Table 1, it is seen that the 1-iighest similarity coefficients 
observed in the uniform and log-normal distributions reoresent 
probabilities of 90-97% for difference from zero with 8 chemical 
parameters, and 76-96% for 4 parameters. Accordingly, the 95% 
probability level was selected as a limiting level for statistical 
significance; this level is exceeded by not more than 3% of the 
similarity coefficients in any random data matrix constructed from 
the uniform or log-normal distributions. 

Although maximum similarity coefficients tend to decrease with 
decreasing number of sampling stations, this is merely a function 
of the decreased similarity coefficient matrix size, because the 
use of a subset of a given set of samples will merely produce the 
corresponding subset of similarity coefficients from the original 
Z matrix; in addition, mean similarity coefficients in Table 1 
are generally independent of the number of samples. 

As the use of the parametric Z test assumes that the similarity 
coefficients are normally distributed, this hypothesis was tested 
using the Kolmogorov Smirnov test (IBM, 19 70; Poulton and Palmer, 
1973). Typical probabilities, 0.557 for several sets of random 
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samples, and 0.464 for coefficients derived from Anders' data and 
Hamilton Harbour data, strongly support the assumption of normality. 

To provide a comparison with Anders' method, his examples were 
re-analysed using the method presented here. Anders' sediment 
analysis data are presented in Table 2. The resulting similarity 
coefficient matrices are given in Tables 3, 5 and 7, and the 
probabilities derived from the Z test are given in Tables 4, 6 
and 8. For example 1, only one of the two relationships found by 
Anders was significant by this method. If the relationship between 
the stations (Figure 2) is accurate, this conclusion is reasonable, 
as sediment from the river (location 42 91) could be carried through 
the lake to the exit (4307). The two methods are in agreement for 
examples 2 and 3; in the latter, the most significant correlations 
are again between the deep water stations. The occurrence of 
several less significant correlations (in agreement with Anders) 
may be due to geological similarities as these samples were collected 
in an area remote from any heavy metal discharges. In this case 
Anders said that a somewhat higher significance level for randomness 
should be used; however, no statistical validation of this state- 
ment was offered. 

Results and Discussion of Hamilton Harbour Data 

Analysis of sediment samples collected in 1972 are presented in 
Table 9. Sampling locations are given in Figure 1. Unfortunately, 
not all samples were analyzed for cadmium and mercury; consequently, 
the data were analyzed twice, using 4 parameters for all stations, 
and using 6 parameters for stations 20, 256, 257, 262 and 269. The 
similarity coefficients and Z-tt^ p jjLabilities are given in Tables 
10-13. In the case of the comparison for 4 parameters, none of the 
similarity coefficients exceeded the 95% limiting level, which is 
0.882. The probability associated with stations 20 and 264 was just 
below 95% (Table 11) . As these stations are at the opposite ends of 
the harbour, it is difficult to justify this finding. As station 
264 could be associated with solids derived from the Dundas WPCP , 
and station 2 with the Hamilton WPCP, a relationship might be ex- 
pected, but solids from the Dundas WPCP should be mainly deposited 
in Cootes Paradise, (immediately west of Hamilton Harbour) and have 
little effect on the bay itself. In addition, solids from the 
Hamilton WPCP should affect station 268 more than station 20. The 
lack of a significant relationship between these stations suggests 
that solids from the Hamilton WPCP are largely deposited within the 
Windermere Basin. 

In the 6-parameter case (Tables 12-13) , the 95% limiting level for 
significance is 0.753. This case reveals the existence of signifi- 
cant relationships between stations 20 and 269, and between 257 and 
269. As the numerical model (MOE , 1974, section E) indicates major 
oscillatory current activity along the eastern and northeastern 
shore, it is reasonable to expect sediment movement within an arc 
extending between stations 20, 269 and 257. 

The use of a statistical criterion dependent upon the number of 
chemical parameters is thus supported by the Hamilton Harbour results, 
as well as bv the random number data presented earlier. The result 
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of the 4 parameter survey indicates the difficulty in analyzing 
inadeauate data sets. It is concluded that a minimum of 5 or 6 
chemical parameters should be analyzed in order to provide a 
reasonable basis for determining the relationships between various 
points in the harbour. Increasing the number of sampling stations 
does not affect the significance of a similarity coefficient com- 
paring any two stations; it merely increases the sensitivity of the 
analysis by providing more points for intercomparison . 

Furthermore, the chemical parameters chosen should all be those 
which are known to be discharged to the harbour as compared to 
background substances such as Si02/ Ca or Mg which are naturally 
present in large quantities. The use of the latter could be ex- 
pected to increase the value of all similarity coefficients (such 
as in the random data sets with small relative standard deviation) , 
effectively destroying the basis for establishing statistical sig- 
nificance. 

It is possible that the nature of the sediment sample would intro- 
duce a variable factor into this method. Significant variations 
in metal concentrations may occur for relatively small changes in 
sampling position or depth. The effect of this factor on the sen- 
sitivity of the technique should be investigated by taking several 
samples from points separated by a few meters, and by analyzing 
samples taken from different depths of a core. 

In 1975, it is recommended that sediment samples be taken from at 
least ten locations within the harbour, and possibly at as many as 
twenty locations (the 17 location grid used for chemical sampling 
in 1972, plus stations 264, 268 and 1030). Sampling should prefer- 
ably be done both before and after the proposed mixing program is 
undertaken. All samples should be analyzed for at least five of 
the six chemical parameters used in 1972, plus organic carbon, COD, 
BOD and Kjeldahl nitrogen. The latter four parameters would allow 
for the presence of organic discharges, as well as being of predic- 
tive value in the oxygen uptake studies currently in progress. (See 
Section E, Sediment Dissolved Oxygen Uptake) . 



n-6 



Acknowledgement 

Many thanks are due to the assistance of Mr. B. Noble, who 
performed most of the computer programming associated with this 
study. 



B-7 



R eferences 

Anders, O.U. , unpublished. Ratio matching: a statistical 
aid for discovering generic relationships among 
samples. Dow Chemical Co. , Midland, Mich. 

IBM, 1970. System 360 Scientific Subroutine Package, 
Version III, Programmer's Manual. 

Mendenhall, W. and Scheaffer, R.L. 1973. Mathematical statistics 
with applications, pp. 419-423. Duxbury Press, 
North Scituate, Mass. 

MOE, 1974. Hamilton Harbour Study. Water Quality Branch, 
Toronto, Ont. 

Poulton, D.J. and Palmer, M.D., 1973. Water chemistry data 
for the St. Clair River at Corunna, Ontario, as 
determined by continuous water quality monitor. 
Proc. 16th Conf. Great Lakes Res., p. 309-320. 



B-8 



TABLE 1 

SIMILARITY COEFFICIENTS CALCULATED FROM 
RANDOM DATA SETS 



Characteristic 


No. of 


No. of 


of Data Set 


Stn. 
13 


Parameters 


Normal, each par- 


8 


cuneter had separ- 


8 


8 


ate mean in range 






0-100 ppm. 






Normal, mean 


13 


8 


10 ppm, SD=2 ppm 






Uniform distribu- 


16 


8 


tion in range 


12 


8 


0-100 ppm 


10 


8 




8 


8 




14 


6 




12 


6 




10 


6 




S 


6 




14 


4 




12 


4 




10 


4 




8 


4 


Log-normal dis- 


16 


8 


tribution , In 


12 


8 


mean = 1, In Std. 


10 


8 


Dev. = 1 


8 


8 




14 


6 




12 


6 




10 


6 




8 


6 




14 


4 




12 


4 




10 


4 




8 


4 



Probability of 
Similarity Coefficients difference 
Std. from Zero for 

Mean Dev. Max. Maximum 



0.6606 0.0682 0.8253 
0.7065 0.0707 0.8277 



0.8709 0.0333 0.9381 



0.5262 
0.4526 
0.4650 
0.5063 
0.5288 
0.4834 
0.5294 
0.4880 
0.5391 
0.4981 
V.4 785 
0.4607 

0.4804 
0.4543 
0.4933 
0.4725 
0.4710 
0.4889 
0.4874 
0.5048 
0.4528 
0.4668 
0.5041 
0.4314 



0.0788 
0.0909 
0.0807 
0.0750 
0.0971 
0.0932 
0.0713 
0.0706 
0.1281 
0.1331 
0.1021 
0.1045 

0.0585 
0.0523 
0.0602 
0.0607 
0.0834 
0.0824 
0.0810 
0.0570 
0.1225 
0.1123 
0.1207 
0.0900 



0.7222 
0.6594 
0.6818 
0.6765 
0.7332 
0.7586 
0.6998 
0.6508 
0.9042 
0.8696 
0.7372 
0.7132 

0.6672 
0.5881 
0.6367 
0.5924 
0.8012 
0.7501 
0.6813 
0.6110 
0.8730 
0.7968 
0.8384 
0.6811 



0.991 
0.992 



0.998 



0.974 
0.948 
C.958 
0.956 
0.939 
0.953 
0.917 
0.879 
0.965 
0.940 
0.818 
0.812 

0.952 
0.902 
0.935 
0.905 
0.972 
0.948 
0.904 
0.845 
0.943 
0.877 
0.914 
0.760 
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TABLE 2 
DATA ANALYZED BY ANDERS 



Example Sample As Sb C u Zn Au Cd H^ C^ 



4293 


3.5 


0.2 


12 


67 


.013 


.88 


.72 


24 


4291 


8.3 


1.7 


298 


386 


.075 


7.8 


.86 


31 


4307 


12 


0.75 


124 


182 


.040 


4.0 


.41 


64 



II 



4298 


3.9 


.263 


34.3 


72.9 


.025 


1.26 


1 


.24 


15.1 


4305 


7.8 


1.4 


29 


211 


.022 


2.52 




.26 


83 


4303 


6.4 


.497 


83.3 


267 


.012 


7.14 




.40 


50.8 


4399 


5.92 


.296 


7.07 


39.8 


.0101 


.645 




.16 3 


17.8 



III 



21 


3.75 


.053 


1.77 


10 


.00075 


.070 


.014 


3.74 


22 


2.62 


.12 


7.40 


34.3 


.00122 


.27 


.009 


10.64 


23 


6.3 


.32 


35.0 


112.0 


.00361 


.873 


.148 


45.50 


24 


6.68 


.095 


1.83 


16.3 


.00054 


.030 


.021 


5.01 


25 


4.73 


.110 


13.1 


48,5 


.00047 


.412 


.076 


14.30 


26 


3.50 


.041 


0.65 


5.0 


.00005 


.030 


.008 


1.92 


27 


7.25 


.491 


25.9 


109.0 


.00374 


.727 


.122 


38.90 


28 


12.60 


.696 


52.0 


206. 


.00600 


1.660 


.234 


71.70 


29 


4.70 


.016 


1.90 


11.5 


.00600 


.100 


.009 


0.017 


30 


2.74 


.036 


1.83 


9.6 


.00016 


.064 


.008 


3.30 



NOTE: All concentrations are in ppm. 
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Table "^ 

MATRIX OF SIMILARITY COEFFICIENTS 
EACH COEFFICIENT IS BASED ON A SCALE FROM TO I 



1 


4293 


4291 1 


4307 


I 


42931 


I. 0000 


0.5141 1 


0.6225 


1 


42911 




1.0000 ) 


0.6971 


1 


43071 


1 


I 


1.0000 


1 
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Table '! 

MATRIX OF PROBABILITY STATISTICS FOR SIMILARITY COEFFICIFNTS 



I 4293 I 4291 I 4307 I 



42931 1.0000 1 0,8361 | 


0,9259 1 


42911 I 1,0000 1 


0.9652 1 


43071 1 1 


1.0000 1 
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Table 5 

MATRIX OF SIMILARITY COEFFICIENTS 
EACH COEFFICIENT IS BASED ON A SCALE FROM TO 1 



1 


4298 


1 4305 


4303 


4399 1 


4298 1 


I. 0000 


1 0.5109 


0.5292 


1 0.5459 1 


43051 




1.0000 1 


0.5645 


0.6647 1 


4303 1 




1 


1.0000 


0.5255 [ 


43991 1 1 


1.0000 1 
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Table ^^ 

MATRIX CF PROBABILITY STATISTICS FOR SIMILARITY COEFFICIENTS 



1 


4298 


1 4305 


4303 


4399 1 


42981 


I. 0000 


1 0.8328 


0.8509 


0.8665 1 


4305 1 




1 1.0000 


0.8827 


0.9503 1 


43031 




1 


I. 0000 


0.8473 1 


4399 1 1 




1.0000 I 
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If 



l[ 



u 



If 



If 



If 



Table 7 ^ 

matrix of similarity coefficients 
each coefficient is based on a scale from to 1 



I 





1 21 


I 


22 




23 


1 24 


1 25 


1 26 


1 27 




28 


I 29 


1 30 1 


21 


1 I. 0000 


1 


0.5962 




0.6070 


1 0.632 9 


1 0.5629 


1 0.6355 


1 0.6725 




0.6182 


1 0.4604 


1 0.7142 I 


22 




1 


I. 0000 




0.6924 


1 0.5167 


0.6064 


1 0.4749 


1 0.7411 




0.7449 


1 0.3832 


1 0.6296 1 


23 




1 






1.0000 


1 0.5317 


1 0.6874 


1 0.5115 


1 0.8102 




0.8742 


1 0.4098 


f 0.6463 I 


24 




1 








1 1.0000 


0.5059 


1 0.6553 


1 0.5840 




0.5487 


0.3789 


1 0.6293 ( 


25 


till 


1.0000 


1 0=5779 


0.6594 




0.6889 


1 0.4140 


0.6338 1 


26 1 


1 1 1 




1.0000 


1 0.5419 




0.5199 


0.3838 


0.6453 1 


27 1 




1 




1 


1 


1 


1 




1.0000 




0.8616 


0.3895 


0.6974 1 


28 i 


1 1 1 j 










1.0000 


0.4059 


0.6914 1 


29 1 




I 




1 


1 


1 


1 




1 1 


1.0000 1 


0.4581 1 


30 1 




1 




1 


1 


1 


1 




1 1 


1 


1.0000 1 



Table ^. ^ 

MATRIX OF PRC3ABUITY STATISTICS FOR SIMILARITY COEFFICIENTS 



to 
1 





21 


22 


1 23 1 


24 




25 


1 26 1 


27 1 


28 1 


29- 1 


30 1 


21 1 


1.0000 


\ 0.9077 


I 0.9154 1 


0.9590 




0.8814 


i 0.9340 1 


0.9542 1 


0.9231 1 


0.7773 1 


C.9718 1 


2Z 1 




1 i.aooo 


1 C.9632 1 


0.8387 




0.9150 


1 0.7941 1 


0.9804 1 


0.9815 1 


0.6774 I 


0.9304 1 


23 1 


1 1 1.0000 


0.8533 




0.9611 


1 0.8334 


0.9943 1 


0.9991 1 


0.7138 


0.9404 1 


24 


1 i 


1.0000 




0.8277 


1 0.9454 


0.8985 


0.8690 1 


0.6713 


0.9302 I 


25 








1.0000 


1 0.8937 


0.9476 


1 0.9617 1 


0.7194 


0.9330 1 


26 




1 1 1 1.0000 


0.8628 


1 0.8419 


0.6783 


0.9398 1 


27 




1 t 1 


1 1.0000 


1 0.9985 


0.6863 


1 0.9653 1 


28 






1 




•1 




1 I. 0000 


1 0.7085 


1 0.9628 1 


29 




i t 1 






1 I. 0000 


1 0.7746 1 


30 




1 I 1 








1 1.0000 1 



I 



TABLE 9 

ANALYSES OF SEDIMENT, HAMILTON HARBOUR 
July 19, 1972 



Stat;ion 



264 15 



262 



257 255 256 269 1030 267 20 



268 



PCB (ppm) 
Lead (ppm) 
Chromium 

(ppm) 
Cadmium 

(ppm) 
Mercury 

(ppm) 
Iron % 



.80 2.0 
100 420 
40 210 



3.0 
930 
150 



1.7 



.58 
170 
120 



18.5 8 



.81 



.05 

25 

27 



.17 
55 
4 7 

5.0 



1.3 
140 
97 



21 .7 



.22 

60 

13 



.26 

60 

41 



.66 10.1 
80 340 
43 295 

5.0 

.4 



2.86.1 9.4 4.4 1.6 3.3 3.8 1.4 2.3 2.4 3.8 



I 
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Table T * 

MATRIX OF SIMILARITY COEFFICIENTS 
EACH COEFFICIENT IS BASED ON A SCALE FROM TO I 





1 


2 64 


1 


15 1 


262 


1 257 


! 255 1 


236 I 


269 1 


1030 1 


267 1 


20 j 


268 1 




264 


1.0000 


1 


0.6775 1 


0.7240 


I 0.6143 


1 0.4888 1 


0.5638 1 


0.7800 


0.7015 1 


0.6392 


0.R745 1 


0.4748 I 




15 


1 1.0000 1 


0.6727 


1 0.6763 


1 0.4607 


0.5255 


0.7405 


0.6100 


0.6636 


0.7150 1 


0.5173 1 




262 1 


1 1 


1.0000 


1 0.5769 


1 0.4861 


0.5524 1 


0.6209 


0.7157 1 


0.5777 


0.6676 


0.4917 1 


t3 


257 1 


1 1 




1 1.0000 


1 0.5790 


0.6983 1 


0.7301 


0.6788 1 


0.8081 


0.6554 1 


0.4609 1 


00 


255 


1 


I ! 1.0000 


0.7612 


0.4938 


0.4706 


0.6241 ! 


0.5252 


0.3817 I 




256 


1 


1 i 


1.0000 


0.5451 


1 0.5594 


0.7121 


0.5624 


0.4243 I 




269 


1 1 


[ ) 1 


1 


1.0000 


1 0.6200 


0.6799 


0.8401 


0.5210 1 




1030 


1 1 


1 1 


1 




1 1.0000 


1 0.6421 


1 0.6718 


1 0.4144 1 




267 


1 1 


1 i 1 


1 






1 I. 0000 


1 0.7122 


1 0.4357 1 




20 


1 1 




J 


t 


1 




1 


1 I. 0000 


1 0.4807 1 




268 


1 ( 


1 1 1 1 


1 1 


1 




1 1.0000 1 



Table 11 

MATRIX OF PRCBABILITY STATISTICS FOR SIMILARITY COEFFICIFNTS 



1 264 I 15 1 262 1 257 I 255 | 256 1 269 I 1030 I 267 | 20 I 268 I 


264 1 1,0000 1 0.7564 | 0.8048 I 0.6892 I 0.5502 I 0.6334 | 0.8607 | 0.7816 I 0.7155 1 0.9441 | 0.5346 | 


15 1 1 1.0000 1 0.7513 1 0.7551 | 0.5190 | 0.5910 | 0.8217 | 0.6839 I 0.7416 | 0.7955 | 0.5820 1 


262 1 1 1 1.0000 1 0.6477 I 0.5472 1 0.6208 | 0.6958 I 0.7962 I 0.6487 | 0.7459 | 0.5535 1 


257 ill! 1.0000 I 0.6501 | 0.7783 | 0.8111 I 0.7578 1 0.8873 | 0.7329 | 0.5191 I 


o255 1 1 I 1 1 1. 0000 1 0.8 ^23 | 0.5559 | 0.5300 I 0.6993 I 0.5907 | 0.4303 | 


S256 1 1 1 1 1 1 1.0000 1 0.6127 | 0.6285 1 0.7925 I 0.6319 \ 0.4782 I 


269 1 1 1 i 1 i 1 I. 0000 1 0.6948 j 0.7589 | 0.9159 | 0.5860 | 


10301 1 1 1 i 1 1 1 1.0000 1 0.7186 | 0.7504 | 0.4671 | 
267 1 t 1 1 1 1 1 1 i 1-0000 1 0.7926 1 0.4910 | 
20 1 1 1 1 I 1 1 1 1 1 1.0000 1 0.5412 | 


268 i 1 1 1 1 I 1 1 1 1 1 1.0000 1 



THE PL/1 STEP OF M0ES6D HAS EXECUTED COMPLETELY 



Table P 

tlAI3iLIflN-tlABfiOilE-i£DIfi!£blI-DAIA-z=-6-£AEAtlEI£Ki 
MATRIX OF SIMILARITY COEFFICIENTS 
EACH COEFFICIENT IS BASED ON A SCALE FROM TO I 





262 


1 


257 


256 


269 


10 1 


262 


1.0000 


1 


0,6458 


0.5648 


0.6937 


0.7510 1 


257 


1 I. 0000 


0.6882 


0.7905 


0.7223 1 


2 56 


1 


I. 0000 


0.5860 


0.6180 1 


2 69 


1 




1.0000 


0.8106 1 


20 




1 


1 






1.0000 1 
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Table 1^ 

MATRIX OF PROBABILITY STATISTICS FOR SIMILARITY COEFFICIENTS 





262 


257 


256 


269 




20 1 


262 


I. 0000 


0.8755 


0.7993 


0.9128 




0.9489 1 


2 57 




I. 0000 


0,9088 


0.9681 




0.9320 1 


2 56 






UOOOO 


0.8208 




0.8511 1 


269 








1.0000 




0.9760 1 


20 








1 1.0000 1 



B-21 



79° 50' 



03 
I 



43° 16' 





i 


loop 500 Q 




1000 


2000 METRES 










\ 






^ 


.»*".••"■■•• ... 1 


, 














■'■'-"" v^'""^ ••'-•■"■■■ 




V 


Oy^O .050 








/^ ©»' 






Vv 


^r; 


^ 






1 


X* '"'"^■■■■/" 






■^ ^ 


A 


\ 






y 


/^.^s- 




1 

- 


.. o** 




^ \ 






• /^'• 


:--^ V ■■■ 










\ \ 




1 X • '*"•♦ ' ■ 


;^^" 








0"' 


1- 


A\ 


y 


I 264 // 


t 

; 6 


f 


.*' 


i / 


.yv^^^'^®" 


i 


Y A 


i 




: 




C/ 


A 




r\ \ 










I 


y 


r Ij 


/ 


/j)/K ^\ 


^ ' ^\^ 



FIGURED HAMILTON HARBOUR SEDIMENT SAMPLING LOCATIONS 

JULY 19, 1972 



! 



\ 



9W«. »"^ 54*» 




FIG. 2 SITES OF SAMPLES COMPARED IN ANDERS' EXAMPLE I 



APPENDIX I - COMPUTER PROGRAM 

A program for the evaluation of similarity coefficients has been 
developed and is named MOESED. The job consists of two steps, 
the first written in PL/1 computes the coefficients, and the 
latter is a FORTRAN version of the Kolmogorov-Smirnov Test from 
the IBM System/360 Scientific Subroutine Package. 

The coefficients generated by the PL/1 step are both printed 
out and temporarily stored on disk to be used as input for the 
FORTRAN step where they are tested for normality. Matrices of 
Z statistics and probability statistics are also printed. 

The PL/1 step can accommodate an infinite number of independent 
data sets. It is intended that as more and more sediment data 
becomes available, their data decks will be added to the. prograiri 
so that the maximum number of similarity coefficients can be 
used in the Kolmogorov-Smirnov Test. 

The following pages include a discussion of the program flow and 
restrictions, followed by an explanation of the required input 
(control cards and data) . Finally, the complete source listing 
of the PL/1 program is presented. 
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PROGRAM FLOW 

1. Read in title of run 

2. Read in: a) Indicators as to whether or not to 

print the computed X and Y matrices. 

b) critical probability value (use a nega- 
tive number if no particular value is of 
interest ) 

. c) number of measuring stations and chemical 
parameters . 

d) names of chemical parameters and stations 

3. Read In data for one station. 

4. Print input data. 

5. Compute X matrix. 

6. Output X matrix (if indicated). 

7- Go back to step 3 until data from all stations are 
read in. 

8. Compute Y matrix from X. and X. =, 

9. Output Y matrix (if indicated). 



10. Go back to step 6 until the Y matrix is formed 

between the Xj^ and X^ matrices (n Is the subscript 
of the last X matrix ) . 



11. Compute a matrix of similarity coefficients, store 
and print . 

12. Compute Z test statistic for similarity coefficients 
derived from Y matrices and print. 

13. Compute probability statistic and print. Check for 
significant probabilities as determined by crit-Z 
and underscore. 

1^. Return to step 1 if end of file condition does not 
exist . 

15. Execute step B - Kologmorov-Smirnov Test using the 
similarity coefficients previously stored on disk. 
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RESTRICTIONS : 

1. The title for each independent data set run is restricted 
to a character string of maximum length 80. 

2. The names of the measuring stations are restricted to 
character strings of maximum length 4. Any character on 
the keypunch can be used. 

3. The names for the parameters measured are restricted 
to character strings of maximum length 4. 

4. The number of stations and parameters measured are limited 
by the core available for execution and secondly the line 
width of the printer. The results are printed out in the 
form of matrices; hence if the line width (132 spaces) is 
exceeded the output will be incomprehensible to one 
unfamiliar with the program. To avoid this the number of 
stations should not exceed 14 and the number of parameters 
should not exceed 12 without altering the output formats. 
Note the size of the arrays are dynamically determined by 
program inputs. 

5. There must be an equal number of data items available 
for each station. For example, consider the following 
hypothetical data:- 

Station A Station B 

Zn 5 ppm Zn 3.2 ppm 

Cu 3 ppm Cu 4 . 8 ppm 

Cd 5.2 ppm 

The above data is insufficient and the program will fail 
unless a measurement for Cd is taken at Station A. 

6. The FORTRAN Step is limited to 700 similarity coefficients. 






PROGRAM DOCUMENTATION 

Documentation occurs within the program in the form 
/* comment */. Extra documentation has been included 
for parts which differ from FORTRAN. 
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INPUT/OUTPUT 



STN-NO 
PARAM 

STNS (1: STN-NO) 

PARAMS (1: PARAM) 

X (1: STN-NO, 1: PARAM, 
1: PARAM) 



Y (1:X, 1: PARAM, 
1: PARAM) 



INSTRUCT 1, INSTRUCT 2 

MASTER 
NAME 

COEFF (1:X) 

CRIT-Z 



number of data measuring stations 

number of parameters measured at 
each station 

array of names of stations 

array of names of attributes 

a 3-dlmenslonal array storing the 
computed concentration ratio 
matrices for all stations 

a 3-dlmenslonal array storing the 
computed Y ratio matrices of two 
samples N.B. X = (num*<*2-num)/2+num 

character strings of a maximum 
length of 3 for Instruction of 
whether or not to output X and Y 
matrices respectively. 

disk file for temporary storage of 
all similarity coefficients computed 

character string of a maximum 
length of 80 to label output for 
each Independent run of data sets 

a 1-dlmenslonal array used for 
similarity coefficients, Z-values 
and probability statistics 

critical probability statistic to 
be used to underscore significant 
values (N.B. a negative dummy value 
may be substituted) 
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CONTROL CARDS 

Either the PL/1 F or optimizing compilers can be used. The 
optimizing version is preferable because it is more advanced 
than its predecessor F. The program is set up to use the 
PL/1 F compiler. 

The following is the arrangement of the job and the job control 
language: 

// MOESED JOB (P, 0123456, 789) , 'JIR, NAME* ,CLASS«A 

// A EXEC PL1LFCLG,PARM.PLIL='SM= (2,80,1) ,ATR' 

// PLIL.SYSIN DD * 

PL/1 PROGRAM 
// GO. MASTER DD UNIT=DISK, DISP= ( ,PASS) , SPACE= (TRK, (10 , 11) , RLSE) 
// GO. SYS IN DD * 

PL/1 DATA DECK 
// B EXEC FORTGCG 
// FORT. SYS IN DD * 

FORTRAN PROGRAM 
// GO.FTOIFOOI DD DSN=* .A. GO .MASTER, DISP* (OLD , DELETE ) ,UNIT=DTSK 
// GO.SYSIN DD * 
// 
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DATA_OUTPUT 

The output is readily comprehensibln providinq the rostrictions 
previously notrd are observed. Tn tlie event the restrictions 
are deliberately ignored a single row of a given matrix will 
appear over two or more lines of printing. (Also refer to data 
input for output options.) 

DATA IiJPUT 

Unlike FORTRA:j, data input does not require any special format 
in PL/1. For certain inputs, formats may be advantageous, 
however, the program SED does not use any. In light of this, 
data can be entered in any position on the card so long as it 
is separated from a successive item by either a comma or a blank 
or both. 
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DATA CARDS 

DESCRIPTION 

1. On the first data card, starting In column 1, enter the 
title of the data set (o.lphanumeri c of maximum lenr^th 
80 characters ) . 

2. On the second data card stating in column 1, punch 
— YES — to print the X matrix (else punch — NO — ). 
Leave one blank space and repeat with either --YES — 
or --NO — for control of printing of the Y matrix. 
Leave one blank space, and then punch a critical r 
value. Note that the critical F value results in heavy 
underscoring of significant probability statistics. A 
dummy negative value may be used. 

3. On the third data card enter the number of surveyinr 
stations, followed by a comma, and then followed by the 
number of parameters examined. 

k. On the fourth data card, starting in column 1, enter 
the surveying station names, each name followed by a 
blank space. Note each name must be a string of 4 
characters, thus if the name is less than h characters 
add blanks to the end. 

5» On the fifth data card, starting in column 1, enter the 
names of the measured parameters, each name followed by 
a blank space. Note each name must be a string of 4 
characters. 

6. On the sixth (and additional data cards if required) 
punch in the data, station by station, each number 
followed by a comma. 

7. Repeat steps 1 through 6 for another data set. 

NOTE: No data input is required for the FORTRAN step. 
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S£D:»R0CEOURF OPT lONSl MAINJ ; 



000000 10 



PAGP 



<;ED:PR0CEDU»E aPTIOSS(MAIN) : 00000010 

/****^»»»»«;^««**«*#ft»«*****3r*** *****************♦*********#»********** /00000020 






/* 
/* 
/* 

/* 
/* 

/* 
/* 
/* 
/* 
/* 
/* 
/* 
/* 
/* 

/* 
/* 

/♦ 
/* 
/« 
/* 
/* 
/# 
/* 
/* 

/» 

/* 

./■* 

/* 
/* 
/• 
/• 
/« 
/• 
/# 
/» 



^CESEO IS A TWO STEP JOB- THE FI«ST STEP IN PL/1 CALCULATES THE 
X AND Y RATIO MATPICES, THE SIMILARITY COEFFICIENT, THE Z TEST 
VALUE. AND LASTLV the PR09A?ILITY 3^ 5 I GN I ^^ I C ANC E . THE SIMILAR- 
ITY COEFFICIENTS tSE ''LACED CN DISK, TO ^E USED AS INPUT pO' THE 
FOR THE SECOND STEd- thE IBM SCIENTIFIC SUBROUTINE PACKAGE 
VERSION OF TH'P KOL MQGOROiV- SM I RNOV Tt-^T I IN FQiRTPAK ), 



THE PROGRAM IS 
TO IMPROVE THE 
IS ADVISABLE NOT 
TO ADD ON TO IT. 
^ND THE FORMAT |S 



DESIGNED TO HANDLE AS M.ANY DATA SETS AS AVAILABLE 
RESULTS OF THE iKOLMOGOROV-S^IRNQV TEST, HSNCE IT 
TO REMOVE THE PREVIOUS DATA DECK, BUT INSTEAD 



THUS, ONLY 
AS FQLLOtaS 



THE PL/1 STEP REQUIRE": DATA IN^UT 



ON THE FIRST DATA CARD, 
THE TITLE OF ThE DATA Si 
80 CHARACTERS) . 



STAOTING IN COL. 1, ^NTER THE 

T (ALPHANUMERIC CF mAXIMUN LENGTH 



i cni, 1, PUNCH 

E PUNCH — NO — ) . 
ilTH EITHER — YES — 



*/00000030 
*/G0000040 
*/000'^0050 
*/OOC00060 
*/00000070 
*/000'^0080 
*/ 00000090 
»/00000100 
* 700 0001 10 
*/ 00000 120 
*/ 00100 130 
*/ 00 000 140 
*/ 00000 150 
*/00000160 
*/0000ni70 

*/oooo'^igo 

*/00000190 
*/ 00000 ?00 
*/00000?10 
*/0000C220 
*/0O0O0230 



ON THE SECOND DATA CARD STATING 
• — YES — TO POINT THE X MATRIX (E 

LEAVE ONE BLANK SPACE AND REPEAT ^ITH EITHER — YES— OR 
-NO — FOR CONTROL OF PRINTING :f the V MfiTRIX. LEAVE ONE*/00000240 
BLANK SPACE, AND THEN PUNCH A C^^ITICAL PRCBABILITV VALUE. */O0OCG25O 
NOTE THAT The CRITICAL 2 VALUE -=SJLTS IN H^AVY UNDER- 
SCOPING OF SIGNIFICANT P'OBABILITV STATISTICS. A DUMMY 
NEGATIVE VALUE MAY 3E USED. 



4, 



6. 



ON The Third data card enter the number of surveying 

STATIONS, FOLLOWED 3Y A COMMA, AND Thfn FCLLOWEO BY THE 
.NUMiBER OF PARAMETERS EXAMINED , 

OiN The fourth data capo, starting in col. i, enter the 

SURVEYING STATION NAMES, EACH NAM:E FOLLOWED BT A BLANK 
space. note =ACH NAME MUST BE A STRING OF 4 CHARACTERS, 
THUS IF THE NAME IS LESS THAN 4 CHARACTERS ADD BLANKS TO 
THE END. 

ON THE FIFTH OATA CAOD, STARTING IN COL, 1» ENTER THE 
NAMES OF THE MEASURED PAOAMFTER.S, =ACH NAME FOLLOWED BY 
A BLANK SPACE. NOTF EACH NAME MUST BE A STRING OF 4 
CHARACTERS. 

ON THE SIXTH (AND ADDITIONAL DATA CAPOS IF REOUIPED) 
PUNCH IN THE DATA, STATION BY STATION, EACH NUMBER FOL- 
LOWED BY A COMMA. 



*/00000260 
*/00000?70 
*/0O00O28O 
*/000002Q0 
*/ 00000 3 00 

*/oonon3io 

*/00OC03 20 
*/00O00 330 
*/ 00000 3 40 
*/00OO0 350 
*/00000?60 
*/00000?T0 
*/00000380 
*/000003O0 
*/00000400 
*/00n00410 
*/00000420 
*/0000n430 
*/00000440 
*/0OO00450 
*/00000460 
♦/00000470 
*/00000480 



/********«**«**♦***#*****«**********!>**«****♦*$*♦*** ♦*********<E*******/Q()0004 90 

00000 500 
00000510 



^crj.oBTQPQysc :pt lONSt MiiN) : 



00000^10 



p tr, c 



■ I 



ft 
m 
m 
m 

11 

31 

ST 

m 



4# 
4i 






4f 



1^ 
tl 

ft 
61 



hZ 



ii 



PUT c ILEi: SYSPRINT) SKIP; 
00 1=1 TO PAR AMI BY I ; 

SFf FILE (SV SIM) LIST IXOE FF •( I 1) : . 

PUT FIiLE( SYSPRINTi SKIP EDIT ICOEFFIII," PPm , PAftaHS ( [ 
f Xt 57) ,F«9,3) ,&,xn l,AJ : 
EMO: 

DC 1=1 Tn PAkAM BY I ; 
DO Ksl TC ! BY 1 ; 

X{ J,I .O^COEFFd )/roPFF(K) ; 

cNO; 

IF IMSTeuCTl=' vES' TriEfl 

DO; 

• PUT F IL^l, SYSP='F^lT^ SKIPI2) cCH (' • CC-^put F «, 

• CONCENTS ATI ONi RailO Ni&TSIX FQ» ST JlT | Qtl' , STN S '! J 
I Xll I ,A,,A,Xi( I) ,A»: 
PUT F ILE t SYSPRINTI SKIP: 
DIJT c ULc ( SYSPP INT] SKIP EDIT ( ' | ■ , t PAS AMS tl ) , • 1 * 

DC 1=1 TT PARA'^) MX (51 , A, (^A^^^v) ix J 3) ,A,X(?) , A) 
I ; 
00 1 = 1 ^0 PARAVi BY I : 
H = '>ARAM-I ; 
PUT FILE (SYSOP INT ) 5KIP EOI"^ !(•-' 00 G= I TO F)) 

(XII), IF) A) : 
PUT FILE(5YSPP[NTj SKI^ =OIT ( o as A-»^S ( I) , « I ' . ( i< U 
,K),'|' no K = l TT J.CI' CC G=l TO H))i|xm 
, a ,, A , U ) ( r ( q , 4 1 , -^ ■ , 1 H ) ( X ( q I , A ) ) : 
PND : 

OUT F1L= ISYSP'o I'NiTI ^v.\o C3IT (('•-• DC G= I TO' F)) 
t K (1 ) , ( F J A ) ; 
END; 
eNO; 

/* CCM^uTE THE CCNCENTR AT tON t^ATfj «ATRIX FCO thE CQM- 

/* D&p.ISG^J Qc "'WO STATIONS AND 1= INST^^UCTEC tq, PRINT IT 

LI»«IT=0; 

DO 1=1 TO STN_NO BY I; 

DO J=l TO STN_NiO BY 1; 
LI»1IT = lIHIT+1 I 
00 K=l TO PARAM BY 1 ; 
00 L=l TO K BY 1 : 

vj L I M I T , K , L) = K ( I ,, K , L J / X ( J , K , L I ; 
END ; 
END; 

[F 1NSTRUCT2=' YFS' THEN 
DO; 

PUT FILEfSYSPPINT) SKIP(5) EOIT ('COMPARISON ', 
•RATIO MATRIX FOR STATIONS ',STNS(n,» VS. ' 
STNS(J) J(X (l) , A, A, At A, A); 
PUT FILEISYSPRINT) SKIP EDIT { • 1 ' . ( PAR AMS ( 6) , ' I ' 
DO G=l TO PAftAM) ) (X (5} , A, (PARAM)(X(3J, A,X( ?) 
A) J ; 
DO K=l Tc PAP AM BY 1; 
iH = PARAM-K; 
PUT FILEISYSPRINT) SKIP EOIT (("-• 00 &«i TO 

MX tn , (f}a) r 

PUT FILEISYSPRINT) SKIP EDIT IPABAMS (K ) t ' P » < 



00001090 

00001 100 
O'OOOIUO 
) lOOOOl 120 
OOOf^l 1 ?0 
OOOO1K.0 
00001150 
Ononii60 
00iOll"'O 

ooooiieo 

OOOOll'^O 
00001200 
00001210 
00001220 
)»000C12 30 
00001240 

00 001250 
0O0fil26P 
000ni?70 
00001280 
0000129" 
00001300 
0000131'" 
0OD01320 

, 100001330 
00001340 
000013 50 
00001 ?6n 
000013TO 
000013 90 
00001390 
00001400 
00001410 
*/COOO1420 
*/000"1430 
01001440 
00001450 
00001460 
03001470 
00001480 

00001 49'0 
00001500 
00001510 
000015 20 
00001530 
00001540 
00001550 

, 00001560 
00001570 
00001580 

, 00001590 
00001600 
00001610 
00001620 

F J 000016 30 
00001640 
00001650 



SEO:Ps;CCEiDURF O'PT I QNS f 'NA I N) : 



OiOOO'OOlO 
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6 7 
63 

6 9 

71 



7 7 
11 
74 
75 
75 

77 

7Q 

79 
an 

:P1 
P? 

S3 

3 4 

^^ 
5 7 

6 a' 
e,9 

11 



V {L IMI T.KtL) t • I ' OG L»l TO K 1 , ( " ! » 'DO' L = l 
TC H ) M X ( 1 I , A , a , ( K I I F ( 9 , 4 K A ) ♦ « H H X ( 9 ) , A J 

) ; 

F'MD; 

PUT PILEISYSO^i [NT| SKIP EDIT !!(■-• OH G=l TC 'P ) > 

f X I 1 ) , t F 1 A 1 ; 



E m ; 
EN'D; 



/* CALCULATE SlWlLASITV CQEpFICIENT AMD po I M IT */ 
CALL SIM(¥,CCEPP} ; 

CODE P= 11: 

CALL TITLE IfaA^tCtCODEP) I 
C0.DEP = 4; 

CALL P'^INT (CO£FP,STi^S,CPIT_P,CQiO'PB'l : 
/* CALCULATE Z_TEST VALUE mD PPIM IT */ 
CALL Z.TESTjCOEFF.PflPAW) ; 
CODER =2; 

CALL TITLE { N£'*1p ,C0DE3 ) ; 
C0iD'PP = 4; 

CALL P'BIMT(CDECc,STN,S,fB'IT_3,C0DE3|;. 

/* CALCULATE P9CBABILITY STATISTIC A^jD p^JMT IT */ 
CALL P_?TATlCnEFF| ; 
COiOEO = 3; 

CALL TrTLE<NA«£,CODEP) 1 

CALL P^INiT tCCEFP,SThjS,Ca!T_p,,,C0:D'EO) I 
EHD D STORE 5 

CND- ' ' 

FINALtOOt 

PASSER. = '?«;9999.0; 
WeiTE FILC |P*'4STEF| FRCP' (PASSFe); 
F NO ? 



00 001 660' 
OOO'^niSTO 

ooO'Oiseo 

0'OOi-9L6'90 

0'Q^0'01720 
0'aO'01T30 
00001740 
00001750 

00001760 
00001 770 
00001780 
OOC01790 

ooO'Oiaoo 
00001 g 10 

000018 20 
00001830 

0009 IP 40 
OOOniFSO 
0'0001860 
00001870 

oooniaeo 

0001° 00 

0O00191C 
00001920 

000019 30 
0000 ic 40 

000010 50 
00 001960 
000019 70 

O0'00i9eo 

00001°<?n 
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OO'OOOilOi 



PA.SF "7 



loa 



109 



in 


Hi 


1 1 3 


114 


115 


116 


117 


1,19 
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120 


121 
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129 


130 
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/* 
/* 

/* 



^ ^^L-^ ****;** 



*3^** ************ **:**♦ *:*t*^^* ******************** ******»****/ 



THE IMT[-:5NAL SUERCUTINE o^INT lABEl? 4\'i OUTPUTS fi TC|iMGuLflP 
MATHIX CCNT^aFO IN A 1 i2 [^^ftPA^T£R P'lN'^ F I r L 3 . '^H'- Pa^^i^FrEa 
MATRIX I '^. ft DNF DIMf'NSIONAL APRAY. IF A MATRI/ 0"= ^'POBAeiLITr 
STAT!<5TICb IS i^eiNG PRINTfO JINnifAT^n gv Si^TTinC ^hp paPAMf^T^=t 
C0DER1 = 11 Thf-N SIGNIFICANT PESIJLT^^ HILL ^^E UNirR^COrtfO IF A 
LlGNiriCAMT f5lT_P VALUE IS MASSED. 

-t*a*±*fta*#s:t***¥^**jea*J?******4±*w ****************************** <•"***" 



V-! :Di<i>c:'"j!?t ( MAT'S IX, LOCAL E, f' = !"^_'-' ,crOt-^l ) 



OeCLARt MATPIXI*J FLQiAT, 
LOCALE 4*1 CrtftiB:ii:4Ji 
COOEPl FIXED, 
CB'IT_P FLCiiAT, 
(SP.JM,.) FIXED, 



VAO'YrNIG, /• STATn-hi iKAMFS */ 

/• TITL* OiF IKATiBJX */ 
/* C-RITICAL P VAIU? */ 
/♦ CEMTII^RING VARIABLES */' 



IF,G,Hi„l, 12, J, J2,KfiR,,SI FIXED 13, 01; /• FOR liM^DEXING */ 

/* CO'NOUTC NCCESSART SPACI'MG "tQ CENTER '»<'ATq i; i . hCTE, IF THE S I 7 : -*=/ 
/* OF 54ATRIX EXCEiOS 12 X 12 THE OUTPU^ QF EACH RCH .^ILL APOEA=' -*/ 

/* ON ONS OP nO^i PRINTING LINES. */ 

Ji N* 113 2- %*m jU^O ( L OC 4i E f 1 } -% I /2 EO' J 

ic JN<Q THEIN JN=0'; 

■OUT FILE{Sy:50'.Q'Ihj|T» SKlPiSl S'DIT »« 1 • » < LCCAiL^ I I I , • | ' '''O 1=1 ^0 

HBOyND 'I LOCALE, 1)1} (XI JW+4 1 ,A, I HflOUMO I L3C AL E , 11 J I X I 2) , A, X ( 2 1 . A ) ) : 
iK = HB0UNOi('LOCALE»ll ; 
I>1| 
I2-n 

DO' J=l TC IHBQUMD! LOCALE, 11 nv I i 
R*K-I+1; 
H = HiB0UNO I iOC AL E , 1 1 ~» I 

OUT F ILE«SYSP«tlNT} SXIP 
P U'T F I L f t: S t SP« I ftiT Ji SK I P 



EO I T 

EDIT 



itMITP'III SI ,' :| ' DO 
IF(7,4I ,Xll},A,» || I 
IF C^IT_'p>ni ThE'Ni 

IF COOES 1 = 3 THEN 

SP-JN+4*iHI*9j *'■ 

00 Ji2 = l:2 TQi -K BY 11 

IF li^:.ITRIXI J21 >=CRIT 

90.; 



I Fl A) ; 

TO HJ , 

ssi TO KJ Hx i[:j!>jir,a,A, tH) ma JtAJ, iR » 



I LOCALE UK 'f 



Tn F} 1 I WUN) , 

, ( • I ' DO L = I 



¥ i, »iiATeiXl J2,Kl,'10'00l 



EOIT I"" 





PUT FILEtSfSPRIHT* S« HP 11^1 




SPli.iJi 




PUT iFILfiSYSPRIMTli SKlPfQl 




SPJ.AIM 




SP=SP*9| 




END; 


im ; 




6»| 1 





EDIT t 



THEN 



I=K*1| 



Q0'i0?390 

noaopAOQ 

00G0?410 

(noo?4?o 

"0002430 
00002440 
nOQr^2450 

n 002 4 60 
0ino:'470 
0002480 

j:o:i''^o 

OOO'^^'JOO 
00002510 
000T^5?0 
00002530 
00002540 
O'**002 550 
00002560 
00002570 
0^002530 
T ? 5 9 
00002'jOO 
OOOr^2&lO 

00032^40 
"'0255'^ 
000^2660 
0')002fc7C 
00002tP0 
3^769'^ 
■^0002700 
0000?71'i 
000-^2720 
00002730 
0^01^2740 
0T00.Z r5C 
00002760 
on 002 7 70 
00002780 
00002790 

oooo2eor 

O0J0281O 
00002820 
OTO02S3O 

000"2P40 
00 )0''S50 

oooo?efcO 

O0QO?e70 
0T002880 
O000?£*90 
00002900 
0000 2^10 
000n?O20 
00002930 
00002940 
000029 50 
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J 



L. 
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QiOQiOOOlO' 



PAGE 



CD 

I 

CO 

CO 



142 



143 



144 
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14? 
149 
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164 
165 

166 



/* */ 

/* THE INTEPNAL SUBBOUTIfslE SIM COMPUTES THE SrHILARlTY COEFFICIENT, */ 

/* FRQiM THE PARAMlETER Y„ THE MATRIX OF CQiEFFIC I EMiTS ARE RETURNED'- */ 

/* IP4 THE S_,COEFF PARftMETeR'. In 40DITI0iNi, THE VALUE QP THE COElFF. */ 

/* IS STORED WITHOUT FORMATTIN'G ONTOi DISK, F'On USE I Isi THE FORTRAN */ 

/♦ STEP. */ 

/* */ 



SIM:PP0CEDUPE 'C Yt S_'C0EFF 1 ; 



DE'C'LflSiE Y '(*,-*,*) FLOAT, /* Y M,ATRIX */ 

PASSEPill :100.Ji REAL FLOAT INITIAL (I'O'EOJf /* Fhr RECORDING */ 
S^C01EFFI*J PLOAT, /* SIMILARITY C'C EFF I C I EiNiT */ 

CHOOSE FLOAT; /* N CHOO:SF 2 */ 

/* CO'MP'UTE SliMlLA'PITY COEFFICIE^LT */ 
S_COEFF=0; 
CH0CSE^HBOUiNlDI¥,2l' ; 

CHQiOi'^.E = lFO'/« «XHOO'SE**2-CHOCSE)/2=0'» t 
GC 1 = 1 TO HBOUNDtS_CQEFFf 1 » B^Y l| 
DO J=l TO HBQUNOIYfZ} BY 1; 

DO K^l TO' tJ-U iBY l; . 

S_C GiEF F J I )' = S^COEFiF ( I J + CHQO'SF / I ABS I LQ^G I Y M ,, J, K 1 J) •♦ IBO » ; 
ENDI 
END'? 

/- RECORD SIMILARITY CDEPFICIENIT ONi DISK */ 
Ji = l; 

5ECD:D0 I^J TO K iBY 1 MHILE I I<'=HBO'UMD « S_CGEFF , l )•)' ; 

PASSiEPI n = S_COEFF'in ; 
S\0 RECOi; 

WJTE FILE 'IHASTERJ FROW |P'ASSER)| 
IF « I < I H BOUiNO' I S_C C E F F , 1 J M THEN 
DO; 

K= 1+100; 
J= I :; 

GO TO RE CO'; 
END; 

ENOi SI-^; ' 



OQ!niO3.02'0' 
00003030 
0'0'003040 
0'Q"003 050 
00'0'0'3060 
0000'3'070 
'001003080 

000'030'^'n 

0000310'0 
00003110 

00'0031,20 
00003130 
00003140 
00003150 
000031 '60 
001003170 

oocc3ian 

000031 '°0 

00003200 

nin0032lO 

00003270 ■ 

00003230 

00003240 

00'0'03250 

0D0fl3,?'6'0 

000032TO 

'0003? SO' 

00003? 90 

00003300 

00003? 10 

0000332 C 

00003 3 30 

00003340 

n'0'003350 

00'0i03160 

0O'0O3?70 

00003 3S0 

00003390 

00003400 

00003410 

00003420 

00003430 

01003440 

000'n.3450 

00003460 



^ '^ im 4. 4. J« im im U'U im ^ UU 



-...,^^^^' #*- '1^^^^ ■■— — *:|^^^^ 



hm 



SFO'*P^CX=0'U«.E OPTIQNS!'Mi&lN| ; 0000^0019 PAGE 11, 

/* */000n37l0 

/* IP_ST4T IS AN INTERNAL SUBPO^UT I iN;^ WHICH CaP'PUTES «• MATRIX D'F */0300?72C 

/* PROBABILITY STATISTICS, ^hi S 'f^ATO'ix QVERWOITES Tl-E PliiRAMETPO */0')003730 

/* ?_V&LUE. • */nor)0^740 

/* */0300?7^0 

^ 3')00?770 

J^ 0000?78n 

O 17^ P STAT:p:RQCEiDUPi tZ_VAiUEl; '" 0000^790 

01003800 
0000? P 10 

175 DECLARE Z_¥ALU€,t*) FLOAT, OUT FLOAT; OOOO^E^O 

oono3e30 

176 DC J=l TC HBQUNDI Z_VALUEflJ PY 1; ' OOOO^rAO 

177 CALL WTP lI_VfiLUE(Jl ,0UIT| •. • ' 000a?e50 
17y Z_VftLUEt JJ^OUT; - 0000^660 
179 cisiD: 0000 36 70 

OOOT^^SO 
190 ='*JiO P_ST6T; " -" 000 J* 690 






S.cQi-piJQiCEO'U^F OPTIQNiSIMaiNJ S OOOO'OOIO PAGE 12 

/* */0'0'00'3910 

/* THE INTERNAL FUrC.TFON NOTR IS A PLl TRflNSLAT I OUsl CF THE FQi^TRaN: */000fl3920 

/* SU8R0UTINIE OF THP S&HE NAMilF, LISTED IN THE IBH SYSTPMi/360i */0000193C 

/* SClfNTIFIC SUBCIO'UTINE PACKAGE 'PUBLICATION. THF CALCULATED iP^.OB- */0000394C 

/* A9ILITV IS RETURMEO TiO THE POINT WHE^E THE PROiCFCUPE IS INVOKED. */000^0'395C 

/* */ 00^0039 60' 

CD • 90003 9 30 

L 000a399'^' 

^ 181 HOTf- I PROCEDURE lEXfP'l; OOOO^C'On 

00004010 

0^00040 20. 

iO'2 OECLABE IEKfDEN,AX,P,T) FLCAT; 00004030 

00004040 
18 3 AX=AESfEX}; 10004050 

184 T = lf 0/1 IE0+0..2316419E0'^AX) ; 00004060 

1P5 'aEN=C..3989423E0*EXPI-EX*EX/2E0) ; 0000407C 

1&6 0' = lEO-DE!N*Tftll I i( 1, 330274E0*T-1„821256E0»*T + 1,.781478E0)*T 0nn04080 

-0.3565638E0}*'T + 0. 319381 5E 1 ; 00004090 

1B7 '■ P = IP-C,5E0}*2E0l 000^041 nO 

•■ 00004110 

la^e E'NiD NDTP: ; 000041,20 

0004130 

00004140 

. 139 ENH' SED; ' 00'004150 
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HAMILTON HARBOUR STUDY 
SECTION C 
DISPERSION 



SUMMAR Y 

A CLUSTER OF THREE DROGUES WERE TRACKED FOR THREE HOURS 
UNDER DIFFERENT WIND CONDITIONS IN THE EASTERN PORTION 
OF THE HARBOUR AT A DEPTH OF 2m BELOW THE WATER 

SURFACE. Measured effective dispersion coefficient 

RANGED FROM 1.7 TO 10^ CM^/SEC. ThE MEASURED EFFECTIVE 

dispersion coefficients as a function of time behaved 
IN the same manner as other Great Lakes studies. The 

EFFECTIVE DISPERSION COEFFICIENT OF 1.2 X lO'^ CM /SEC 
USED IN THE NUMERICAL MODEL AFTER 1974 REPORT FOR 
W HOURS WAS FOUND TO BE REASONABLE. 



HAMILTON HARBOUR STUDY 
DISPERSION (MULTIPLE DROGUES) 



INTRODUCTION 



In the application of the two dimensional model it was 
necessary to use a dispersion eoeffloient. The dispersion 
coefficient used consisted of two parts, namely a constant 
part due to wind stress, and a variable part related to 
local current speed, depth and bottom roughness. Sensi- 
tivity runs were carried out on the model using values of 
and 9.3 X 10' cm^/sec for the wind generated part of the 
coefficient. In general the lower value of wind generated 
dispersion caused numerical instability problems while the 
larger value produced more meaningful results when the model 
was compared to water quality survey results. Dispersion 
coefficients ranged from to^ 1.2 x 10''* cm,^/sec In the 
numerical model. The upper value is within the limits deter- 
mined by dye injection experiments near the lakeshore. It 
Is also compatible with the apparent horizontal dispersion 
coefficients in the range 10** to 10®cm^/sec for length 
scales of 1-10 km and model runs for 48 hours to 25 days, 
used by Simons (1972) for his Lake Ontario model. To ensure 
that the absolute value for u ;.:^^pi.'i-3lon did not adversely 
affect a num,erlcal so'lutio^n thc^ same value for the disper- 
sion coefficient was used in the comparison of alternatives. 
While it is known that the wind generated portion of the 
dispersion coefficient should not be assumed constant 
throughout the numerical grid or constant with respect to 
time it was not practical to evaluate the dispersion coeffi- 
cients at numerous locations In the Harbour simultaneously. 
Consequently these assumptions were necessary. Whether the 
values assumed for the dispersion coefficients were reason- 
able and whether the dispersion coefficient behaves in an 
expected manner as a function of time could be checked In a 
qualitative manner. This work was undertaken in 197^. 

Three similar drogues (2 aluminum sheets 1 m^ at right 
angles suspended from a small surface float) were set at 2 m 
below the surface and released from a common point at the 
same time. The positions of the drogues were determined at 
various times after release. The separation of the drogues 
fro'm, their center of mass (centrold of fixed locations) was 
determined graphically for each time that the drogue loca- 
tions were known. The variance (a-) of these distances was 
then an estimate of the dispersion or the apparent horizontal 
dispersion coefficient e = a^/2t where t = time. The drogue 
release was repeated at well off-shore locations (water 
depth approximately 12 m) at the eastern and western ends of 
the Harbour on different days in July. Wind records during 
the drogue releases are presented in Table I. 
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RESULTS AND DISCUSSION 

Computations of tno :ipparent hor'-lsontal dinper:;iori c.oeV- 

ficientf and vari-ince:3 for c^^c■Y^ nf the nuiin are pr 'Seriled In i 

Table ?^ V'loit^ ^ f the lc>f, vrjr'Jaruje versuL^. lor,, time aiui lop; I 

dispersion coe I'f 1 t:ient yfvsMi: It^;; t,ii7ie are preLK-nted In ^ 

I'Ugures 1 and ? with t[jr ;nan^ /t Icn 1 .1 y ci^t m^uiI rK^(i Lter-t-fH ij 

raraiglit liner.. Tlie rl.ope of t.nc i o^'; varinnctf v^ ]';>.*. 1 oj;; "■"I 

Lime Is ?.U^ wh-'ch is similar te the results oV Okubo arid ^| 

Farlow (196W) who used multiple drogues (30 to ^iO ) at depth: |^ 

uf !.h to 6 m on lakes Erie and Michigan and found a slope 

of 2.0. Well it is obvious i^h:^i. estimating; dispersion 

coefficients by tracking the reparation of three drorues ean 

at best be a crude r-ipproxirnation (subsequently drcir.ue 

experiments used ^1 and h dr'"\p;ne '!lusi,ers), ttie j-esul ts f^r 

tlie variance compare' favourahlv vm t h i^ther more t- k tensi ve 

experiments. And it ai^p-ars that hori/,ontal dispersion a.; 

measured by drogue clusters mi li^rnlllon Harbour i-xhib!t the 

same growtii rate with time as op.en lake experiment.-.. It is 

also notewn-'thy thcjt Okiiho rom:^(^r\lr^ thnt "the f-xnet'ted 

tendency oP the group sJ/,e t^« \,)cr'C.ju- with idiia w.-i.s not 

always noticed". A similar JbuGtiun k' appar-ent Ui this 

llrndted experiment in Table ?. 

An apparent dispersion coefficient of 1.? x 10"cm^/sec is. 
not attained until a time of 50 hours on the fitted line and 
it falls within the 95? confidence intervals at k.H liours 
according to the drogue experiments, HoweveFj the longest 
period of drogue tracking was 2.5 hours s.o that dispersion 
coefficients of 1..? x 10'''cm^/sec is not occuring until well 
beyond the limits of the experiment. Okubo and Parlow^ 
(1967) whose experiments extended for nearly 5 hours, found 
the log variance to grow with time at the more rapid rate of 
3 in the period 1.5 to 5 hours. Consequently one would 
expect 1.2 X 10'*cm^/sec to occur after shorter perlo-ds of 
time approximately I6 hours (or 2 hours for the 95% confi- 
dence intervals). It appears that for model runs of 20 
hours and more in the Hamilton Harbour report ^ (May 197^ ) ^_ a 
horizontal dispersion coefficient of 1.2 x 10** cm' '/sec Is 
no't unreasO'nable . 
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TABLE 1 

HAMILTON HARBOfi 197^ 
MINDS 



Winds. 



Date 



July 5 



July 9 



July 10 



Tim.e 

0600 
0700 
0800 
0900 
1000 
1100 
1200 
1300 

0600 
0700 
0800 
0900 
1000 
1100 
IPOO 
1300 

0800 

0900 

1000 

1100 

1200 

1300 

1^00 

1500. 

160Q 



Dir. 

From 

225 

270 
270 
270 
270 
270 
270 
270 

225 

270 
270 
225 
270 
270 
270 







45 


45 
• 



Speed 


Drogue 


m/sec 


Tracking 


3.57 




H,H6 




5.36 




5-80 


start 


6.211 




6.69 




6.24 


stop 


5.36 




3.. 24 




3.13 




3.13 




3.13 




3 . 57 


start 


^.46 




6.2H 


stop 


6,24 


- 


3.57 




2.68 




3.57 




4.46 




3.57 


start 


2.24 




2.68 


■•! 


3.57 


stop 


4.01 





i 
i 
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TABLE 2 







HAMILTO'N 


HARBOR 1974 








MULTIPLE 


DROQUES 








t 




a'- 








Tlfne 




¥arlanee 


e = (JV2t 


Date 




■ sees. X 10 ' 




cm,^ X 10^ 


cml'^/see 


July 


5 


1.98 
5.58 
8.82 




1.12 

11.3 ■ 
3.2 


283 

1015 

182 


July 


9 


1.71 
3.6 

5.7 




0.07 
1,41 
5.0 


.^60 

1953 
4390 


July 


10 


2.1 
4.86 

8.88 




0,49 
3.16 

6.39 


117 

30' 5 
360 



n" 



c-h 



I# 



I0« 



I 

Ml 



ro* 



no* 



S5% CONFIDEMCt INTERVAL 




LEeEND 

iLO'S CT-^ = (;t..04i LOe t) + O.OTtS 

STANDARD ERROR - t.S 

r a 0.83 



IIQ'8' 



10- 



f TIM'i (SECSl 



10* 



FIG. 1 - MULTIPLE DROSUES JULY IST4 



C-b 



F 



10* 



Z io» 



O 
O 



Q 10 = 



SS% CO'MFUDEiCC INTERVAL^ 




Li JEW 

L9« € - |l.04il LOS tl ♦ '0-S3 
STANDARID iRROR = §-3 
r = 0.f I 



10'' 
10* 



lo- 



10' 



I 0^ 



t TIME IS ECS) 
FIG. 2 - HULTlPLf PiOiUiS JULY IS74 
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HAMILTON HARBOUR STUDY 

SECTION D 

PERIODICITIES AND MASS EXCHANGE 
THROUGH THE BURLINGTON SHIP CANAL 

sumARY 
Periodicities related to semi-diurnal and lake oscillation 

EFFECTS were OBSERVED IN THE RECORDS OF WATER CHEMISTRY 
AND CURRENT METER DATA OBTAINED IN THE CANAL AND HARBOUR 

DURING 197^. These periodicities ranged from 12 hours to 

20 MINUTES AND INFLUENCED THE INTERPRETATION OF WATER 
quality DATA COLLECTED BY CONVENTIONAL SURVEY TECHNIQUES, 

Analysis of current data obtained in the canal for a 
10 day period in the absence of thermal stratification 
indicated a net exchange to the lake of approximately 
10 percent of the harbour volume. a similar exchange rate 
towards the harbour was observed during a period of 
thermal stratification; however/ the presence of a complex 
flow regime frequently including three layers of different 
temperature and conductivity increases the uncertainty of 
this estimate. nevertheless. these estimates indicate the 
obvious importance of these exchange mechanisms in main- 
taining harbour water quality. 
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HAMILTON HARBOUR PERIODICITIES AND MASS 
EXCHANGE THROUGH THE BURLINGTON SHIP CANAL 

INTRODUCTION 

As part of the Hamilton Harbour Study (MOE, 1974), recording 
water chemistry and current meters were operated at two 
depths in the Burlington ship canal during 1972 and 1973. 
The data obtained indicated that water motion in the harbour 
and ship canal is a complex phenomenon. In the ship canal 
non-stratified periods were characterized by alternating 
flows of harbour water to the lake and lake water into the 
harbour with periods corresponding ta semi-diurnal and lake 
oscillation effects. In times of stratification, these 
motions were combined with simultaneous outflow of warm 
harbour water above inflowing lake water and even occasional 
inflow of lake water above harbour water. In addition, high 
frequency variations of periods 0.1 to 0.2 hr, likely due to 
harbour oscillations, have been observed. 

Within the harbour, currents and water quality were also 
found to be influenced by both harbour and lake oscillations. 
These motions produced temporary displacement of the thermo- 
cline of up to several meters, rendering the interpretation 
of grab sample data difficult and necessitating the use of 
numerical averaging techniques in determining stocks of 
nutrients, dissolved oxygen, and other chemical substances. 

In 1974, the instrumentation program was extended to include 
six pairs of recording current and water quality meters. 
The purpose of this program was to provide a synoptic data 
set suitable for calibration and verification of a three- 
dimensional numerical model. In addition, an improved 
estimate of mass exchange and periodicities in the ship 
canal was desired. This section discusses the results of 
this latter phase. 

EXPERIMENTAL PROGRAM 

In 1974, six installations, each consisting of one current 
meter and one water chemistry meter were operated in June 
and July. These installations occupied four distinct 
geographical locations, shown in Figure 1. Two of these 
(the ship canal. Ml and M2, and the offshore location, M6 
and M7) consisted of installations at two different depths. 
As in 1972-73, shore-mounted Schneider robot monitors were 
operated on the ship canal (installations Ml and M2 ) . These 
units measured temperature, conductivity, dissolved oxygen 
and pH, and- were supplied by submerged pumps located 3.6 and 
6.7 meters from the bottom. Recording current meters were 
operated at the same depth and locations. The depth separa- 
tion was less than desired to measure the mass exchange 
between the harbour and the lake in a two layer system, but 
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was dictated by the wave activity In the ship eanal. The 
instrument separation did meaeure the two-layered system In 
the canal on some occasions and wa,s supplemented by depth 
profiling. 

Submerged water quality and current meters were -also used in 
the Ottawa St. Slip (installation M5) and the aoutheast 
corner (Windermere Basin, installation M:8) in order to 
document the variation of water quality from the major 
Industrial discharges and the Hamilton Water Pollution 
Control Plant,, for input to the numerical model. 

In addition, water chem,istry was measured at two depths In 
the eastern" end of the harbour (installations M6 and M7)s 
using NERA recording water quality meters. These Instruments 
consist of submerged sondes employing Hydrolab sensors- which 
mO'nltor the same parameters recorded by the Schneiders ^ as 
well as depth and oxidation-reduction potential. The signals 
are amplified and recorded on a surface-moiunted digital 
recording package. Two such units were supported from.: a 
Bralncon plank-on-edge-buoy . 

RESULTS 

As the methods of data analysis are identical to those used 
in 1972-73 (MOE, 1974), they will not be repeated here. 

Means and standard deviations of water chemistry data 
obtained in 197'4 are summarised in Table 1. Unfortunately,, 
not all meters funGtioned during any one time Interval. The 
intervals for which most meters functioned were June 25 to 
July 5 and July l6 to 22. Almost complete data was obtained 
from the ship canal for both Intervals. On the other hand, 
complete data were obtained only for the first interval from 
Windermere Basin^ and the second interval at the offshore 
location. The meter at the Ottawa St. Slip produced only 2? 
hours' of usable data Immediately after installation. 

A statistical summary of 1974 eurrent meter results is given 
in Table 2,, and current meter measurements are compared with 
drogue trackings, in Table 3. At the Ottawa St. slip, all 
meter speed readings were lower than the drogue readings by 
a factor of 2 to 5- This was due to slime and oil debris 
fouling the I'ctor, despite several cleanings. The time 
average speed is thus expected to be between 10 and 20 
cm/sec. The Windermere Basin record is characterized by 
greater than 50^ negligible ( <0.3 cm/sec) currents. Although 
the results do not correspond exactly to current meter 
reading, the drogue trackings indicates low speeds (less 
Chan 2 cm/sec, except for July 2, when drogue results may 
have been influenced by a 20-2"^ mi/hr wind from the west). 
Drogue comparisons for the upper offshore location (M6) were 
fairly good (speed readings on the meter tended to be 
so:newhat faster than those observed during drogue tracking). 
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No record was obtained from the lower current meter at the 
Bame location (M7). All other current meters performed 
satisfactorily In 197^. These results emphasize the difficulty 
of generating a synoptic (instantaneous spatial water quality 
deecrlptlon) data set by operating a large number of instru- 
ments sdmultaneously . 

Means of the along channel component of current data from 
the ship canal for one month periods varied from 1,5 to 8 
cm/sec; corresponding standard deviations were 2 tO' 11 
cm/sec (1.0 to 1.5 of mean). The large relative standard 
deviations reflected the alternating directions "of flow 
within the canal. 

Figures 2 and 3 show examples of the depth distribution of 
temperature and conductivity measured at different times 
with bathythermograph casts and conductivity meters. Three 
layers of differing temperatures were frequently observed 
during periods of thermal stratification. This could possibly 
represent outflow of two layers of stratified bay water plus 
lake water inflow, or other similar combinations , however, 
with the instrumentation used, a detailed explanation is 
impossible. The complexity of the depth distribution in the 
canal cannot be adequately described with only two recording 
chemistry meters during periods of strong stratification and 
unsettled conditions. Howev'-^r, it v/as possible to use the 
current meter records to estirk!.te narbour-lake exchange 
rates by avoiding very active periods (high winds and veloci- 
ties) and stratified periods when the chemistry meters did 
not sample both the eplllmnion and hypollmnion. 

The significant (95 percent confidence level) periods appear- 
ing In the auto- variance density spectra have been summarized 
In Tables H and 5. The data in Table H are a summary of the 
Epectral densities observed using all data sets obtained in 
1973 and 197^'. Not every data set exhibited all these 
periodicities at the 95 percent confidence level; however, 
these periodicities were present in so'me data sets at the 95 
percent confidence level and in most data sets at the 80 
percent (confidence level. Periods longer than one hour were 
not always resolved for a measuring interval of 10 min, but 
were obtained upon analysis, of the data averaged over 30 m.in 
Intervals. A typical example Illustrating the strong period- 
icity observed in the ship canal appears for dissolved 
oxygen in Figure ^, and for currents in Figure 5- Typical 
varlaiice density plots for dissolved oxygen and currents at 
an offshore harbour location appear In Figures 6 and 7 
respectively. As fewer records are available in the harbour 
for shorter intervals of time, all spectral periods occurring 
at the 80 percent confidence level -irc presented in Table 3 
to facilitate comparison with canal data. 

In the previous report (MOE, 197^), the canal periodicities 
were compared with theoretical lake O'SclllatlO'n perlO'ds 
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calculated by Rockwell (I966). Since then, Rao and Schwab 
(1974), and Freeman, Hamblin, and Murty (197^) presented 
theoretical calculations of periods of Lake Ontario and the 
Helmholtg mode for Hamilton Harbour, The periodicities 
calculated by Rao and Schwab shown In Table 4 compare very 
well with the observed data. Unfortunately, the ■■presence of 
the 2.0 hr period is indicated by only a small number of 
data sets, and does not really confirm Its existence. Some- 
shorter periods are visible; although these may be higher 
longitudinal modes or transverse modes of oscillation. It is 
impossible to verify their nature. 

In addition, a period of 2.5 hr (recalculated as 2 . 3 hr 
using average 1973-74 water levels) has been calculated for 
the Helmoholtz or pumpin.o; mode for Hamilton Harbour (Freeman, 
Hamblln and Murty 197'0. This mode arises from the presence 
of the canal, which i::^ an opening in the harbour with a 
width mjueh less than the harbour dimensions. Such an openJnr; 
Induces an additional gravitational mode of a period much 
longer than the f undametal node . This gravitational mode 
Increases to infinity as the width of the openlnp; is decroa^c^d 
to zero. 

The observed period of 2.4 hr (Tables 4 and 5) may be due to 
either the 'third free oscillation of Lake O'ntario or the 
Helmholtz mode of Hamilton Harbour. However, it Is possible 
that the period of 2.0 hr observed at both depths In the 
offshore harbour location (Table 5)3 but only seen very 
rarely In the 1973 chemistry data in the canal (Table 4), 
may be due to the Helmholtz mode. 

Significant cross correlations between water chemistry and 
currents in the canal and the offshore locatlo^n are given ^ in 
Table 6, and significant cross correlations of water chemistry 
between" upper and lower levels at the above locations are 
presented in Table 7. The presence of a large number of 
signiflGant cross correlations between chemistry and currents 
attests to the common mechanism responsible for their period- 
icities. However, the reason for the lack of many signifi- 
cant coherences between chemistry and currents at the offshore 
location is not known. 

Chemistry is generally not correlated between the two levels 

at the offshore location due to the strong thermal stratifi- 
cation that was present. On the other hand, the complex 
flow regime in the canal, combined with the instrument 
spacing used (Figure 3)* produced significant coherences 
between the two levels in the canal. 

CALCULATION OF MASS EXCHANGE 

The current meter records were examined in order to estimate 
a value for mass exchange in the canal. During some periods 
of stratification the complex nature of the flow regime 
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preveoted the estimation of any meaningful value for mass 
exchange. 

During a period of no thermal stratification (Figure 2) but 
periodic flows (September 26 to October 6, 1973) using" the 
10 minute readings of velocity ^ one finds a net mass flow of 
2.7 X lO"'' kg/sec towards the lake. Mean flows In each 
direction were approximately 10.6 x 10** kg/sec persisting 
for 3.2 K 10^' sec towards the harbour and 5... 3 x 10^' sec 
towards the lake. This represents a mass flow of 3,H x 10-° 
kg (3.4 X 10^ m') of _ water towards the harbour and 5.7 x 
10 kg (5.7 X 10^ m^) of water towards the lake in 10 ■days. 
The net exchange to the lake in these 10 days is approximately 
10 percent of the harbour volume indicating" the importance 
of the ship canal in maintaining harbour water quality. The 
exchange dilutes chemical constituents, in the harbour and 
provides dissolved oxygen for the botto^m sediments and 
chemical oxygen demand discharges. A similar analysis of 
mass exchange for June 13 to July 22, 197^ indicates a mean 
mass flow to the lake of 3-5 x 10** kg/sec and to the harbour 
of 6.1 X 10** kg/sec. The net exchange Is 2.6 x ID'"* kg/ sec 
to the harbour^ or 1 percent per day of the harbour volume. 
However, as this estimate was obtained during a period of 
strong thermal stratification, it is subject "to a larger 
amount of uncertainty than the first estimate. At any rate, 
these e'Stimates indicate the magnitude of exchange which 
does occur through the shli.) "irial . 

The important -observation Is that actual r(3cords i-ire r'eqalred 
in a numerical model for the h.arbour (MOK, 197^, section E) 
If it is to be representative of the processes. VJhlle the 
annual throughput Is a p p r o x 1 m^; t e 1 y 8 p e r c c- n t per y e a t ■ ^ t h e 
periodic exchange may be of the order of 10 percent in LO 
days. The magnitude of these exchanr.e mechanisms makes It 
difficult to statistically detcr-mine reductions or increases 
in harbour stocks of material. An exnriiple of the 1r:ipcrtance 
of the effects Is indicated in Section [■: , 01 s -solved Oxygen 
Uptake where it Is shown "that th^; amount of dlL-.solved oxygen 
entering Hamilton Harbour throafii tiie jijlp cnnal ii^ important 
In satisfying chemical oxygen dej-inni lo:;.d Ln-^^, and L-,ediment 
oxygen demands. 

Considering the importance of the ship canal in maintaining 
harbour water quality, it is important that any future 
development of an outer harbour in Lake Ontario be constructed 
in such a manner that the exchange through the canal will 
not be significantly reduced. Prior to any work being 
undertaken, it is desirable that proposed projects be tested 
in some exchange model and the results be incorporated into 
an environmental impact statement demonstrating that the 
developments will not reduce the harbour -lake exchange. 

The 197^' survey program attempted to generate a synoptic 
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data set at depth by opei-ating six pairs of water quality 
and current meters concurrently. These six pairs of data 
sets would provide a time history of the principal Inputs to 
the model (the two largest waste discharges and the harbour- 
lake exchange), as well as an offshore location in the 
harbour, for model verification purposes- Presently It is 
difficult to obtain a concurrent record at six locations 
with all sensors working in a calibrated condition when 
checked against In situ samples. Until this technology Is 
perfectedlt is unlikely thai; meaningful predictions of 
water chemistry can bo made for the coastal regions. 
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TABLE 1: Means and Standard Deviations of 1974 Water Chemistry in Hamilton Harbour as measured by 



Recording Water Quality 


Meters 




























Location 


Time Period 


# of 
Readings 


Tempi 
Mean 


srature 
SD 


Dissolved 
mg/1 
Mean SD 


Oxygen 

% Sat 
Mean 


• 

SD 


Mean 


pH 


SD 


Conduct 
(\OThos/cm 
Mean 


ivity 
at 25°C) 
SD 


Upper Canal 
{6.7 m from 
bottom) (M) 


June 2 7- July 4 
July 16-22 


1013 
844 


15.7 
12.5 




3.1 

2.3 


7.8 
9.6 




2.3 

1.7 


78 
89 


21 
13 


7.9 

8.0 




0.3 

0.2 


455 

401 


56 

37 


Lower Canal 
(6.7 m from 
bottCTi) (M2) 


June 29-July 12 
July 16-22 


1896 
845 


9,9 
8.9 




3.1 
1.7 


10.2 
11.6 




1.1 
1.1 


90 
99 


10 
9 


8.1 
8.1 




0.2 
0.1 


367 
354 


34 
22 


Upper Offshore 
a (15.7 m from 
^ bottom) (M6) 


July 11-22 


Sll 


22.2 




1.2 


7.6 




1.6 


87 


19 


7.9 




0.3 


518 


9 



Lower Offshore June 25-July 5 
(6.6 m from July 11-22 
bottom) (M7) 



Ottawa St. 
Slip (M5) 

Windermere 
Basin (MS) 



§11 



14,1 
12.6 



0.8 
0.7 



1.3 0,9 12 



7.4 0.1 



June 25-26 



11 28.3 1.0 



5*0 



1,2 64 15 



June 25-July 9 2007 21.0 1,6 4.1 2.0 46 23 



501 
461 



11 
12 



677 10 



603 46 
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TABLE 2 



Summary of Current Meter Results 



PERIOD 



Lo^cation 



JUNE 197^ 



Ml 



M2 



M5 



M6 



M8 



JULY 197^ 



Ml 



M2 



M5 



M6 



M8 



I 



Resultant Speed (cm/sec) 
(Veetor average) 

Resultant Direction covering 
from (degrees) 

Period of Negligible** Speed 
\% of recording time) 

Average Speed (cm/sec) 

Maximum Speed (cm/sec) 

Persistence 

Total No, of readings 

Interval between readings 
(minutes) 



1.48 



82 



8.13 

43 

0.18 

2448 
10 



6.70 6.41 2.89 



10.28 

41 

0.65 

3102 

10 



6.87 

27 

0.93 

4874 

5 



3.77 

13 

0.77 

4888 

5 



0.97 



74 190 159 166 



8* 26* 60 



1.04 

9 

0.93 

2304 

10 



6.66 10.85 4.46 1.77 1.09 



66 



7.80 

40 

0.85 

2448 

10 



33 190 



12.54 

41 

0.87 

310 3 

10 



35* 

4.69 

18 

0.95 

5057 

5 



16 1 
19^ 

3.57 

13 

0.50 

6247 
5 



169 
48 

1.39 
11 

0.79 

3099 

10 



* <2,0 cm/sec 
** < 0.30 cm/sec 



TftBLE 3 

DROGUl MIASURE'MINTS 
HAMILTON HARBOUR, 1974 



Location 



Date 



Time 



Drogue:_ 



Meter 











Direction 


Speed 


Diirection 


Speed 










(going to) 


Ccm/eec) 


(going to) 


(cm/sec) 


Ottawa St. 


J'une 


26 


15.20 


20° 


10-12 


15-35° 


1.9 


Slip (M5) 


Jime 


27 


12.15 


20° 


10-11 


28° 


2,7 




July 


2 


11.45 


25° 


11-12 


11° 


4.2 




July 


9 


11.00 


25° 


11-12 


10° 


4.7 




July 


10 


13.30 


29° 


6-10 


5-30 


1.6 




July 


16 


11.00 


29° 


8-9 


15-30° 


5.0 


Upper offshore 


June 


26 


13.45 


342° 


3-5 


327° 


7.3 


location (M6) 


June 


27 


11.10 


337° 


3-5 


325° 


2.2 




July 


9 


11.55 


15° 


3 


30° 

330 


4.7 




July 


10 


14.30 


332° 


4 


6.7 




July 


16 


12.00 


115° 


0.5 


3.0 


Windermere 


Jujie 


26 


15.00 


iri 





308° 


0.4 


Basin (MB) 


June 


27 


12.00 


357° 


0.7 


10° 

Oi 


0.2 




July 


2 


12.15 


106° 


2-6 


200-205 


1.0 




July 


9 


11,30 


49° 


a 


- 


0, 




July 


10 


13.45 


340° 


1 


** 


■"Qn 
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TABLE 4 

THEOHETICAL AND OBSERVED (1973-4) 
PERIODICITIES IN HOORS IN THE 
BURLINGTON SHIP CANAL (95% confidence) 



OBSERVED 



Theoretical 
Lake OntariO' 
(Rao & Sehwab 
1974) 



5.11 
3.11 
2.31 
1,87 
1,78 
1,46 



Chemistry 
(Dissolved oxygen.* pH 

and conductivity) 



Upper 



Lower 



16 

4,6 to 5.3 
3 . 2 

2.2 to 2,4 

1.95 to 2,0' 

1.65 to 1.75 

1.4 to 1..45 

1.0 
0,80 to 0.85 



5.0 
3.2 

2,3 



12 

to 5.3 

to 3.4 
to 2,5 



1.65 to 1.75 

1.4 to 1.45 
1.25 

1.0 
0.75 to 0.85 



Currents (along canal) 



Upper 



12 
5.0 
3.2 
2.4 

1.7 
1,35 tO' 1.4 
1.25 
1.0 
0.86 
0.67 
0,56 



0.35* 



0.65 to 0.70 0.63 to 0.67 

0.53 to 0.59 0.52 to 0.57 

0,43 to 0.47 0.44 to 0.46 0.43 to 0.47 

0,38 to 0.39 0.38 to 0.41 0.37 

0.35 to 0.36 0.35 to 0.38 0.35 



Lower 



12 

5.0 
3.2 

2.2 to 2.4 

1.7 
1.35 to 1.45 

1.25 

1.0 

0.84 

0.67 

0.55 

0.47 
0.37 to 0.3B 

0.35 



NOTE:* Predicted first mode for H^amilton Harbour, using Merian 
equation. Periods less th.an 0.5 hr. , were observed by 
.analyzing unsiioothed data. 

Not all periods were observed at the 95% confidence level 

in any one data set; however, spectral density concentrations 

at the 80% level or lower were generally observied in each 

data set for most periods greater than 3 times the measur- 
ing interval . 






D-ll 



I 



Theoretical 
Lake O'ntario 

(Rao and 
Schwab, 1974} 

5.11 
3.11 
2.31 
1.87 
1.78 
1.46 



Ottawa 
St. Slip 

(MS) 
Currents 

5.5 
3.2 

2.4 to 2.7 



TABLE 5 



OBSERVED PERIDOICITIES IN HOURS 

HAMILTON BARBOUR, 1974 (80% Confidence) 



Offshore Locatioin 



Upper (Mg) 



Currents 



Chemistry 



Lower CM7I 
Chemistry 



4.8 to 5.3 4.8 to 5.3 

3.4 3.2 to 3.4 

2,4 to 2.8 2.5 to 2.7 

1.9 to 2.0' 

« - 1.8 to 1.85 

1.3 to 1.45 1.4 to 1..55 
1.25 

1.1 to 1.15 1.1 to 1.15 



Note: Current iieter peaks are a composite of peaks obse.rved 
-in the nor'th-south and east-west directions. 



Windermere Basin (M8) 



Currents 



2 . 6 
2 . 
1.7 

1.38 
1.18 
1.08 



Chemistry 



3 . 



1.6 5 to 1.8 



1.11 
0.81 
0.68 

0.50 
0.44 



"^^ !■ Ini' " "Mt^~ im^ mm 



I ^_ J 



1 i^J 



J 



TABLE 6 : Significant Cross Correlations between Water Chemistry and Currents, 
Hamilton Harbour, 1974. 



Location 



Date 



Temperature Dissolved Oxygen pH 



Conductivity # of Data Points 



Burlington Ship Canal June 27- 
Upper (6.7 m from July 4 
bottom) (Ml) 



1.2-10. 

*3.3(25) 
*2.5(18) 
*1.67(40) 
*1.43(1) 



1.25-5.0 
*4.0(-157) 
*2. 5(173) 
*1.67(-165) 



0.95-1.00 (-4) 1.25(166) 
0.83-0.87 (-24) 1.00(103) 
0.48(-153) 



2.5-5.0 
*5,0(-161) 
*2.9(-171) 

1.67(174) 



1.2-2.0 
*5.0(15) 
*3.3(8) 
*2.2(-l) 
*1.67(7) 
* 1.43 (-27) 
*1.25(-32) 
1.00 (-73) 
0.83(-82) 



1019 



i 



Burlington Ship Canal 
Upper (6.7 m from 
bottom) (Ml) 



July 16-22 



Burlington Ship Canal 
lower (3.6 ro from 
bottom (M2) 



June 
July 



29- 
12 



2.2-9. 
*4.5(9) 
*3.0(19) 

*2.2(14) 
1.78 (-18) 
1.45(5) 

1.2-6.0 
*4,0(41) 
*2.2(16) 
*1.71(2) 
*1.41(10) 

0.86 (-29) 



2,7-13. 
*13(-154) 
*5.3(-173) 

*3.0(-168) 



1. -12. 
*4.8(-146) 
*2.2(-176) 
*1.71(-171) 

1.04(143) 
0.86(142) 



none 



4.8(-134) 
2.4(173) 



2.0-13. 
*4.4(10) 
*3.0(18) 

*2.2(-17) 

1-78 (-29) 



1.3-12. 
*12(1) 
*4.0(19) 
*2.4(2) 
*1.71(3) 
*1.41(-24) 
0.98(-38) 
0,83(-41) 



1896 



Burlington Ship Canal 
lower (6.7 m from 
bottom) (M2) 



July 16-22 



3.0-5.3 
*4.5(14) 
*3.3(19) 



4.5(-170) 



4.5(-166) 



4.5(16) 



^5 



Upper Offshore Location July 11-22 
(15.7 m from bottom (M6) 



none 



none 



14. (-151) 

1.13(164) 



1,71(172) 



513 



NOTES: Along channel component is used for ship canal locations. North component is used for 

offshore location. Indicated figures are the periods in hours for 5 percent significant 
coherences (Panofsky and Brier 1968) , followed in parentheses by the phase angle in degrees. 



♦Asterisks indicate periods of maximum coherence within the period (stated immediately above) 
in which coherence is significant throughout. 



TABLE 7 : Significant Cross Correlations of 1974 Water Chemistry between Opper and 

Lower Levels, Burlington Ship Canal and Offshore Location , Hamilton Harbour 



I 



Location 


Date 


Temperature 


Ship 


June 29- 


1.18-6.7 


Canal 


July 5 


*2.86(-6) 
* 1.67 (-9) 
*1.43(14) 

1.00 

0.80(62) 


Ship 


July 16-22 


3.0-6.0 


Canal 




*4,8(22) 



Dissolved Oxygen 



Ifflf 



Conductivity # of Data Points 



2.2-6.7 
* 5,0 (-15) 

1.25-1.82 
*1.43(-21) 



4.0-6.0 
*4.8(20) 



2.2-3.3 

*2.9(0) 



4.8-'6.0(-2) 



1.25-6.7 
*2.9(-2) 
* 1.67 (-10) 
*1.33(-20) 



3.0-6.0 
*4.8(23) 



755 



845 



Offshore July 11-22 none 



4. 8 (-35) 



none 



none 



wm 



Note: Indicated figures are the periods in hours for 5 percent significant coherences 
(Panof sky and Brier 1968) followed by the phase angle in degrees . 



* Asterisks indicate periods of maximum coherence within the period 

(stated immediately above) in which coherence is significant throughout. 
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CONVENT I OM A L SAMPLING STATIONS 

® METER LOCATIONS 

1-8 WATER DEPTHS IN METRES 



FIGURED HAMILTON HARBOUR SAMPLING AND METER LOCATIONS 
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FIGURE ^2 DEPTH PROFILES OF TEMPERATURE AND CONDUCTIVITY 
BURLINGTON SHIP CANAL, JULY II AND OCTOBER 1,1973 
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FIGURE -3 DEPTH PROFILES OF TEMPERATURE BURLINGTON SHIP CANAL. 

JULY 16 AND 22, 1974 
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FIQURE4^AUT0 VARIANCE DENSITY SPECTRA FOR 
DISSOLVED OXYGEN CMG/Ll SHIP CANAL STm 
ABOVE BOTTOM FOR JUNE 25 - JULY 5, 1974 
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FIGURE 5 = VARIANCE DENSITY SPECTRA FOR 
SHIP CANAL CURRENTS JULY, 1974 
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U FIGUii e - VARIANCE DiNSITY SPECTRA FOR DISSOLVED OXYQEN 

(M6/L) LOCATION 117, JULY 11-2.2/74 
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FIGURE. 7 - VARIANCE D'ENSITY SPECTRA FO'R CURRENTS 
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HA^iILTON HARBOUR STUDY 

SECTION E 

SEDiriENT DISSOLVED OXYGEN UPTAKE 

SUM MARY , 



A REVIEW OF LrTERATURE DATA ON IN SITU AND LABORATORY 
METHODS OF MEASURING SEDIMENT OXYGEN UPTAKE IS PRE- 
SENTED. EFFECTS ON OXYGEN UPTAKE RATE OF ORGANIC MAT- 
TER CONTENT. SEDIMENT DEPTH. SEASON AT WHICH SEDIMENT 
WAS COLLECTED . BACTER I AL COUNT . MACRO I N VERTEBRATE 
POPULATION/ AND DISSOLVED OXYGEN CONCENTRATION OF THE 
OVERLYING WATER ARE DISCUSSED. 

THE DISSOLVED OXYGEN DEPLETION CURVE WITH TIME IN 1972 
IS DISCUSSED IN TERMS OF UPTAKE RESULTS MEASURED ELSE- 
WHERE i FOR A 109 DAY PERIOD. AN ESTIMATED VALUE OF 
OXYGEN UPTAKE WAS 3 X 10^ KG. AS COMPARED TO AN 
OBSERVED DEFICIT (COMPUTED FROM OBSERVED AND SATURA- 
TION DISSOLVED OXYGEN STOCKS) OF 1.5 X 10^ KG. A 
MAJOR REASON FOR THE DIFFERENCE BETWEEN THESE FIGURES 
IS BELIEVED TO BE MASS EXCHANGE THROUGH THE BURLINGTON 
SHIP CANAL. ESTIMATES OF THE MAGNITUDE OF DISSOLVED 
OXYGEN SUPPLIED TO THE HARBOUR RANGE FROM 4 X 10^ KG. 

DAILY DISCHARGES OF BOD AND COD TO THE HARBOUR ARE 
ABOUT THREE TIMES THE ESTIMATED SEDIMENT OXYGEN UP- 
TAKE. THIS DISCREPENCY IS DISCUSSED IN TERMS OF 
EPLIMNETIC REAERATION AND PHOTOSYNTHETIC OXYGENATION. 



SEDIMENT OXYGEN UPTAKE IN 
HAMILTON HARBOUR 

A.; part of the 1972 Hamilton. Harbour ntudy (MOE, 197^), It 
waj found that dlr-sol/ed ■■)xv^^en levelL^ below S mg/1 existed 
in thp hypo I Imr J Ion from July to Sept ember 5 and in many 
plac Gr>, anoxic conditions were ap)proached. In July 197^-U 
averaf';e dissolved oxygen concentrations recorded in the 
liypolimnion of the east end of the harbour remained below 
2.^3 mp,/l for 11 days, and were frequently below 1.0 mg/1 
(Section A). It was estimated that about tv/o-thi.rds of the 
deficit was caused by discharp;es and sediment demand and 
one-third by recent phytoplankton decay in Auf^ust. In the 
repoT't it was recom.mended that chemical oxygen demand (COD) 
loadings be reduced by approximately 6l) percent, in order to 
reduce oxygen deficits. 

Ai' tlu^ oxyr^en demand of the sediments is an important factor 
In determining what Improvement in dissolved oxygen deficit 
mig,ht be expected from reductions in CC^D loadings, this 
factor should be measured to clarify the effects of existing 
bottom sediments. I'hls section begins with a review of 
existing literature data on sediment oxygen demand. The 
situation in Hamilton Harbour is then discussed with respect 
to the literature data, and suggestions for experimental 
work are made . 

Fvquipment used In O'xygen uptake Studies 

Three basic types of experiments have been used in oxygen _ 
uptake measurem,ents . These are: (a) In situ studies, (b) 
laboratory 'studies on quiescent system,s, and (c) laboratory 
studies using well agitated samples in the Warburg apparatus. 

In Situ Studies 

.^^teln and Denn i son (19b7 ) , and larer, Edberi-; and Hofsten 
(197^,) us.cd cylindical plexlgaiss units for in situ oxyp^en 
ufitakf^ s.tuciies. Tb.e cylintier used by Kdberg and hofsten was 
][]{) mm In diameter an^l .-h-'O mm in heiglit, and bore an outer 
flange -'10 mm from the \owev edge of the cylinder. This 
fLang,e rests e^n tlie mud surface, and the entire unit is 
anchorc>d hy weights or steel rods. Tlie r;xy'^en uptake is 
measured with an oxygen elect, T-ode ^Yellow Springs Instruments 
hr-ckman), with a stirring motor being employed to provide an 
adioquato water ^J.ow past the electrode without stir'rlng up 
the sediment. Although Stein and Dennison used divers to 
deploy their ec^uipment and remove any trapped surface water 
by fcMnning, Edberg and Hofsten employed divers only in the 
case of a rough bottom. The latter authors used a movable 
lid which closes on the bottom to ensure that representative 
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bottom water was trap'ped. A similar apparatus was also used 
'by Hunter et al (1973). 

Ag a portion of project Hypo, an intensive survey of the 
Lake Erie central basin hypol imnion conducted in 1^)70, laicas, 
3.rid Thomas '"197.-'), as well as Blanton and Winklhofor (l'j7.:-'), 
neasured Lake Erie sediment 0Kyf7;en demand with in situ 
plexiglass chambers. The former authors used a trianr:ular 
pr Lsm-shapec3 chamber vvfiich covered 0.2'3 ni^ of bott^nn -iwd 
enclosed 13 £ of water^. A submersible fMimp cirtnilatod tfie 
water with minimum sedirent di si:urbar^ce . T!ie latter autht:)!^;; 
eiiployed ;i box ivhich i-^nc io.-^ed J ^7 ^ of wateT=, coveiM rif'; 0.rM)b 
m'' of bottom. Both typc^s of chamber were equipped witti J'ln:> 
to prevent them from s, inklnp: into the soft Lake Kr'\c s^ediment 
a.jid In both cases, dissolved oxyp;en and temperatur-e were 
rr' o n i t o r e d c o n t .! n u o u s J_ y by r e c o r d in;'; met e r* s^ y s t c m ; : . 

James (197^) used two types of in situ resplrometers to 

measure sediment oxygen demand. The first consisted O'f a 
metal box which was very similar tO' the chambers already 
described. The other was a plastic tunnel 30 m long by 1 m 
wide which could be ancho^red in a stream bed. Oxyr,en deple- 
tion was measured by comparing the dissolved oxygen concen- 
tration at the entrance and exit ends of the tunnel. 

Laboratory Studies - Q,ulescent Samples 

.Several types oC appar^atus have b.-'oii used in the laboratory 
to attempt to simulate ciulescent natural conditions, and at 
rdie same time study the effect:", '"if variables such as tempera- 
tures, Initial oxygen concentrations and the presence of 
ot'ier ions such as nitrate on t\ie v^^xygen uptake rate. The 
u'l situ studies of Stein and Den 1 son (ld67) and fildberg and 
Ik'fsten (IQ'73) were acct^mps.ni ed s,v laboratory investigations. 

Hanes and Irwin (1969) used a Pyrex r'eactor 8.9 cm deep and 
16. 5 cm inside diameter covered with a plexiglass lid 
through which a thermistor and oxygen electrode were inserted. 
Thio water wa,s circulated without stirring up the mud by 
!iioans of a magnetic stirrer placed in a petri dish at the 
bottom of the reactor. A similar setup, employing wooden 
aquaria, v/as used by Stein -and iienlson (1967); they determin- 
ed oxygen by the Winklei' method on siphoned water samples. 

In their laboratory studies , Edberg and Hofsten (197 3) used 
plexiglass tubes 500' mm lo^ng and 70 mm in diam,eter, equipped 
with a magnetic stirrer and a Yellow Springs oxygen, electrode . 
Similar tubes 19 or 29 cm In length were used in a series of 
investigatlO'ns at the Water Pollution Research Laboratory in 
England (Knowles, Edwards and Briggs, 1962; Edwards and 
Rolley, 1965; Rolley and Owens, 1967)- These Invest IgatlO'ns 
i^'eacted six core samples in a cluster o^f tubes mounted, on a 
rotating platfo^rm. Oxygen was measured with a, dropping 
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mercury electrode which was inserted into each tube in turn; 
the water In each tube was stirred by a magnetic stirring 
bar mO'Unted on a perfo^rated disk and driven by a rotating 
magnet mounted in the center O'f the tube cluster. 

Jenkln tubes were us^ed by Hayes, and McAuley il'^^59) . These 
consist of a glass tube equlppea with a spring loaded brass 
disk and rubber pad at the bottom, and are used with a 
spec:ial apparatus (Mortimer, 19^1^'?) for In situ sampling of 
mud cores. Oxygen was determined by the Winkler method. 
Some screw capped bottles were also used by these authors. 

Sedlm,ent samples were collected either dlrectlj,'' In the 
reaction tubes, or by use of an Ekman dredge. In the case 
of dredge samples, the reaction tube was loaded with sediment 
by plunging directly into the dredge. Lake water, was then 
carefully Introduced intO' the tube, in order not to disturb 
the sediment, and aerated for a period of time before begin- 
ning the experiment. 

Laboratory Studies-Agitated Samples 

Several authors (Liu, 197 3; Gardner and Lee; I965) have used 
a manom,etric technique using the Warburg resplrometer , 
Sediment and water are' shaken in 125 ml flasks. Carbon 
dioxide is absorbed using 10^ KOH solutions and oxygen 
uptake is measured m;anom,etrlc.i M y . This apparatus Is well 
tuiited for measuring maximum uptai;.:^ ratss, with results 
Indicating: a factor of approximately ten tim,es more oxygen 
absorbed wltli shaking co'mpared to without. However, for our 
purpO'Ses, this measurement is not desired; thus these tests 
will not be further discussed.. 

Results of Oxygen Uptake Studies 



Some results obtained using. ' ■i''^ fMuipmetit .'j I ready described 
are given .in Table 1. For o /r o^ comparison, all rates 
have been expressed la g0 2/''!r[^'-hr' ■-ina rcca1 cul -"ited to lb C 
using the relation (Edwards and Hoi ley, 196'j) 



t = 
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Ci 'ice tjie oxyron upt. livc-^r. at l e-^mperatures T2 
r, tlie temperature' ct'effMiMont (asr.umed = Qi.O?). 
s terTU)e?Mture coefficient correr^pondr. well to that measured 
Edberg and ([ofsten (197-!) r-Ptweeri 10 and ?0'^C\ however, 
y found a larger coef f Ic i ent at lowr^T' temperatures. In 
w of tlic difficulty in obtaining accurnte oxygen uptake 
ults as well as tiie changes that may occur prior to 
forming laboratory experiments, this temperature coeffi- 
nt can be regarded as a satisfactory approximation over 
experimental range of t^) 
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Oxygen uptake results have been compared according: tO' a 

number of factO'Ps. Among these are organic (o^r volatile) 

matter, season of collectionj, bacteriological counts, s^edlment 
depths presence of ions such as nitrate, addition of macroin- 

vertebrateSj etc. 

(a) Organic Matter 

Most authors have attempted to relate oxygen uptake ratns to 
organic matter (or carbon) content of sediments. Most 
autl'iors have measured organic matter content as loss on 
Ignition of dried sediment at 500 - 600 C. Rolley and Owens 
(196 7) and Edwards and Rolley (1965) measured organic carbon 
by converting CaCOa to CaSOs by a 2k hour pretreatment witli 
SO2 gas 5 followed by oxidation of organic carbon to CO 2 at 
600°C. Organic carbon measured in this fashion correlated 
with loss in ignition with a correlation coefficient of 
0,923, indicating a nighly significant relationship. 

In general J oxygen uotakc increases with increasing organic 
content of sediments • This can be illustrated by the data 
obtained on Swedish lakes (except as indicated) by Edberg 
and Hofsten (1973), summarized as follows: 

Oxygen Uptake Rate 

(g/m^-hr at 15 C) 

O1.02 

0*0^1 - 0.15 

0,1 

Baltic Sea, 11.^1 - 17-7 - ' O.OM - 0.13 

Pulpmill 

effluent site 

Eutrophlc Lake 30 - -^0 ^_ 0.15 - O.I8 

Although no significant difference occurs between the oxygen 
uptakes of the samples with 8 to I8I organic matter, a 
marked effect of organic content is Indicated by the samples 
with l."^ and 30-40^ organic matter. For celluoslc fiber 
sludge deposits with an organic matter content of 80-951 in 
the top 1-2 cm, and 10 5^' beneath the top layer. Stein and 
Denison (1967) obtained an uptake rate at 15' of 0.20 g/m^- 
hr . Due to the extreme gradient In o^rganlc co'ntentj this 
result is difficult' to assess, but it is obviously suggestive 
of an increased uptake caused by the high s^urface organic 
co'ntent. 

In reasonable agreem,ent with Edberg and Hofsten , Rolley and 

Owens (1967) obtained oxygen uptakes of 0.05 - O.O6 g/m^-hr 
for river muds of average 111 'loss on ignition, and 0.10- 
O'.l^l g/m^'-hr for river muds (affected by sewage) of average 
15-21'% loss on Ignition. Edwards and Rolley (1965) present 



K-H 





Organic Matter 


Source 


(% dry weight) 


Lake 


1.3 


Rivers 


8-11 


Lake 


12.2 



a graph in which oxygen uptake Is clearly independent of 
organic carbon over a range of H to 13I C. This represents 
an approximate range of organic matter of 8 to 261,, somewhat 
greater than the similar range of no significant difference 
evident from- the results of Edberg and Hofsten (1973). 

Recently, Hunter (1973) reported a comprehensive study of 
factors determining oxygen uptake rates of benthal stream 

deposits. The organic content of the muds was determined In 
term,s of total organic carbon^ COD^ BOD (5 day and 20 day), 
hexane extractable, hexane-benzene extractable, and total 
volatile solids. Oxygen uptake rates were 0'.05 to 0.5 g/m^ = 
hr and total volatile solids ranged from, 5 to 15^. The 
significance of the various param,eters in determining 
benthal oxygen demand was explored by multiple regression. 
However, m,ost of the coefficients calculated were unexplaln- 
able; some were even negative, and all eight together 
explained less t,han 701 of the variation. Even "the variable 
with the highest degree of correlation (000) explained only 
37^ of the variations. 

These results show that although large variations of organic 
content exert an influence on oxygen uptake rates, other 
factors play important roles In producing the observed large 
variations in uptake rate. 

(b) Sedl,ment depth 

The effect of sediment depth on oxygen uptake, as m,easured 

in the laboratory, can provide an Indication as to how deep 
beneath the sur,face oxygen can be taken up in a natural 
system in the absence of scouring action. ' 

In the laboratory portion of their study In celluloslc fiber 
sludge. Stein and .Denison (1967) found no effect of sediment 
depth on oxygen uptake using sediment depths of 2,5 and 5.1 
cm,. Likewise, ,Edwards and Rolley (1965) fou,nd no effect" of 
depth on oxygen uptake in sediment depths between 2.1 and 
10.2 cm.. 

Prom experiments with sewage sludges. Baity (I938), and 
Pair, Moore and Thomas (194l) suggested that oxygen uptake 
was approximately proportional to the square root of the 
depth for depths between 0.1 and 10 cm. The former authors 
combined the effect of organic matter concentratlo,ns in 
using t,he areal concentration of organic matter as a parameter 
,Although the organic content of their sludges was similar to 
t,hat used by Stein and Denison, the sludges of the former 
authors were settled only briefly prior to the experiment, 
and moreover were subject to considerable agitation due to 
anaerobic decomposition and consequent release of gas bubbles. 
On the other hand, the river and lake sediments used by most 
aut,hors have been accumulated over long periods of settling; 
deeper mud tends to exert less oxygen demand as Indicated by 



E-5 



some "dummy Jenkins"' experiments (liayeB and McAuley, ;1959), 
and Is more resistant to bacterial decomposition than surface 
mud (Aliens, Grendley and Bro^okSj cited by EdwardB. and Rolley, 
1965). 

In a situation like Hamilton Harbour, tlie effect of bottom 
scour, which Is prevented by the apparatus in any In situ 
study J would probably be greater than that of mud depth. On 
the other hand, depth tests might prove interest Ing if 

laboratory work were undertaken. 

(c) Season at which sediment was collected 

The effect of seasons has only been Investigated in laboratory 
studies (Rolley and Owens, 1967; Edwards and Rolley, 1965). 
Although the former authors found that oxygen uptake in the 
summer was only very slif:':htly greater than that in the 

winter, the latter investigators found a definite seasonal 
effect with maximum uptake obtained from samples taken in 
the spring and summer and minimum uptake from winter samples. 
ni^^e reason foT' the difference is suggested to the seasonal 
cdianges In ^iopul;ition of nilcro and macro invertebrates. 

(d) Bacterial count 

While Hayes and MacAuley (19b9) reported a linear relation- 
ship between bacterial count and oxygen consumption in lake 
muds, Edwards and Rolley (19*^^.0 found no dependence of 
oxygen uptake on bacterial count. Possibly other organisms 
or' chemical substances might be important in the sediments 
studied by Edwards and Rolley. Perhaps the effect of mud 
sterilization siiould be investigated in a laboratory study; 
howe\^er, the process of ster J 1 ization would nrobably result 
in chemical changes In the mud. Perhaps potassium cynlde or 
formalin could be used to Inliibit biological activity and 
differentiate chemical from biological oxidation (Liu, 
1973) . 

(e) Macroinvertebrate population 

As a result of tbe lack o^ significant correlation between 
sediment chemistry and oxygen uptake. Hunter et al (197 3) 
extended the study tD include the effect of macroinvertebrate 
population. A very significant relationship between oxygen 
demand and tub! fie id (tub if ex) worm population was obtained, 
with the oxygen demand Increasing from 0.07 to . ^s 8 g/m -hr 
as the tubi'fex worm population Increased from '1300 to 28^,000/m 
After allowing a correction for oxygen demand of sediments 
devoid of m.acroinvertebrates , a demand per worm was calculated 
'IS 0.032 mg/worrn/day , As this was about five times the res- 
piration rate of the wc^i'ms , other factors such as their 
burrowing action bringing oT'ganlcs to the surface may be 
important. In addition, the se<llment contained 1700 to 
13^,000 bloodworms/m^ J which were not considered in the 
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calGulatlO'D but should be significant.. Although their 
er.timatiO'n of a demand per wo'rm may thus be IncO'rrect, their 
resultn illustrate the r>lp.nificant effect of th,e macro;] nvert- 
rhrate cornimunlty in tlotorml nlrif; oxyf;en uptake, ^ 

The 197i? bottom fauna survey In HamlltO'n Harbour shO'Wed the 

presence of more than 11 ,,000 ollgochaetes (tublf lcidae)/i7i^- 
throughout mO'St of the harbour ^ with very little area contain- 

■ing less than IjlOO/m^. Considering the effect oT these 
organism,s on sediment oxygen uptake as Illustrated by Hunter 
et al (1973) J, any experimental pro'gram must be accompanied . 
by a bottom, fauna enumeration. If laboratory studies are 

done the uptake of sediments containing macroinvertebrates 
shO'Uld be compared with that of sterilized sediments, 

(f) O'ther variables 

Edwards and Rolley (1965) studied the effect of nitrate con- 
centration in the ovei^lying water on o>;ygen uptake. They 
found that increased nil rate concentration up to 20 mg/1 N 
had no effect on oxygen uptake. 

rieveral iriver,t tgat I on.-> con/, i dt^rtnJ ttie depondencf^ of uptake 
p.-ite. on oxyg,en concentrn t i.ori , unln>^' an equation oV t,he ft:rm 
R=aC ^ where 1\ in, uptake rate, '" i r. oxygen concentrate on ^ and 
a and b are empirical conntantL: . i-.dwards and Rolley (19'^">5) 
found that the oxygon uptake: i■'^a Pit thl.r. equation; Edberg 
and Hofsten (1973) stated their d ix.a oDtaaned from laboratory 
measurements did not fit th:is equation, although they present 
a graph which should apparently fit It. Theli' in situ data, 
obtained for a shorter time at DO concentrations above 7 
mg/1 5 on the other hand, give a rate independent of oxygen 
co^ncentration. 

Lucas and Thomao (1972) also presented borne data which show 
an effect of oxygen concentrat ien on uptake rate. Oxygen 
uptake rates, corrected to 1S"C, ciecreased from 0.08-0.14 
g/m^-hr at DO concentrations r^f 7-9 mg/l, to 0.05-0.06 g/m^- 
hr at DO concentrations of 3 mr,/l. A further decrease in DO 
content, however, did not result in -my further drop in the 
oxygen uptake raLe. 

The effect of sunligjit c-^u.-jn^' i^h.oto .ynthet Lc oxygen production 
was noted by Lucas and Thomas (1972) in Lakt? Krie. They 
meas^ured oxygen uptake rales in s. Itu usinr; clear and black 
plexiglass clKimbers in Aur;ust r-)7 a,t various times of the 
day. Although tne rate in dark chambers was constant at 
0.05-0.06 g/m^-hr (15^0), the rate in light chambers varied 
dlurnally, being similar to uaat in dark chambers at night 
but as lov; as zero during the da,y1lg:it hours. At that time 
of the year, there was sufficient sunlight penetration to 
the lake bottom to permit algae ro f^j'ow there, or allow 
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surface algae to remain viable after settling tO' the bO'ttom, 

The cO'nslderably lO'wer uptake rate measured by Blanton and 
Mlnklhofer (1972 ), also in Lake Erie as a part of Project 
HypO', must have been Influenced by algal photosynthesis of 
oxygen. Although algae were stated to be p,rowlnK only 
outside the box, microscopic algae may have been growing _ 
within the box. Blanton and Winkelhofer^ in stating that 
their rate is comparable to the "Integrated daily rate" of 
Lucas and Thomas (which includes photosynthetlc effects) s 
admit that their rate was influenced by photosynthesis. 
Certainly, these results indicate how widely variable the_ 
observed oxygen uptake rates can be under apparently similar 
conditions . 

Comparison of in situ and laboratory results 

While Stejn and Denlson (I967) found a hi gher ^ oxygen uptake 
In laboratory experiments as compared to in situ results, 
Edberg and Ilofsten (1973) found a lower uptake in the laboratory 
Results of the lalt-.-r, when corrected for temperature, 
tended to be lower by a factor of as much as two or throe 
than the in situ results. As well, their laboratory results 
Increased by a factor of about two if the sediments^ were 
stored under aerated water, T})eir explanat:ion of this 
difference is that respiring organisms multiply or become 
better adapted to the new envl rr^nmc^nt during storage, or 
that burrowing organisms may increase active surface area. 
Certain! V, removal o? the r.odlment has a decided effect on 
Its biological activity, and this effect could be different 
rvir sediments of d 1 :'f'eT'f>nt nnturr-, as is the case in the two 
in situ 1 n ve s t i ga 1 1 o n s . 

Th.e obr.ervod differences in Lake Krie sediment oxygen uptake 
results: (Li:cas and Thomas, 1972; Rlanton and Winklhofer, 
1972) may also be viewed as partly due to modification of 
the natural conditions produced ny Inserting closed boxes of 
different ^rimensions. Although it Is true that gelatinous 
l^y-rs or mlcrozones may be disturbed during deployment of 
in' situ eouipment, it is felt that the In situ procedure is 
the most direct and .■•iccurat- meth.od of measuring oxygen 
c o n s u m p t i o n s o f s c d 1 rri e n t s . 

Sediment Oxygen Uptake in Hamilton Harbour - - ' 

In the 1972 study of the oxygen budget in Hamilton Harbour 
' (MOE, 1974), dissolved oxygen concentatlons were measured at 
eight aampllng stations twice in the spring, and then weekly 
from raid-July to mid-November. Locations of sampling stations 
used are shown in Figure 1. Dissolved oxygen stocks were 
calculated by vertical integration of the DO-depth prof Hep. 
at each station followed by summing over the entire harbour, 
A plot of i.he 'lirsolved oxygen stock thus obtained against 
sampling date is reproduced in Figure 2. Error estimates 
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and a curve of saturation DO stock have also been Included 
In Figure 2. This curve was analyzed in terms of the oxygen 
depletion processes. Computed from observed and saturation 

DO' stocks with appropriate tem,perature considerations, the 
oxyKen demand was ' approximately O'.Y x I'O* kg In Aprils and 
l..'"5 *x lO'* kg in August. If the oxygen demand In April Is 
assumed to be a result of waste discharges and bottom sedi- 
ments, and that the temperature coefficient of these processes 
is assumed to be (1.0^6)" , where at 
between to be bottom, temperatures In 



estlm 

the t 



ated that about 1.0' x 10' ' kg of 
otal demand, are used in August 



bottom sediments. This estimate may 
estim 



Is the difference 

August and April, It Is 
stock, or two=thirds of 
by waste discharges and 

nO'W be compared with an 



ate based upon oxygen uptake data already discussed. 

Under conditions of 7 to 9 mg/l DO at IS^C, Lucas and Thomas 
(1972) determined oxygen uptake rates of 0.08 to O.IM g/m - 
hr for the central basin hypoliinnl-on of Lake Erie. For 
similar conditions in Swedish lakes and rlver:^, containing o 
to \2% organic matter, PJdberp: and [lofcten ; 1V^7 3 



be 



oxygen uptake to 

rorit;iining ^1 to 

(I965) obtained an uptak 



13% 



0.0^1 to 
organic 



1 g/m"-hr 



u 



caT'l-on 



rate o: 



Thl-S rate was proportional to the 

concentration .in 

Harbour contained S to 1*3'^ 



With 
, Kd wards and 
0.0 3"^ to 0. U' 
. ^S power^ of 



determi ned 
river muds 
Ro 1 1 e y 
g/m^ ~hr . 
the DO 
the overlylnr water. As most oT Hamilton 

H-ilc matter In 1}64-S (OWRC, I966), 
ther.e figures wore selected a,-- ot'^-/. nit lug conditions 
-.imllnr to IJamilton Harbour. iJndiM' this ru^sumptlons , the 
following ei3tlmate may be niado of oxygen uptake for the 
iH^riod April 28 - August Ik, 197? (109 days). 
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izlng. that tlie oxygen uptake rate is a Function of 
omperature and dissolved oxygen concentration, the 
period was divided into tliree intervals: April 28-Iviay 
ne 1-30 and July 1-August l4. For the first interval, 
be assumed that the average bottom temperature and DO 
ilton Harbour are 8-10 C and 7-9 mg/l, respectively. 
Lng to Lucas and Thomas (1972), an appivximate uptake 
nde'r these conditions is 0.07 g/m^-hr, 
otal uptake for the harbour of 1-3 x 10 "^ kg, m 



wh.lch results 

6 



For June, if one assumes that the bottom DO has dropped to 2 
or 3 mg./l and the temperature is 12 C, the oxygen consumptions 
may be"estim,ated as 0.05 g/m^-hr, using Edwards and Rolley's 
(1965) function of overlying DO concentration. This gives 
an uptake of 0.8 x 10'^' kg for the harbour in 30 days. 

temperature 



average DO and 
and 14°C, respectively. (The 



In July and August, the bottom 

were estimated to^ be 1 mg/l 

measured bottom DO average for August 6-7 was I.3 mg/l for 

th,e entire bay (18 stations) and 0.5 mg/l for a subset of 

ttie 6 deepest stations). From ^ the data of Edwards and 
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Rolley, oxygen uptake is 0.04 g/m^-hr, or 0.9 x 10^ kr. for 
4i( days. 

This gives a total estimate of oxygen uptake of 3 x 10^ kg 
for the 109 day period^ which is certainly only an order of 
magnitude estimation due to the assumptions involved in the 
estimation, as well as the wide range of observed figures in 
the studies cited. This estimate can now be compared with 
the measured figures of DO stock and deficit, exchange 
through the canal and COD loadings to the harbour. 

As indicated in Figure 2, the dissolved oxygen stock in 
Hamilton Harbour was 3 x 10^ kg in April 1972. Were it not 
for photosynthesis, reaeratlon, and mass exchange through 
the Burlington Canal, the entire stock would be depleted by 
August. That the observed deficit (compared to saturation) 
in August was 1,5 x 10^ kg, of which 0.5 x 10^ kg was assumed 
due to phytoplankton decay Indicates the Importance of 
oxygen restoring mechanisms in maintaining harbour dissolved 
oxygen stock. 

A major source of oxygen to the harbour is mass exchange 
through the canal. Recording water chemistry and current 
meter data obtained in the canal were examined in an attempt 
to estimate values for mass exchange in the canal (See 
section D, mass exchange). Two estimates of dissolved 
oxygen mass exchange in the canal are as follows: 

(a) Mass exchange in the absence of thermal stratifica- 
tion, but with periodic flows present: A 10 day period of 
current meter data recorded at 10 minute intervals indicated 
volume flows of 3.M x 10^ m^ of water toward the harbour and 
5.7 X 10"^ m^ of water toward the lake. If harbour water 
contains an average of 6 mg/1 DO and lake water 11 mg/1 DO 
this represents a flow in 109 days of ^.1 x 10® kg DO to the 
harbour and 3.7 x 10^ kg DO to the lake. The net flow fs 

0.4 X 10® kg DO to the harbour. 

(b) Mass exchange with thermal stratification present: 
A similar analysis of 39 days of recording current meter 
data indicated a volume flow of 20.6 x 10 m^ of water to 
the harbour and 11.8 x 10^ m^ of water to the lake. Using 
the above estimated DO concentrations, estimates of DO mass 
flows of 6.3 X 10^ kg to the harbour and 2.0 x 10^ kg to the 
lake are obtained. The net flow is thus 4.3 x 10^ kg to the 
harbour. 

These estimates represent volume flows of water of approxi- 
mately 1 percent of the harbour volume per day in opposite 
direction. Thus the actual flow of dissolved oxygen to the 
harbour would be between these extremes, but likely towards 
the lower estimate, as the net direction of water movement 
Is towards the lake. Moreover, as dissolved oxygen concontrat 1 (^n; 
measured in the canal varied between extremes of less than 
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1 tO' 11 mg/1 for harbour water, and 9 to^ 1^1 mg/1 for lake 
water (with the lowest figures occurring In mid to late _ 
summer), actual masr, flow^ will vary widely with time. In 
addition, bathythermograph data indicated the presence of a 
three-layer system at times within the canal. Although the 
instruments used thus provide only a rough estimate of DO 
mass flow, the:y do indicate that the DO exchange with the 
lake can be large and is a factor in determining the dissolved 
oxygen stock. 

Another comparison of interest can be made with the Industrial 
and municipal COD and BOD loadings to the harbour. As 
enumerated In 197 2, daily loadings are 1.0 x 10^ kg/day. By 
comparison, the estimated daily sediment oxygen uptake for 
May 1972 is 0.37 x :i ^ kg (calculated from estimates used 
previously). Thus the dall>' discharges as about three times 
the estimated sediment oxygen uptake as estimated for a time 
when the bottom DO is relatively high. Several reasons for 
this discrepancy may be enumerated. Firstly, only a portion 
of the COD will be particulate and be deposited in the 
sediments The remainder will be assimilated by epllimnetic 
oxygen from ohotosynthesis or reaeratlon, or emitted to the 
lake through the canal. As most of the Industrial discharges 
are in the epillmnion, th.ls dissolved COD will not contriDur.e 
to the critical hypolimnetic depletion observed during ^ 
stratification. When the pronosed destratif I cation project 
is undertaken, this dissolved COD wJ '11 be distributed through- 
out the harbour and may conrributc to hypolimnetic oxygen 
deficits. However, as the low DO wat-r will be continuously 
circulated to the surface, n.ucli of this demand should be 
satisfied by atmospheric reaeratlon. In addition, if the 
particle size of suspended matter is large, much of the bUD 
will be deposited in the vicinity of the outfalls. J his may 
result In the accumulation of sediment to depths greater 
than that which will be readily oxidised. It should be 
recalled that sediment oxygen uptake is independent of 
depths greater than about 2 cm in laboratory tests (Stein 
and Denison, 1967; Edwards and Rolley, 19^^"?)- 

Although these estimates present a rough picture of the 
sediment oxygen uptake situaLion In Hamilton rlarbour, there 
are too many variables in the chemically and bLologically 
active systems studied to apply data obtaine.j elsewhere 
accurately to the situation existing in Hamilton Harbour. 
It is therefore proposed that sediment oxygen upt:^ke expcrl = 
ments be oerformed during 1975 in Hamilton Harbour, Although 
most attention should be paid to locations near the major 
discharges, all portions of the Harbour should be considered 
as the numerical model has IndlcHted the possibility of 
particulates being carried in currents around the northeastern 
corner to the north shore of the bay. Purthemore statistical 
analysis of sediment chemistry (section b) has indicated 
correlations between stations 20 and 26^, and between 269 
and 257. In order to interpret the oxygen uptaKe data, 
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sediment samples' should be analyzed for o^rgarilc carbO'nj COD,, 

B'OD and KJeldahl nitrogen,, and the bottom fauna prer,ent 

should be enumerated. 

To Inves^tlgate the posslbllltj of sediments at f';reater 

depths exe,rtlng O'xygen demand,, core samples shO'Uld be taken 
and oxygen uptake o^f Individual sections studied in the 
laboratory. This is pa„rtlcularly ' Important near the location 
of the proposed, dlffuser. 

It must be recognized, that in situ oxygen uptake experiments 
present considerable experimental difficulty and variability, 
as. the examples from the literature illustrate. However, 
actual data should be m,ore reliable than estimates based, on 
O'ther water bodies,. If these data can be compared with 
measured harbour DO' stocks a,nd correlated with analyses of 
industrial wastes, it Is hoped that an improved picture of 
dissolved O'Xygen status in Hamilton Harbour, and the effect 
of reductions of Industrial leadings, cO'Uld be obtained. 
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TABLE 1 : TYPICAL OXYGEN UPTAI^ RESULTS 

Sample Condltioins 



Ri'vers-8-11% org,an,ic 
matter 

Lake - 1.3% organic 
Lake - 12,2% organic 
Lake - 30-40% organic 
(eutrophic) 
Baltic Sea-PulpiTiill 
effluent site - 11.4 - 
17.7% organic 

Cellulosic fiber 

sludge. Strait of 

Juan de Fuco. 
Control area, Strait 
of Juan de Fuca 



Paper-waste sludge 



(a) :River muds 

(b) Mver muds affected 
by sewage 

Mver muds 

Org. C=4 to 13% 



in situ 0-10 C 



in situ 14 _C 
in situ 17°C 

in situ 5-7°C 



in situ 13-15 C 



in situ 11 "C 



in situ 11 C 



Oxygen Uptake (g/m -hr) 
Cat temperature (reduced to 
of measurenient) 15°C) 



0.013-01.06 0.037-0.15 



0. 


,02 




0. 


.01 




0. 


075-0 


.1 


0. 


04-0. 


13 



0.15-0.18 



0, 15 



0.046 



0.20 



0.060 



Author- 



Edberg S Hofsten^. 
1973. 

Edberg & Hofsten, 
1973 



Stein i Den i son, 
1967. 

Stein S Denlson, 
1967 





Lab. 


25 C 




0.13 0.06 


Hanes & Irvine, 






20°C 
15°C 




, 07 . 05 


1969 . 








0,07 






Lab. 


15°C 


Ca) 


0.05 winter, 
0.06 summer 


molley & 'Owens r 
1967, 




DO = 


7 mg/1 


(b) 


0.10 winter 
0.10-0.14 suimner 




(a) 


Lab. 

n 


20°C 


fa) 


0.035-0,14 0.05-0.20 


Edwards & Rolley, 
1965 


(b) 


20°, 


high speed 


Cb) 


1.2 0.85 






stirring , SS 










'^ 4500 ppm 









Cominents 



% organic matter by 

dry weight 



Laib. values 
lower . 



80-95% volatile matter 
in top 1-2 cm. 

Sludge depth >2.5 cm 
has no effect on upt^ake 
rate (det'd by lab test) 
Lab values higher,, 

"^58% volatile solids. 



Ca) Org. C '^ 6% dry 

LOI =11% 
Cb) Org. C - 7-12% 

LOI ■= 15-21% 

Seasonal effect - 
max:. O2 uptake in spring, 
S'uiirner, min. winter. 
40% O2 cons'umption due 
to macroinvertebrates . 
O2 consumption approxi- 
mately proportional to 
(DO) 0.45. 



TABLE 1: Continued 



Sample 



Conditions 



Oxygen Uptake (g/m'-hr) 
Cat temperature (reduced to 
of measurement) 15 C) 



Author 



Coraients 



River muds (a) organi- 
cally polluted (b) less 
polluted . 

Lakes Montague, South- 
port,, PEI 
A nranber of o^'ther 
Maritime lakes. 



Lab. 



Lab (Jenkin) 



L-ab Wank in, 
dumny J, bottle) 
all ■'^' 11°C 



(a) 0'.3~0,9 
Cb) 0.01 



0.03-0.04 
0,00-6- 0.02 



0,. 04-0. 05 
0.008-0.026 



Knowl e s ,,• Edwar d s 
S Briggs, 1962. 



Hayes & MacAuley^ 
1959 



Temperature not stated. 



Dunany Jenkins results 

lower. 
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FIGURE I ■' REGIONS USED IN CALCULATING DISSOLVED OXYGEN STOCK, 8 STATION GRID 
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HAMILTON HARBOUR STUDY 

SECTION F 

THE PHYTOPLANKTON 



SUMMARY 
The PHYTOPLANKTON OF HAMILTON HaRBOUR WAS SURVEYED ON 

six occasions during 1974. six sites were sampled 
through the euphotic zone. phytoplankton were 
identified to genera/ and enumerated to determine 
population composition/ abundance and distribution. 
Algal populations were low despite high nutrient 

LOADINGS. CyANOPHYCEAE WERE UNCOMMONLY SCARCE/ AND 
ChLOROPHYCEAE dominated all SAMPLES. GeNERA SUCH AS 
EUGLENA AND PhACUS WHICH ARE COMMON TO POLLUTED 
WATERS/ AND OTHER GENERA SUCH AS CrYPTQMONAS AND 
RhODOMONAS / INDICATIVE OF ORGANIC POLLUTION/ WERE 
ALWAYS PRESENT. DiATOMS WERE NOT ABUNDANT/ DESPITE 
THE FREQUENT OCCURRENCE OF StEPHANQDISCUS . REASONS 
FOR LOW PRODUCTION ARE DISCUSSED/ AND POSSIBLE 
EFFECTS OF ARTIFICIAL MIXING ARE PRESENTED. PHYTO- 
PLANKTON PRODUCTION CONTRIBUTED LITTLE TO THE TOTAL 
ORGANIC CONTENT OF THE HARBOUR WATERS. 



HAMILTON HARB€'UR 
THE PHYTOPLANKTO'N 



Despite Intensive chemical s.urveys O'f Ham,llton Harbour since 
1972s there has been no complementary study of the phyto- 
plankton. During, 197'^ six stations, numbers 20:, 262/2565, 
265 1, 25^' 3 and 253 were monltO'red for phytoplanktO'n co'mposi- 
tlO'n and abundance. Samples were collected by P'.asslng a 32 
oz bottle through the euphotic zone. Euphotle depth was 
considered to be twice the depth of the Secohl disc ohser- 
vatlon. 

Phytoplankton were identified to genera and standing crop 
determined in Areal Standard Units (one A.lS. U. = 40'0 square 
microns phytoplankton surface area). This is not an accept- 
able measure of blomasSj but does serve as a basis to deter- 
mine relative changes in phytoplankton populations. Future 
surveys should repres^ent blomass as algal volume. 

Only two stations (253 and 20) were monitored continuously 
from June 25, 1974 to October 2^ 1974. Because of the sampl- 
ing routines and the method of enumeration^ spatial distribu- 
tions of phytoplankton in 1/ - harbour could not be determined 

Only relative temporal vari ati^^ns can be commented on, with 
general comparisons between stations. Table 1 outlines the 
phytoplankton standing crops in Hamilton Harbour. During 
197^ station 253 tended to have the higher standing crops 
reaching a maximum on Aur,ust 21^ of 3052 A.S.U. (% 20 Uf^^",' 
Chlorophyll - a). This is extremely low ccjimpared to chloro- 
phyll - a levels exceeding 3,0 x 10^ ug-l~" recorded in 
Cootes Paradise by Harris (197^)- It is not known if the 
standing crop at station 253 was the result of higher 
product ion 5 or was an artifact of the mosaic of transport 
mechanlsm,s operating within the harbour. 

Station 20 had low phytoplankton populations (see Table 1, 
and Figure 1), but did reach a maximum of 35,200 A.S.,U. in" 
June, 1974. Similar to the other stations, Chlorophyceaen 
species dominated the algal component of the seston at this 
,slte. 

The highly eutrophlc nature of the harbour is reflected in 
the species composition of the phytoplankton. Genera such 
as Euglena , Phacus , Crypt omonas ^ and Rhodomonas reflect the 

high level of o^rganle pollution in the harbour. The large 
centric diatom Stephanodiscus is common in eutrop>iic bodies 
of water and was present at all times of the year in the 

harbour. 
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There was a conspicuous lack of blue-p;reen ali"^ae (Cyano- 
phyceae) in all samples. As Hamilton Harboui' has a rela- 
tively low pH (7-0 - 8.0) and relatively moderate alkalinity 
(80 - 110 mp;. 1 ^ as CaCOg) there would be abundant CO2 to 
support both green and blue-green a.lp:ae. The low production 
of blue-green populations might be a response to the low 
thermal stability of the water column in the harbour, ri.s th^r- 
therrnocline usually stays at about four metres depth d-^^spitc^ 
increasing, radiant energy. This lack of stability is quite 
likely a result of large horizontal flov/s throughout the 
basirii and disturbance of the water column by shipping. 
Many Cyanophycoae are adapted to stabilized water columns, 
and the lack of such might limit their production. Gem^r-n 
of the Chlorophyceae such as Scenedesmus , Coe la strum and 
Oocyst is are not adversely affected by low i'.tabi In ty j and 
demonstrate a growth not en tin] in tlie liart'O r' (r^ee I'M ^;ui'es 
1,2 and 3) • 

The phytoplankton populations outlined in Table 1 are low 
consideT'ing the hypertrophic nature of the harbour. Nutr"!eut 
levels in Hamilton Harbour (Section A; Table 1) suggest h/ig^h 
nitrate availability, but surprisingly low o^tiiophosphate 
levels. These phosphate results are spui'^ious. , as soluble 
phosphate levels are near the minimum level for' detection, 
and increase directly with chlorophyll - a. Lf soluble 
phosphate was lim.iting an inverse I'elationsh ip would be 
expeci^ed. Sam.ple collections were not frequent enougii to 
comment en the phosphate - production dynamics oT the harbour- 
Compared with total organic carbon_levels of pr-evious years 
(M.O.E., 197^0, of 6.0 - 8.0 mg. 1 ^. the phytoplankton com- ■ 
ponent is relatively small (1 mg . 1~M- It would appear 
that much of the organic content of the harbour is from 
waste discharge into the epiJimnion, Possible effects of the 
phytoplankton on the total d1. ^solved oxygen budget (Section 
}-:) ai'o not known, as measuremt.-nts of primar'y productivity 
(light and dark bottle method) have not been made. It mur.t 
also be established if there were any localized g.rowths 
within the hai^bour, and a realistic determination of total 
phytoplankton production made. 

Zooplankton data is scarce 3, but net collections on June 25, 

197'^^ suggested moderate zooplankton levels.' The hauls were 
just through the euphotlc zone,, and are not indicative of 
the water column. It would appear that zooplankton popu- 
lations were not sufficiently high to maintain such low 
algal standing crops. (U'hlmann, 1971). There is a cons- 
picuous absence of Da.phnla in all colleetlons. Further 
study is required to determiine reasons for this .•i,bso'ncc% l>ii'i 
ctiinhhwi wMli I nw rilgn/t pi-ndnic 1, 1 on, , t,iixlctl..y iiiii;".t nn't. h.f 
vjv. f I'iuUcd , A -•■ccu'iid |H>s,i,;I b le lliiiltliig. s,ourco Is. Lh,e 
LurbldJiv oi' the harbour watei's. Secciil disc o^bservatlons 
iiKiicati' low transparency which must be a result of high 
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The vertical light attenuation coefficient of Hamilton 
Harbour Is approximately 1.5, but submersible photocells, are 
required for more accurate determinations. Such a high 
attenuation, coefficient would suggest there Is very limited ■ 
phO'tosynthetic available radiation penetrating the water 
column. This will severely limit phytoplankton production 
(see Section G on artificial mixing) in moderately deep 
water columns. 

Ar'tlf'tclaL riilxlng, of llamllton illatH)OUF could rcrailt in a 
cliani';e in tJie abundance and composition of the phytoplankton, 
and will in turn affect the zooplankton. Mixing can result 
in increased dissO'lved carbon dioxide in the euphotlc zone, 
and comhined with low stability of a mixed water column 
often results in a shift of bluegreen algae to green algae,, 
Hamilton Harbour, however, is presently dominated by members 
of the Chlorophyceae . If there is a change in the algal 
composition it might well favour diatoms, and genera such as 
Stephanodiscus , Pragilaria , Aster ionella , Melosira , and 
Cyclotella would Increase in abundance as they are adapted 
to vertical mixing in deep water (see Lund, 1971). 

It Is not possible to pr^eciict the effects nV mixing on the 
algae. Physical and ciienilcal ch3n;':er; eKeT't numen-nis .in- 
direct effecti- on th^^ pirH]kton, anii tin? re^n'^^nse ±r. tn the 
summation of the rf'fecti.. \\ i.^. not unlikely that tl:e 
abundance? oV r'l^een alg,ae iv i i . ilr-.' if the turbidity oT uJie 
harbour remnins hi^-h, and t h.o ::i}^:;il ct'lls are mixed tiirougb 
the watur column. The rM-LL;:. '^ i L 1 not obtain sufficient 
f)hotosynt"not 1 c enc\'r.y to 'j :>s 1 in 1 1 n to carbon under these 
c L r c; u m s tan c t^' r, . i ' i r: u r c s i , ..' , a nd ^ i 1 1 u s t r a t o a a c eve it s e i n 
piiy toplanktori aburulanct.' in thn puLuinn of l'?7^*. This Is 
assocl:ited with thr brciKini" u^.^wn of' thor'rnaL stratification 
and the j^esultant increar.f In tiie deptn of vertical mixing. 
Production In Hamilton Harbour* can be limited by the physical 
pi-^o c e s s e ;,' vi h i c i i i n f 1 uc n c e t h c ni i. x 1 ng vi e p t h . 

During 1975 an_ Intensive biological survey will be initiated., 
Phytoplankton composition, abundance and distribution will 
be determ,ined. Primary production experiments will be set 
up, and blo,assays used to assess the production potential of 
the harbour water. Relative abundance and composition of 
zO'Oplankton .and fish popula.tlons will be determined. This 
survey will Include an investigation in changes of bottom 
fauna and b,acteri,al populatio^ns ,as there is suitable back= 
ground data to determine the effects of mixing on these 
parameters. 
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TABLE 1 
Phytop'lanktO'n Abundance 197'4 
All Figures are in Areal Standard Units ml'^ 
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Date 20 262 256 265 25^ 253 
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FIG. I SEASONAL VARIATION OF PHYTOPLANKTON IN HAMILTON HARBOUR 1974, 

(COMPOSITES OF EUPHOTIC ZONE, STATION 20). 
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FIG. 2 SEASONAL VARIATION OF PHYTOPLANKTON IN HAMILTON HARBOUR 1974, 
(COMPOSITES OF EUPHOTIC ZONE, STATION 254). 
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FIG. 3 SEASONAL VARIATION OF PHYTOPLANKTON IN HAMILTON HARBOUR 1974, 
(COMPOSITES OF EUPHOTIC ZONE, STATION 253). 
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Artificial M 1 x L n i ^, a n d Wa t e r 
t^juailty ill llaiiilltoii llai''bom' 



Introduction 

Artificial mixing has been used to improve water quality and 
manage biological production in aquatic systems. Anoxic 
conditions J, such as the hypollmnetic depletion of oxygen in 
Hamilton Harbour, can be ameliorated by the correct mixing 
process. VariO'Us methods of mixing are common (Past^ I968). 
The mO'St frequently used and mO'St economical form of mixing 
Is aeration. This technique Is often applied to GituatiO'ns 
of low oxygen levBls (Imhoff, 1969; Patrlarche, 19t3l) to 
restore dissolved oxygen to the system,. Mechanical pumping 
of water Is a more expensive alternative to aeration, but 
does Increase the effective exchange area of the wateP'-alr 
Interface. 

Hooper, Ball and Tanner (1952) were the first to report the ■ 
use of artificial mixing In North Am,erlGa. Complete destra- 
tlficatlon was not achieved by the mechanical pumping of 

hypollmnetic water Into the eplllmnlon, but phytoplankton 
prO'duc t ion was s t Imulat ed . 

Reasons for Implementing art i P 1 ,- II. I nilxing vary according to 
water use (KPA, 1973)- ImpounJiru-nt s foi' fros^nwater supply 
use artificial mixing to prevent hypolimnetic anoxia, thereby 
Increasing the volume of potable water (Ridley, Cooley and 
Steel, 1966; Ridley, 1964; Rl(ii.5y, 1970). Natural lakes are 
mixed to prevent winterkill of fioh (Burdick^ 1959; Halsey, 
1968; Schmitz and Hasler, I'J'jS), but not always successfully 
(Patrlarche, 196I). YouriL (IQpb) I'alled to reduce the 
levels of nutrients in hypertrof)!' 1 c lakoG In l^ltDrida by 
aeration, possibly as a result o? a poorly established 
oxidized microzone at the mud-wai:or Interface. Improvement 
of polluted lakes by aeration b.as been reported (heath, 
1961), but in most case.s niixiji^'' 1 1- used as a management 
process, not as a corrective ni' iako renewal mc^chanism. 
There is little known work on the-- ofTects cjf mixing on 
bodies of water which have larr:e i n^iur^trlal discharges and 
intakes. High COD's are ucuall;, treatea before discharge In 
many circumstances, and artificial mixing would not be as 
effective as pretreatment processes. 

Artificial mixing ir> oCton a,~,ed \o [;]-event aj' destroy tiiermal 
stratification. Independent m. Ixlnt'-; of i-ae opllmnlon (Bernhardt, 
1967) or the hypolimnion (Hpeece, l-'70) is T)osslble, but is 
used in Isolated circumstances as a result of inefficient 
mixing techniques or [ilgn cos^ts (see Dunst, 197^)- Hypollm- 
netic aeration has a projiii sin,-:^ fucui-e should appropi'iate 
methods be found. 
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A major review of the prlnclp 
clal mixing has been made by 
compiled an extensive list of 
efficiency of vertical artlfl 
reservoirs. The American Wat 
reviewed the results of art if 
and noted that 90% of mixing 
A review of mixing comparing 
with types of mixing and prob 
the end of this section. 



les and techniques of artlfl- 
Fast (1968). Symons (1969) 

references on the use and 
clal circulation In lakes and 
er Works Association (1970) 
icially induced circulation^, 
trials were considered successful 
depth, area and volume of lakes 
lems encountered is appended at 



Re_sult3 of Artificial Destra_tif ication 

General improvements in water quality associated with artifi- 
cial mixing are: 

1. Elimination of taste^ odour and colour 
(A.W.W.A. 1970') 

2. Improved transparency. (A.M.W.A. 1970) 

3,: Larger habitat for aquatic life, (Inland 
Fisheries Branch, 1970). 

H, Increased dissolved oxygen, (Posey and 
De Witt, 1970; Imhoff, I969). 

5. Reduced algal blooms but not necessarily 
implying reduced algal production. 
(B^ernhardt, 1967; SymO'ns, Carswellj and 
Robeck, 1969; Bakus, 1973). 

6. More neutral pH,, as CO2 Is released from 
hypolimnlon. (Riddick, 1957; see, however. 
Lackey, 1972). 

7. Decrease in dissolved metal content, 
especially Fe and Mn. (Symons, Carswell and 
Robeck, 1969; see, however, Lackey, 1972). 

Chemical Results 

Chemical effects of artificial circulation are dependent on 
the efficiency of the mixing process (Pasveer, 1965; Symons 
and Robeck, 1969). Inverse clinograde distributions of Pe , 
Mn, N-NHg, CO2, alkalinity and conductivity associated with 
a discontinuous temperature profile can be corrected. 

Iron and raanganeBe precipitate out as hydroxides in aerobic 
conditions assuming the pH Is sufficiently high. The main- 
tenance of an oxidized m,lcrozone will prevent Tecycling of 
these metals and their salts, iron phosphate complexes In 
particular (O^ntarlo Ministry of the Environment, 197'^'). 
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A possible deleterlo'us point of ai»tiflclal mixing is a 
general increase In to^tal phoBphates as a result "o^f disturb- 
ance O'f the lake floor. As.so^el ated, with this problem are 
decreased sinking rates of orT.anlc and Inorganic particles, 
and extended periods O'f high turbidity. 

The development of an oxidized mlcrozone at the sediment- 
water ^ Interface is an essential feature of successful 
artificial mixing. This microzone might satisTy some 
surface sediment oxygen demand, but does not necessarily 
affect _ the oxygen demand O'f sedlm,ents as determined by 
laboratory analyses, and by the requirement of at least 
Intermittent mixing to maintain aerobic conditions. This 
microzone retards the recycling of nutrients (Mortimer^ 
19'41; 19^2; 1971), and the sediment will act as a nutrient 
sink (see, however, Tessenow, 196^1, 1966), with a gradual 
build up of reduced material below the microzone. "There is 
little evidence that artificial mixing ever satiates the 
sedlm,ent oxygen demand. 

Biological Results 

Changes in bacterial populablonu as a result of artiricial 
circulation are difficult to determine, and little work has 
been done in this field. Anaerobic decomposition of organic 
matter releases NH3-N as an end product. A second source of 
ammonia is the reduction 01' nitrate by bacteria, and hence 
ammonia tends to accumulate jn th. I:ypolimnion. ' Anaerobic 
decomposition Is slow, and c^'ten I'e.-.ulc^; in an accumulation 
of organic material In and zibove the sediment. Artificial 
mixing should favor the pre-.^exicc of nitrifying and heterotro- 
phic bacteria, P0^3ultl^g Jn .Incr^oased decomposition of 
organic matter, particulai- M y Tt the sediment-water inter- 
face, and a decline in the level of ammonia. There is a 
possibility of the accumulation of nitrates as a consequence 
of continued artificial mlxlnr;. 

Dissolved cat'bon d.i oxide and alkalinity are to a large 
extent interdependent on the pH and plant communities of the 
lake. ^C02 tends to accumulate in th(-^ hypollmnion of a 
stabilized water column. Art, UMciil mixing can produce an 
orthograde ratio of C02/?-lkaiJ nl by Interdependent on pH. 
Changes^in the available carbor. associated with mixing can 
result in a significant response in phytoplankton t^rowth and 
succession. Composition of the standing crop might be 
altered, as King (1970) noted that high levels of dissolved 
carbon dioxide favored ■^r'ecn flagellates. Rlue-green algae 
can utilize low CO2 levels (common Jn an epilimnion) hence 
low carbon dioxide regimes t^nc] to favoi- their growth in the 
euphotlc zone. Mixing tends to inoj>ease the available 
carbon dioxide in the euphotlc! zone, and green algae are 
favoured. This, howevej^ is 3 broad i'-^eneralizat 1 on, and in 
reality many other factors such as temperature, vertical 
attenuation and modification of raai^int enerry , and mUlng 
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depth can have an equal or greater effect on the phyto- 
plankton response, Hamilton Harbour is already dominated by 
green algae, probably as a result of the cumulative effect 
of all these factors, 

Tiie Gompooitlon J xjer Iodic ity , duration and 3lze of phyto- 
rlankton Gtanding crop can be altered by the mixinfr, proceso. 
Phytoplankton populations are often Interactive and inter- 
dependent on t-he physical-chemical pai'ameter:- wb/icli are 
nf^\-ct'':d by artificial mlxinr. These indirect effects of 
mixinp:, coupled with trie direct effects of mixing; on cells 
(vertically homogeneous alp.al populations), produce shiort 
and lon.p; term effects on the populations, Becau::<e of the 
complexity of the interactions, reservations must be made 
when predicting the effects of m.ixing on phytoplankton and 
zooplankt:on production and distribution (Tailing, 1970; 
r'ast, 19/'l). These in turn can combine with the physical- 
chemical changes and affect bottom fauna and fish population: 
(Murphy, 19 
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While preventing the growth of species adapted to stabilized 
water columns, such as many species of Cyanophyceae , artifi- 
cial destratlf icatlon can stimulate the growth of other 
algae. Green algae (Chlorophyceae ) can (3on;inate when eutro- 
phlc bodies O'f water are mixed, especially if tbe v;ater 
cO'lumn is optically shallo'w. A water column is optically 
shallow when the mixing depth (Zm, ) and the cuphotic depth 
( Z eu ) are s imi lar , 

Diatoms might dominate if the mixed water column is optically 
deep (Zm, >> Z eu)^ and other physical-chemical features are 
favourable . 

Lund (195^Ui 1955 9 19T1) has noted that mixing can alter the 
seasonality and size of populations of MelO'Slra - Islandlca 
(Ehrenb) Kutz. subsp, subarctica 0. Mull. In the English 
Lake District. Ridley (1970) reported delays in the vernal 
maxima of Stephanodiscus astraea (Ehrenb) Grun. which were 
possibly related to prolonged turbidity caused 'by artlfleial 
mixing* 

Haynes (1973) fO'Und mixing resulted in a decrease of 
Aphanlzomenon f los aquae Ralfs,, and Steel (1972) noted that 
vernal growths of phytoplankton In Thames Valley Reservoirs 
we^re inversely related tO' the mixing deptti and vertical 
light attenuation. 

Phytoplankton Productivity 

PhytO'plankton production per unit area in m^' can be formulated 

as : 
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mln 



11, = algal concentration (mg chlorO'phyll a 

_ 3 
m ). 

Pmax - maximum photosynthetlc rate 

vertical attenuation, coefficient of 

_i 

photosynthetlc available radiation (m ). 

Is a correction fo^r spectral modif lcatio.n 

with depth. 

lO' - incident radiation (K er,gs cm sec ) 

Ik = light intensity at onset of photos^ynthetlc 
light saturatlu . (K ergs cm sec ) 



Zm - depth of m,ixed water colum,n (m) 



Such a model (Steel, 1972) of the energy flows associated 
with phytoplankton populations, is advantageous when eutrophic 
bodies of water are mixed. When a e^j^„ Zm is great. 



^ -o £ . Zm 
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Ik 



will tend to zero, and the photosynthetlc profile will be 
contained within the mixed water column. Such, a body of 
water is optically deep ( Zm >> Z eu). Scotch Block Re-servolr 
and Buchanan Lake are optically deep (see Appendix 1, No. 14 
and 16) bodies of water destratified by the Ministry of the 
Envlronm,ent . Valens Reservoir , also artificially mixed, is 
optically shallow with light penetration near or to the 
bottom. 



G=5 



As Hamilton narbou.r is ciar>LO.f led aL> eutroplilc-hypertrophic , 
and. is optically deep, it 1s expected thiat priniai'y production 
will be limited by avaiJabie scjlar radiation. Despite the 
fact that the volume of wnter wh:ch will support aquatic 
life will probably be increaiied hy nrtifj'jial niixinr;, trie 
actual b^iolob';ical pT'oducticui oJ' the hai'bour* will not iiecOL^- 
3 a r i 1 y i n c r e .• 1 1: e . 

E X p e c t e d [ ''. e '3 1^ 1 1 r^ Hami lton iti rhnur 

The-" 'naln real o^ the 'je::;ti'atif ication of Marnilton Hnrbour It? 
\:-j i".^iLCvr tae anoxic condit1(Mis in the I'lyp'^.'^lirnni on . It ai> 
hoped thiat th:'L-; will resuJt in better waten' 'lual'iy, aithoufjh 
chemical una \ iolor.ical re^iponiujs are impopasiblc to predict 
at this tiiiTie. ThiC irregular ^-hape and depth of the harbour, 
the nature of the mud-water interface and th..^ p;eor:r jpiiica 1 
environment of the harbour are some of the many unknown 
factor^', which will influence t?\e responoe of the har'oour to 
artificial rriixinp;. Horizontal transport proco:-;r.es In th<^ 
harbour are relatively larf;e (JVl.O.K., l':)7'\) rind rdiould 
assist the niix.inr'" process, 'out would not necessarily ensure 
aeroui.c conditions at tlit^ iniciHn^one withiout the <'K's is,Lance 
of an aerator. Locating-: the ruei^atoi^ at tlie doepe.'-.t point 
will heir maintain tlie oxidizcai riLLci'ozone tlirourjiout the 
harbour. An a,";ency m.lxint;^ efi'ect of t\](- wi;ui jiilrJit competr?ate 
foT' the j rreg;ularities in the shape of the harbou.r. Cons, Ider- 
in,--^ the mixins" depth of Hamilton Harbour (> 20ni), and that 
the water is op'cica.lly deep, \t is expected that the phyto-* 
plankton composition mip-ht shift towards species adaptive! to 
low light and mixinj^ conditions. Many diatoms are adapted , 
to such conditions., anci It is, pv3S^s>ible there will be an 
Inc r' ea s e in d i rt t o m v. n p u 1 a t i o ns i n t he ha r b o u t' waters. 
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APPENDIX I 

A Compilation of References and Results 
on Artificial Mixing and Reaeration 

Methods 



Turbines equipped with reaeration plates, 
installed on dams. No specifications 
given . 



Results 



Hydroturbine aeration is approx- 
imately 25% efficient. It can 
increase the dissolved oxygen 
of a moving stream 900 Kg. of 



oxygen day , 28 m^ sec" 



?■ 



Q 



(2) Posey, F.H. and DeWitt, 
J. W. 1970. Effects of re- 
servoir impoundment on water 
quality. 

Proc. Amer. Soc. Civ. Eng. 
96: 173-185. 



(A) Reaeration of di-^.charge by: 

1) multilevel penstock intakes 

2) submerged weirs 

3) supplemental high quality spills 

4) turbine aeration 

(B) Destratification by: 

1) air hydraulic guns 

2) perforated air tubing 

3) mechanical pumps 

4) diffused-air-pimiping 

Allatoona Reservoir 4.5 x 10%^ mixed by 
five cross-shaped air diff users. Compres- 
sors (5) supply 250 cfm at 100 psi 

3 - ^ 2 

(7m min at 7 Kg. cm" ) . 



(A) Turbine aeration successful. 



(B) Diffused air destratified 
Allatoona Reservoir. 



Reference 



Methods 



Results 



(3) Symons, J.M. Car swell, 
J.K., and Robeck, G.G. 1969 
Mixing of water supply re- 
servoirs for quality control. 
J. Amer. Water Works Ass. 
61: 322-344. 



Destratify two lakes by diffused air 
Boltz Lake*: 

19 m. (max. depth) 

3. 9x10 ^m^ (area) 

3.6xl0^m^ (vol.) 



Falmouth Lake 

13 m. (max. depth) 
9. 1x10 ^m^ (area) 
5. 7x10 ^m^ (vol). 



Reduction of levels of Fe, Mn, 
SO2 occurred, accompanied with 
a temporary decrease in phyto- 
plankton populations. 
It was also noted that quality 
control should begin before 
any quality deterioration. 



Q 
I 

M 



Both lakes mixed by compressers operating 
1 
3 — 
at 3 . 2m min 

*Mechanical pumping trial used for a 
limited period. 



(4) Speece, E.R. 1970. 
Hypolimnion aeration. 
J. Amer. Water Works Ass, 
72: 6-9. 



Artificially induced aeration of the 
hypolimnion using a 2m diameter pipe. 
Alternative methods are discussed. 
The hypolimnion of Lake Bret was aerated 
by forcing water up the pipe and dis- 
charging it through laterial parts which 
were below the thermocline. 



Advantages of hypolimnetic 
aeration are : 

(1) temperature budget 
is not increased. 

(2) still maintain oxidized 
microzone. 

(3) cold water fish 
propogation. 



The dissolved oxygen of the 
Hypolimnion can be increased 
from 5mq. l" to 10 mg. 1 



). 



f 1 
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Reference 

(5) Imhoff, K.R. 1969. Oxygen 
balance euid artificial re- 
aeration of Lake Baldeney and 
the lower Ruhr. 
Advances in Water Pollution 
Research. Proc. 4th Internat. 
Conf. Prague. Pergamon Press. 
761-774. 

(6) Symons, J.M. , Irwin, W.H., 
Robinson, E.L. and Robeck, 
G.G. 1967. 

Impoundment destratification 
for raw water quality control 
using either mechanical or 
diffused air pumping, j. 
Amer. Water Works Ass. 59 
(10): 1268-1291. 



, Methods 

Lake Baldeney 

10m (max. depth) 

2. 6x10 ^m^ (area) 

was reaerated using various methods in- 
cluding, a floating aerator, compressed 
air aeration system, floating aerator, 
and turbine aeration. 

Same lakes as in reference 3. Boltz Lake 
mixed at 19m depth. Sixteen diffusers 
arranged in a cross pactern as before. A 
1.6 X 10%^ day" mixed-flow pump was 
stationed at the deepest part of the lake 
permitting water drawn from the bottom be 
discharged at the surface. 



Results 



More feasible to apply artificial 
reaeration than to increase re- 
moval rates of sewage plants be- 
yond 85% or to apply tertiary 
treatment . 



Destratification is more effi- 
ciently achieved by aeration, as 
high discharge velocity head and 
hydraulic losses reduce efficiency 
of mechanical mixing. 
Blue green algal populations 
replaced by green, with a decline 
in total standing crop. 



Reference 



Methods 



Results 



(7) Patriarche, M.H. 1961. 
Air induced circulation of two 
shallow Michigan Lakes. 
J. Wildlife Management. 25 
(3): 282-289. 



Aeration with 150m perforated tubing 

1 

3 ~ 

attached to 0.14m min compressor. Two 



shallow lakes 2m and 1.5m depth, with 
areas of 2. 9x10 ^r 
ively were mixed. 



areas of 2. 9x10 ^m^ and 6.8xlo'*m^ respect- 



Depletion of dissolved oxygen 
causing winter kill of fish was 
not ameliorated. Large demand 
for oxygen by decay processes in 
the rich organic sediments was too 
great to be satiated by mixing. 



Artificial mixing resulted in 
large phytoplankton growths . 



Q 
I 



(8) Ridley, J.E. , Cooley, P. 
and Steel, J. A. 1966. Control 
of thermal stratification in 
Thames Valley Reservoirs. 
Proc. Soc. Wat. Treat. Exam. 
15: 225-244. 
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Large reservoirs* 18-25x10 m mixed by 
jet-type inflow, or hydraulic pump 

located on pontoon. Jet type inflow 

_i 

(4.5x10^ - 9. 0x10 ^m^ day" ) or hydraulic 
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pumping (1.8x10 m day 
adequate mixing. 



achieved 



Improved biochemical quality, but 
reservations expressed concerning 
phytoplankton production and 
distribution. 

Increased volume of water avail- 
able for domestic uses. 



,7_3 



*A reservoir of 3.5x10 m was also mixed 
by jet-type inflow in 1972, (author's 
note) . 



Jet inflow allows directional 

mixing within the basin. 
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Reference 

(9) Lackey, R.T. 1972. 
Responses of physical and 
chemical parameters to 
eliminating thermal stratifi- 
cation in a reservoir. Water 
Resources Bulletin 8 (3) : 
589-599. 

(10) Bernhardt, H. 1967, 
Aeration of Wahnbach Reservoir 
without changing the temper- 
ature profile. 
J. Amer. Water Works Ass. 
69: 943-964. 



Methods 

Parvin Lake Reservoir (max. depth 10m, 
area = 4. 0x10 ^m^) was mixed by two 
aerators supplied by 2.1m3min compres- 
sor at 4.2 Kg m" . 



Wahnbach reservoir is 19m deep, with a 
volume of 4. 16x10 ''m . Mixed with one to 
three aerators operating at 2, 4, and 6 
m^min with a 36.5 Kw electrical 
compressor. Aerator designed such that 
water flows through a 2m diameter pipe, 
with lateral discharge apertures. 



Results 

Change in heat budget of the Lake. 
Alkalinity, pH, conductivity, and 
total Fe, Mn, were not signifi- 
cantly affected. 

^ ++ . 

Ca increases slightly, and 

SOi^. decreases . 

Until complete destratification 
was achieved with two aerators, 
Mn increased in concentration 
after full circulation Mn de- 
creased, biological production 
remained unaffected. Also used 
mixing to lower thermocline but 
noticed that O2 demand was still 
high even after complete 
circulation. 



Reference 



Methods 



Results 



(11) Roeber, J. A., Noga j , R. J. , 
Ciabattari, E.J. and Hurwitz, 
E. 1965. Floating aerator. 
Water Works and Wastes Engine- 
ering. 2: 64-66. 



Stream aeration of a shallow river (1 
metre) . Aerator was designed to lift 
water above the general water level and 
spread the water across the surface. It 
was estimated that 8.2x10 Kg oxygen per 
day would be required using 40m of air 
per 0.454Kg of O2. This experiment 
would require 10,000 hp for compression, 
but aerator was only 100 hp. 



Efficiency of aeration depends 
on oxygen tensions existing in 
the water. Energy efficiencies 
varied from .3 Kg to 1.1 Kg 
of oxygen per hp hr. 



o 
I 



(12) Biederman W.J, and 
Fulton, E.E. 1971. Destrati- 
fication using air. J. Amer. 
Water Works Ass. 63: 462-467. 
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Waco Reservoir has a volume 1.28x10 m 
with a maximum depth of 26m. A gasoline- 
powered compressor 25 hp delivering 
_i 

3m min of air was used to destratify 
the reservoir. Replaced plastic 
delivery line with galvamized pipe as a 
result of breakage problems. Diffuser 
was 15m of galvanized pipe with 
^ 3mm holes. 



Solved taste aind odour problems 
associated with stratification 
No algal blooms were reported. 
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Reference 



Methods 



Results 



(13) Haynes, R.C. 1973. 
Some ecological effects of 
aurtificial circulation on a 
small eutrophic lake with 
particular emphasis on 
phytoplankton. Hydro- 
biologia, 43: 463-504. 



Kezar Lake is a small lake, 9m mcix. 

6 3 

depth, with a volume of 2.0x10 m . 
Four compressors of 3 hp each, collect- 



A noxious blue green alga was 
replaced by chlorophycean species 



Ca, Mg, K, Cl, Si02 were not 

_i 

ively discharged air at 3m^min . P.V.C. affected. Increase in Total 
delivery line terminated at a junction of POi+ associated with increased 
four arms each with 3 ceramic diff users. suspended material. 



I 



(14) Ontario Water Resources 
Commission, 1971. 

A.I.D. for lakes program. 
Artificially induced destra- 
tification. Toronto, Ontario, 



A 2 hp compressor d'liiverying 0.3m min 

of air was used to destratify a lake with 

maximum depth of 13.4m, area of 

8.6xl0'*m^, and volume of 4. 3x10 ^m^. 

Turnover time was 4.8 hrs, and vertical 

_i 
velocity due to air is 30 cm/sec" . 



Increased production of phyto-and 
zooplankton . 

Improved water quality in 
regard to levels of H2S, COar 
NH 3 , and Fe . Mcirginal improve- 
ment in bottom fauna amd 
sediment chemistry. 



Reference 



Methods 



Results 



(15) Ontario Ministry of the 
Environment. 1972. 
Artificially induced 
destratification of Valens 
-ake. Toronto. 33pp. 



A 2hp compressor delivering 0.28m min 
destratified a lake 4.0m deep, 
4. 9x10 ^m^ area, 1.2xl0^m^ volume using 
six diffusers of 6 metres length with 
20 holes in each. Velocity of water due 
to air "^ 8.5-12.0 cm sec . 



Iron and hardness of the water 
decreased. Dissolved oxygen 
improved in area of the diffuser, 
Odour problems (H2S) were 
improved, and an extremely high 
density of macroinvertebrates 
were noted. 



Q 
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(16) Ministry of the Environ- 

=ient, 1973. 

Artificial destratification 

of Scotch Block Reservoir. 

A progress report in the 
changes in water chemistry 
through the period May to 

November, 1973. Toronto. 



Scotch Block Reservoir (1.8x10 m with 
cm cirea of 5.3xl0^m^) is "^ 15m deep. 
The reservoir was mixed by a 5hp 
compressor delivering air at a rate of 



0.56m min 



Fifteen metres of neoprene 



pipe connected the aerator-delivery line 
system. Aerator sections {6m length with 
ten 1.5mm apertures) were alternated 
with delivery sections of 3(}m length. 
Six aerators were used. 



Reduced downstream impairment which 

had resulted from stratification 

in the reservoir. Evidence of an 

excessive O2 demand by the 

sediment, which resulted in no 
appreciable increase in dissolved 

oxygen content of the bottom 

waters. Controlled levels of 

H2S, NH3 and CO2. 



(17) Eunpu, F.F. 1973 
Control of reservoir 
eutrophication . 
J, Amer. Water Works Ass. 
65 (4): 268-274. 



Occoquan Reservoir extends over 
6.9xl0*in^ with a volume of 3. 7x10 V^. 

Three rotary compressors with a total 

capacity of 6.3m^min were attached to 

2. 3x10 ^m of perforated tubing, which was 

fed by 7. 6x10 ^m of supply tubing. 



Reduction of H2S and elimination 
of anaerobic conditions. 

Reservoir was isothermal with no 

algal blooms since artificial 

mixing initiated. 
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HAMILTON HARBOUR 



DESCRIPTION OF HARBOUR 



Hamilton Harbour is an enclosed body of fresh water at the 
western end of Lake Ontario approximately 8 km (5 miles) in 
the east-west direction with 4.8 km (3 miles) in the north- 
south direction. It is connected by a ships' canal, 732 m 
(2400 ft) long by 107 m (350 ft) wide and 9 m (30 ft) deep 
to Lake Ontario. The harbour contains approximately 2.8 x 
108 m-^ (9.9 x 109 ft^) of water with a mean depth of 13 m 
(43 ft) . Hamilton Harbour receives drainage from a water- 
shed of 500 km2 (190 mi2) through several small creeks which 
drain into the harbour with a total annual flow of 1.27 x 
108 m3 (2.8 x 104 MIG) consequently, the natural throughput 
is approximately 45 percent per year. 

The industries located on the highly developed south shore 
use 27 m3/sec (510 MIGD) water as an industrial supply, 
returning a similar amount of effluent to the harbour. They 
also depend upon the harbour for shipping of raw materials 
and finished products. 

Municipal utilities also use the harbour as a receiving 
water for 3.2 m^/sec (61 MIGD) of treated wastes and for 
untreated storm sewer overflows estimated at 0.1 m3/sec (1.9 
MIGD) during overflow periods. Together with the natural 
flow, this results in a total annual flow of 2.3 x 108 m^, 
or a throughput of approximately 80 percent per year. 
Presently, the harbour is also used for recreational boating. 

As the harbour shoreline is extensively developed; the 
industries, Hamilton Harbour Commissioners and private 
developers have adopted landfilling as a technique for 
acquiring new land for expansion for many years. Concern 
has been voiced that the landfilling has also become the 
object of public concern lest the water quality of the 
harbour be degraded. Both industry and the municipalities 
recognize the importance of the harbour water quality and 
have installed treatment facilities. While it is unlikely 
that the water quality in the harbour regions directly 
adjacent to the industrial shoreline can be brought to a 
level which meets the criteria for fish, wildlife and 
recreation, these criteria could be achieved in the remain- 
ing portions of the harbour. 

PURPOSE OF THE STUDY 

In 1972, the Ontario Ministry of the Environment undertook a 
water quality study of the harbour. The objectives of the 
study were: to determine existing water quality and its 
trends over recent years; to identify problems and recommend 
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corrective measures; and to determine the effects of land- 
filling on water quality. The water quality in the harbour 
is a result of the interactions of biological, chemical and 
physical processes in what was found to be a very dynamic 
environmental system. Thus, changes in any one process will 
affect the harbour water quality; consequently, it is neces- 
sary to understand the processes and their interactions in 
sufficient detail if intelligent decisions are to be made. 

OUTLINE OF STUDY 

Extensive water quality surveys were carried out during 1972 
and 1973. Various components of the harbour processes were 
examined using the survey data as a base. The study of each 
of these components formed the subject matter of a technical 
bulletin,. A draft of these reports was reviewed by Ministry 
personnel and subsequently amended incorporating their 
comments. 

Each of these bulletins appears as a separate section of 
this report. Some' duplication of figures and tables obviously 
results, but this is necessary if the sections are to be 
self-contained. The sections are as follows: 

A. Harbour Periodicities 

B. Nutrient Trend Analysis 

C. Iron-Phosphate Relationships 

D. Dissolved Oxygen 

E. Numerical Model 

F. Bottom Fauna and Bacteriology 
H. Survey Data 

In this preface, the results of the study are summarized and 
merged to provide comprehensive statements concerning the 
water quality as related to the purposes of the study. 

Several sources were used in the collection of the water 
quality data which form the basis of this report. Water and 
bottom sediment samples were obtained by conventional grab 
sampling techniques. Recording instrumentation (wind, water 
quality and current) was used to obtain time histories of 
physical properties and water quality. Bottom fauna samples 
and fish were collected. Historical water quality data were 
obtained from surveys conducted by the Hamilton Municipal 
Laboratories and the Ontario Water Resources Commission. 

SUMMARY AND FINDINGS 

Trophic State 

The hypertrophic state (mean annual chlorophyll a 20 ug/1 
and total phosphorus 4 to 200 ug/1) of the harbour fails to 
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meet the provincial objectives for fish and wildlife or 
aesthetics and recreation. Annual nitrogen and phosphorus 
loadings are five and thirty times the stock of these 
nutrients in the harbour respectively, indicating the 
importance of chemical and biological removal in preventing 
increases in their concentrations. An analysis of the iron- 
phosphorus relationships in the harbour indicates that 
phosphorus removal will not cause significant increases in 
the iron concentrations. Total phosphorus and ammonia 
concentrations in the epilimnion during stratification are 
significantly higher than those in the hypolimnion, indi- 
cating that destratif ication would be beneficial. 

Dissolved Oxygen 

Dissolved oxygen levels below 5 mg/1 exist in the lower 
portion of the hypolimnion for large regions of the harbour 
from July to September except during high wind conditions. 
These portions of the harbour do not meet the provincial 
criteria of greater than 5 mg/1 for the protection of fish 
and wildlife. The areas in the deep portion, near the 
Hamilton WPCP, along the industrial south shore and at the 
western end have the largest oxygen deficits. It is esti- 
mated that approximately 66 percent of the deficit is caused 
by discharges and sediment demand and 33 percent by phytoplankton 
decay in August. A reduction in chemical oxygen demand 
loadings of approximately 66 percent will result in a 
reduction in deficit by about 5 percent in the first year. 
The sediment demand should be controlled by nutrient concen- 
tration reductions, oxidation and possibly selected sediment 
removal. 

Landf illing 

Landfilling is causing a measurable degradation of water 
quality in the region of the filling. Regions of the 
harbour adjacent to the outfalls will experience higher 
concentrations although only a very small portion will 
exceed the 33 percent increase limit defined as an objective 
for fish and wildlife. This problem can be relieved by 
ceasing to fill, improving treatment at the principal 
discharges and/or utilizing more harbour water for assimi- 
lation through the use of destratification. 

Microbiology 

Bacterial levels of coliforms, fecal coliforms and entero- 
cocci are generally higher than permitted for recreational 
waters. If this use is contemplated the sources of con- 
tamination must be determined and corrected. The storm 
water sewers are a likely source, however, it was not 
possible to correlate rainfall with higher concentrations of 
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total phosphorus or nitrate during the 1972 surveys. (These 
fr^^o^'^i parameters are normally the most sensitive in urban 
areasand have been useful in other studies in identifying 
storm water overflows or bypasses.) It is considered that 
other sources are also contributing to bacterial contamin- 
ation like industrial discharges or wildlife. No changes in 
iDacterial levels was observed in the period 1970 to 1972. 

Bottom Fauna 

The surveys in 1972 and 1973 indicate a slight improvement 
tnr. ^^"^"^ fauna species in density and diversity from the 
1964-65 survey. This may be related to reductions in 
ammonia concentrations. However, the region devoid of 
macroinvertebrate species has shifted with the location of 
tne Hamilton sewage treatment plant discharge. The toxicity 
existing m this region should be identified and corrected. 

Heavy Metals in Sediments 

Compared to other Great Lakes regions, high concentrations 
ot lead, chromium, cadmium, mercury, and iron were found in 
the sediments, particularly, in the Randies Reef area 

u f^""?"" ^^^^ ' '^^^ sources of heavy metal discharges 
should be Identified and some action taken to reduce these 
loadings. 

Pestici des, Mercury and PCB Levels in Fish 

With the exception of one White Bass, mercury levels in fish 
puree were less than 0.5 ppm. Predator fish like White Bass 
had levels of PCB greater than 5 ppm and as high as 11.7 ppm 
^^o^^^ '"^u^"''^'-'^^^^'^ ^"^ unacceptable level. This indicates 
nnS 1 ^"^1 ^?^^ exist and should be identified and controlled. 
^^S3. -^^Y^-^^ ^^ fi^^ puree are lower than the breakdown products 
DDE and DDD indicating that recent DDT discharges are not a 
problem. ^ 

Trends 

Ammonia levels in the harbour have dropped from a mean value 
of 6.9 mg/1 in 1970 to 1.7 mg/1 in 1972 due to control 
measures. The bottom fauna macroinvertebrate community has 
improved slightly in density and diversity between 1964-65 
and 1972-73. Although chlorophyll a concentrations have 
remained approximately constant since 1966, phosphorus 
concentrations have increased slightly since 1968. Oxygen 
demands due to bottom sediments have increased over the 
years and the dissolved oxygen deficits in the hypolimnion 
are presently a problem. 
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RECOMMENDATI ONS 

The following recommendations are based on the premise that 
much of the harbour waters can be brought to a water quality 
standard that will be suitable for fish, wildlife and 
recreation. These recommendations are arranged in order of 
priority to achieve the required water quality enhancement. 

The hypertrophic state and dissolved oxygen deficits are the 
most pressing water quality problems. If the harbour is to 
support fish and aquatic life or meet the provincial standards 
for aesthetics and recreation, these problems must be rectified 
As the problem will probably require 10 to 20 years to 
correct, it is important that these steps be undertaken as 
soon as possible. It is recommended that: 

1. The concentrations of nutrients in the harbour be 
reduced by treatment. 

2. The principal chemical oxygen demand loadings be 
reduced by approximately 6 percent. 

3. Artificial destratif ication be considered as a means of 
reducing nutrient levels in the photic zone, oxidizing 
bottom sediments and preventing phosphorus from recycling. 

4. Aeration of the hypolimnion and selective sediment 
removal be investigated as a means of reducing the 
sediment oxygen demand. 

If the harbour is to support fish and wildlife, it is 
important that the causes of toxicity to bottom fauna in the 
southeast portion of the harbour and the sources of heavy 
metal and PCB contamination be identified and rectified. 

The present levels of bacterial contamination do not meet 
the water quality criteria for recreational use. If this 
use is contemplated for the harbour, the sources of bacterial 
contamination must be identified and rectified. 

Landfilling is causing a local degradation of water quality. 
While the regions of the harbour affected by a degradation 
which violates the provincial criteria for fish and wildlife 
are small, it is advisable to curtail landfilling with the 
existing water quality problems in the harbour. 

FURTHER STUDIES 

Studies will be undertaken in 1974 to further amplify the 
findings of this report. As it is only possible to validate 
a two-dimensional numerical model in a dynamic system 
qualitatively, a three-dimensional model will be developed 
which will be validated by operating concurrently six 
installations, each consisting of one recording current 
meter and one water quality meter. The time histories of 
the major industrial outfalls will also be measured and 
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incorporated into the model and a non-conservative will be 
modelled. The heavy metal sediment results will be analyti- 
cally treated to relate sediment samples. As the oxygen 
demand of the sediments is important, sediment demands will 
be determined both in the laboratory and in the field at key 
locations. Also, additional surveys will be undertaken to 
identify bacti sources. ■ 
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HARBOUR 
STUDY 74 



section I 

Proposed 1975 
water quality surveys 



1975 HAMILTON HARBOUR SURVEY PLAN 

Hamilton Harbour is to be artificially mixed by aeration to relieve hypo- 
limnetic annoxia. A diffuser will be placed at the deepest point of the 
central basin to ensure complete vertical mixing of the water column. 
Combined with wind mixing and horizontal transport, it is likely that 
an oxidized microzone will be established and maintained at the sediment- 
water interface. 

In conjunction with this project, a regular water and sediment quality 
survey of Hamilton Harbour during 1975 is desired. The purposes of 
this survey are as follows: 

T . To determine the extent of mixing in the harbour. 

2. . To examine the effect of mixing on the bottom sediments, i.e. , to 

determine the presence of an aerobic microzone. 

3. To determine changes in the composition, abundance, and dist- 
ribution of fish and plankton. 

4. To determine whether the epilimnetic nutrient concentration has been 
reduced . 

5. To determine the effect of mixing on dissolved oxygen stock. 

(a) WATER QUALITY 

A comprehensive water quality survey is to be undertaken with sampling 
performed at least biweekly from early April until commencement of mixing, 
and then weekly during the period of mixing, and biweekly from cessation 
of mixing until the end of October. Because of Great Lakes commitments, 
samples for bacti should be obtained on Tuesdays. The following stations 
and depths are to be employed . 



Station 


Latitude 


Longitude 


4 


43-17-11 


79-48-00 


252 


43-18-20 


79-48-40 


258 


43-17-20 


79-50-10 


270 


43-16-50 


79-52-20 



Depths to be Sampled (m) 

0.2, 3.0, 6.0, 12.0, 18.0 

0.2, 3.0, 6.0 

0.2, 3.0, 6.0, 12.0, 20.0 

0.2, 3.0, 6.0, 12.0 
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Locations of the sampling stations and the propsed diffuser location are shown 
in Figure 1 . Station ^ was used in the 1967 survey and is in the deep hole immed 
lately west of station 251. Station 270 is immediately northeast of station 265 and 
is selected for a greater available depth. Prior to sampling, the presence of an 
adequate water depth should be ascertained, and if such is not present, the samp 
ling vessel should move to a location which provides sufficient depth. 

Depth profiles of temperature and dissolved oxygen are to be obtained weekly 
at each station. 

The following chemical parameters are to be determined on all samples: 

Conductivity 

pH 

Turbidity 

Ammonia 

Total Kjeldahl Nitrogen 

Nitrite 

Nitrate 

Total Phosphorus 

Soluble Phosphorus 

Dissolved Reactive Silicates 

Alkalinity 

Iron 

Chlorophyll 

The following parameters are to be determined on the surface and bottom 
samples at each station: 

COD 

BOD 

TOC 

Total Conforms 

Fecal Conforms 

Enterococci 

Sulfur Cycle Bacteria 

Nitrifiers 

Total heterotrophic Bacteria 

Pseudomonads 

Sodium (monthly only) 

Potassium (monthly only) 

Calcium (monthly only) 

Magnesium (monthly only) 

Chloride (monthly only) 

Sulphate (monthly only) 
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Total and soluble fractions of the following heavy metals are to be done monthly. 
Two depth composites for each station, one filtered in the field, are to be collected, 
using equal volumes of water from each of the above indicated depths: 

Chromium 

Manganese 

Cobalt 

Nickel 

Copper 

Zinc 

Cadmium 

Mercury 

Molybdenum 

Lead 

Pesticides and PCB's are also to be determined on similar monthly depth 
composite samples . 

For purposes of nutrient and dissolved oxygen trend studies , monthly samples 
will be collected by the Lake Systems group on calm days. At the 17 station 
grid as used in 1972, temperature and dissolved oxygen-depth profiles are to 
be taken. Bathythermograpth casts and dissolved oxygen meters will be used 
if possible, to simplify data collections. Nutrient samples will be taken at the 
surface and bottom from the 8 station subset used intensively in 1972. 

(b) SEDIMENT CHEMISTRY 

Sediment samples are to be taken by dredge once before mixing is initiated and 
at least once after mixing has started . In order to locate areas of greatest oxygen 
demand and discern areas of related sediment quality, a set of 20 stations should 
be sampled: (Figure 2) 

Stations: 4,15,20,252,253,254,255.256,257,258,259,260,261.262,263, 

264,268,269,270 and 1030. 

Stations 4 and 270 are as above and replace stations 251 and 265 respectively. 
The following parameters should be determined on all cores: 

BOD 

COD 

TOC 

Loss on ignition 

Total Kjeldahl N 

Total Phosphorus 

Acid-extractable Phosphorus 
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PCB 

Chromium 

Manganese 

Iron 

Cobalt 

Nickel 

Copper 

Zinc 

Cadmium 

Mercury 

Molybdenum 

Lead 

An attempt should be made to measure redox potential of the sediment as a 
function of sediment depth up to 20 cm below the sediment-water interface. 
Temperature, dissolved oxygen and conductivity of the water immediately 
above the sediment should be measured if possible. 

(c) PHYTOPLANKTON AND ZOOPLANKTON 

A biological survey by a University contract will compliment the Hamilton 
Harbour Study 1975. Phytoplankton enumeration and identification, primary 
productivity, and bioassays will form the bulk of the biological work. The 
four stations will be sampled at three depths to correspond with the water 
quality survey. 

Zooplankton work should be simple, consisting of a vertical net haul through 
the entire water column at all stations. Relative species abundance should be 
determined, and total population size represented as dry weight. Diurnal 
migrations in the vertical and horizontal plane result in difficulties in deter- 
mining spatial distribution of zooplankton; thus more emphasis should be 
placed on seasonal changes within the period of observation . 

Zooplankton grazing experiments and/or analyses of gut contents to deter- 
mine se!ective feeding on the phytoplankton is important. 

The attached sheet outlines the biological survey to be done by contract. 
Station composites should be collected by M.O.E. (100 ml sample) and pre- 
served in Lugol's Iodine. Cell measurements, enumeration and identification 
will be done as a parallel study to determine phytoplankton biomass and species 
diversity. 
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id) BOTTOM FAUNA 

Bottom fauna (composition and abundance) should be collected using an 
Ekman bottom grab sampler. Collections should be made at the same time 
sediment samples are collected for chemical analysis. The should be at the 
20 locations used for the sediment analysis, for comparison with previous 
studies and provide estimates of the effect of bottom fauna on sediment oxygen 
uptake. 

Invertebrates should be separated from the sediment by washing through a 
sieve with 0.65 mm apertures, and preserved in ethanol. 

As background data is available, increased sampling frequency would be 
of value to determine community structure and diversity during 1975. This 
however, is not essential, but would help determine factors affecting the 
bottom fauna composition and abundance. 

(e) FISH SURVEY 

Perhaps one of the first biological responses to aeration will be changes in 
the relative abundance, distribution and composition of the fish populations 
of Hamilton Harbour. At least two nettings, before and after mixing, are 
required at the four sampling stations. Some fish will be collected for flesh 
analyses (heavy metals, pesticides, P.C.B.). It Is hoped that the Ministry 
of Natural Resources will do the nettings. 
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BIOLOGICAL SUKVEY 
HAMILTON HARBOUK 197 5 



Biological responses to mi.<lrip; are to bo dGtermlned bv 
sarapllne at weekly intervals at fo^ur stations. Volume 
samples will be collected ot three depths, representative vV 
J^'^ !°P' middle and bottom of the water column. Phytoplank- 
'^'^i'^ Identification and enuraeratloii will be determliiefl for 
all samples. Chlorophyll = a will also be determined for 
each sample. 

Temperature and dissolved o^xyecn profiler, will be measured 
at the cime of sampling for each station. Vertical net 
hauls through the Mater column using a 60 pm mesh will be 
used to estimate zooplankton crops. 

Primary productivity experiments usinp the light and dark 
bottle technique will be accompanied by nieaGurements of the 
vertical attenuation of solar radiation to form an essential 
component of this study. Combined with simple bloat^say 
studies the production potential of the harbour can be 
determined. 
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FIG. I PROPOSED WATER QUALITY SAMPLING STATIONS, 1975 
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FIG 2 PROPOSED SEDIMENT SAMPLING STATIONS, 1975. 
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